Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,206,1070CA25.3% 165.8 / 12.8 83A308A (GCG→GCTmetCCystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base C (22/40);  new base A (21/0);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 4.82e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGC  >  ACB143/5205972‑5206243
                                                                                                                                       |                                                                                                                                        
gCTCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAATTTAATCg                                                                                                                                         >  5:66911/1‑137 (MQ=255)
  tCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGTCGCGGTCGCCTTTTACTTTCATCGcc                                                                                                                                       <  5:204129/137‑1 (MQ=255)
  tCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCGTTTTACTTTCATCGcc                                                                                                                                       >  3:280579/1‑137 (MQ=255)
  tCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCGGTTTACTTTCATAGcc                                                                                                                                       >  1:92102/1‑137 (MQ=255)
  tCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGcc                                                                                                                                       >  6:277132/1‑137 (MQ=255)
  tCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGcc                                                                                                                                       >  8:362424/1‑137 (MQ=255)
    aGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGccac                                                                                                                                     <  7:132318/137‑1 (MQ=255)
       gAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCCCCCCTCCCAGCGAGTCGTCCAGCACAAAAATCTGCAGCGCCTCCACCCGGCGGTTGAGCGCGGCGCGGTGGCCGTGTAGATTAAACGCcaccac                                                                                                                                  >  1:354758/1‑137 (MQ=255)
        aGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAATTTCATCGCCaccacc                                                                                                                                 >  2:6500/1‑137 (MQ=255)
         gCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCAccg                                                                                                                                <  7:297151/137‑1 (MQ=255)
            gTGCGACATGCTCGATGAGTGTTTGTTCATTCTCTCAAACCCGCCCAGCGAGTCGGACAGCAAAAAAATCTGCAGCGCCTACACCTGTCGGTTGATCGCGGCGCGGTCGCCTTTTACTTTCATCGACACCAcagccc                                                                                                                             <  6:101041/137‑1 (MQ=255)
            gTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATAGCCACCAccgccc                                                                                                                             >  2:261533/1‑137 (MQ=255)
               cGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCAccgccccgc                                                                                                                          >  4:252992/1‑137 (MQ=255)
                 aCATGCTCCAGGAGTGGTTGATCTTGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCa                                                                                                                        <  7:167646/137‑1 (MQ=255)
                 aCATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCa                                                                                                                        <  1:48232/137‑1 (MQ=255)
                 aCATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCa                                                                                                                        <  3:323309/137‑1 (MQ=255)
                        ccAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCa                                                                                                                 <  6:126640/137‑1 (MQ=255)
                          aGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCCCCACAGCGCATc                                                                                                               >  2:72829/1‑137 (MQ=255)
                           ggAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCt                                                                                                              <  6:15532/137‑1 (MQ=255)
                               tGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAATTTAATAGCCCCCCCCGCCCCCCCCCTGCGCCTCTGGCg                                                                                                          >  2:174156/1‑137 (MQ=255)
                                ggTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGt                                                                                                         <  8:269267/137‑1 (MQ=255)
                                 gTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGtt                                                                                                        <  2:19670/137‑1 (MQ=255)
                                   tGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTg                                                                                                      <  7:122013/137‑1 (MQ=255)
                                   tGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTg                                                                                                      <  5:359665/137‑1 (MQ=255)
                                   tGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTg                                                                                                      <  8:180452/137‑1 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGTCGGTTGAGCGCGGCGCGGTCGCCGTTTAATTTGATAGCCACCCACCCCCCCCCCCAGCGCCTATGGCgggttg                                                                                                     >  1:279537/1‑132 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCTCCATTTACTATCATAGCCACCACCCCCCCGCCCAAGCGCAACAGGGGTTTgg                                                                                                     >  2:314628/1‑136 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCGTTTTACTTTAATAGCCACCACCGCCCCGCCACTGCGCCTCTGGCGTTTgg                                                                                                     >  3:158980/1‑136 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTTCTTTAAAAGCCACCACCGCCCCCCCACCGCGCAACAGGCGGTTGc                                                                                                     >  7:98512/1‑137 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAGTTTCATCGCCACCCCCGCCCCCCCACTGCGCATATGGGGTTTgg                                                                                                     >  4:20036/1‑136 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGc                                                                                                     >  8:353848/1‑137 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATAGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGc                                                                                                     >  5:133189/1‑137 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAATTTAATCGCCCCCACCCCCCCCCCACAGCGCAACTGGCGgtttg                                                                                                     >  5:319402/1‑135 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTTCTTTCAAAGCCACCACCCCCCCGCCACAGCGCATCTGGCGTTTgg                                                                                                     >  7:153387/1‑136 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTTATTTCATAGCCACCACCCCCCCGCCACGGCCGATCTGGCGGTTgg                                                                                                     >  2:294865/1‑136 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTACTTTCATCGCCCCCACCCCCCCGCCACTGCGCAACTGGCGTTTGc                                                                                                     >  2:232617/1‑137 (MQ=255)
                                    gATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTAATTTCATAGCCACCACCGCCCCGCCACTGCGCATCTGCCGTTTGc                                                                                                     >  5:139398/1‑137 (MQ=255)
                                         tGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACAGCGCATCTGGCGTTTGCACAGt                                                                                                >  4:135244/1‑137 (MQ=255)
                                              tcCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCGTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCAtt                                                                                           >  5:252981/1‑136 (MQ=255)
                                                 aCCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTACTTTCATCGCCACCACCGCCCCGCCCCTGCGCCTCTGGCGGTTGCACAGTTGAtgggg                                                                                        >  5:94063/1‑134 (MQ=255)
                                                   cccGCCCAGCGAGTCGTCCAGCACAAAAATCTGCAGCGCCTCCCCCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTTCTTTCATCGACACCACCGCCCCGCCACAGCGCATCTGGCGTTTGCCAAGTTGCTGCTGCg                                                                                      >  5:263440/1‑137 (MQ=255)
                                                   cccGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCg                                                                                      >  3:174022/1‑137 (MQ=255)
                                                   cccGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTAATAGCCACCACCGCCCCGCCACTGCGCATATGGCGTTTGCCAAGTTCATGCTGCg                                                                                      >  8:467/1‑137 (MQ=255)
                                                     cGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCggg                                                                                    <  2:160875/137‑1 (MQ=255)
                                                       cccAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTTATTTCATAGCCACCACCGCCCCCCCACAGCGCATCTGGCGGTTGCACAGTTGATGCTGCGGGTg                                                                                  >  5:336988/1‑137 (MQ=255)
                                                        ccAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTgg                                                                                 >  4:126726/1‑137 (MQ=255)
                                                        ccAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTAATTTGATAGCCACCACCGCCCCGCCACAGCGCATCTGGCGTTTGCACAGTTGATGCTGCGGGTgg                                                                                 >  3:168342/1‑137 (MQ=255)
                                                          aGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCt                                                                               <  1:283838/137‑1 (MQ=255)
                                                           gCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTTATTTAATAGCCCCCACCGCCCCGCCCCTGCGCATCTGGCGTTTGGACAGTTGATTCTGCGGGTggggt                                                                              >  3:286351/1‑134 (MQ=255)
                                                             gAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTAGTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTGATGCTGCGGGGGGCtttg                                                                            >  6:333709/1‑136 (MQ=255)
                                                               gTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCa                                                                          <  7:194217/137‑1 (MQ=255)
                                                               gTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCa                                                                          <  7:362424/137‑1 (MQ=255)
                                                               gTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCa                                                                          <  1:314628/137‑1 (MQ=255)
                                                               gTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCa                                                                          <  1:302606/137‑1 (MQ=255)
                                                               gTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTTCTTTAATAGCCACCACCGCCCCGCCCCTGCGCATCTGGCGTTTGCACAGTTACTGCTGCGGGTGGCTTTGCa                                                                          >  4:266052/1‑137 (MQ=255)
                                                                tCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAg                                                                         <  4:174022/137‑1 (MQ=255)
                                                                tCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATAGCCACCACCGCCCCGCCACTGCGCATATGGCCTTTGCACAGTTCATGAAGCGGGTGGCTTTcccg                                                                         >  8:282355/1‑135 (MQ=255)
                                                                 cGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGc                                                                        <  4:71348/137‑1 (MQ=255)
                                                                  ggCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGcc                                                                       <  4:28372/137‑1 (MQ=255)
                                                                     cAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATAGCCACCACCGCCCCGCCACTGCGCCTCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGccccc                                                                    >  5:345034/1‑136 (MQ=255)
                                                                           aaaaaTCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAAt                                                              <  1:174156/137‑1 (MQ=255)
                                                                           aaaaaTCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAAt                                                              <  3:219498/137‑1 (MQ=255)
                                                                                 cTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCt                                                        <  6:68348/137‑1 (MQ=255)
                                                                                 cTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGAGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCt                                                        <  1:80838/137‑1 (MQ=255)
                                                                                      gcgcCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCGTTTAATTTCATAGCCACCACCGCCCCGCCACTGCGCATATGGCGTTTGCACAGTTCATTCTGCGGGTGGCTTTACAGCCCCGGGTAAAACACCTGCTCg                                                   >  7:101008/1‑137 (MQ=255)
                                                                                          cTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCt                                               <  7:353848/137‑1 (MQ=255)
                                                                                            ccaccaGGCGGTTGAGCGCGGCGCGGTCGCGGTTTTATTTGATAGCCACCACCCCCCCGCCCCCGCGCCTCTGGGGTTTTCACAGTTAATTCTGCGGGTGGCGTTGCAGCCCCGGGGTATACACCTGCTCGAccagg                                             >  6:46130/1‑134 (MQ=255)
                                                                                              accaGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTgcgc                                           <  1:232617/137‑1 (MQ=255)
                                                                                              accaGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTgcgc                                           <  1:336540/137‑1 (MQ=255)
                                                                                              accaGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTgcgc                                           <  3:252992/137‑1 (MQ=255)
                                                                                                caGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGt                                         <  8:153387/137‑1 (MQ=255)
                                                                                                    cGGTTGAGCGCGGCGCGGTCGCCGTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGAcc                                     >  3:309471/1‑137 (MQ=255)
                                                                                                      gTTGAGCGCGGCGCGGTCGCGTTTTAAATTAATAGCCCCCACCGCCCCCCCACTGCCCCTCTGGCGTTTTCCCACGTGATGCTGCGGGTGGCTTTACAGCCCCGGGTAATACACCAGCTCGACCAGCGCGGGAccgg                                   >  2:158268/1‑135 (MQ=255)
                                                                                                         gAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCa                                >  7:26231/1‑137 (MQ=255)
                                                                                                             gcggcgcggTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGaa                            <  4:158980/137‑1 (MQ=255)
                                                                                                                gcgcggTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAgcg                         <  6:168723/137‑1 (MQ=255)
                                                                                                                 cgcggTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAgcgc                        <  7:467/137‑1 (MQ=255)
                                                                                                                 cgcggTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAgcgc                        <  6:336988/137‑1 (MQ=255)
                                                                                                                    ggTCGCGTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGGGTTCCAGCCCCGGGGAATAAACCTGCTCGCCCTGCGCGTGACCTTCCAGGAAGCGCGcc                     >  3:66686/1‑137 (MQ=255)
                                                                                                                              tACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAggg           <  4:165695/137‑1 (MQ=255)
                                                                                                                                cTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCTCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCa         <  3:313628/137‑1 (MQ=255)
                                                                                                                                cTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCa         <  1:301767/137‑1 (MQ=255)
                                                                                                                                   tCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTg      <  1:104450/137‑1 (MQ=255)
                                                                                                                                       cGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGgcgc  >  6:71889/1‑137 (MQ=255)
                                                                                                                                       cGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCAACCTGCATGGCATTGgctg  >  8:343557/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GCTCAGCGAGCAGTGCGACATGCTCCAGGAGTGGTTGATCATGCTCTCCACCCCGCCCAGCGAGTCGGCCAGCACAAAAATCTGCAGCGCCTCCACCAGGCGGTTGAGCGCGGCGCGGTCGCCTTTTACTTTCATCGCCACCACCGCCCCGCCACTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGC  >  ACB143/5205972‑5206243

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: