Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1613,5970GC24.5% 350.8 / 44.5 247intergenic (+458/+565)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
Reads supporting (aligned to +/- strand):  ref base G (61/93);  new base C (1/49);  total (69/178)
Fisher's exact test for biased strand distribution p-value = 2.76e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATAC  >  DH1/613448‑613740
                                                                                                                                                       |                                                                                                                                               
cccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                                <  4:362678/151‑1 (MQ=255)
 ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                                >  7:123747/1‑150 (MQ=255)
 ccGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                                                                                                                                                >  1:883501/1‑150 (MQ=255)
  cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTCCTTTTGGAGGGGCAGaa                                                                                                                                               >  1:757427/1‑150 (MQ=255)
  cGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCAGCTTCGGCGGGTTTTTGCTTTTGTAGGGGAAgta                                                                                                                                               >  1:1579962/1‑148 (MQ=255)
   gTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                                                                                                                                             <  4:210363/151‑1 (MQ=255)
    gaCTGAACCCGGTTAGGCGGGTTGAGGCAGTGACCGATTGATGGCCTAGACCCAAGACTCTACCCCACAGAAACCGGCTTATCGGTCAGTTTCACCTGTTTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGa                                                                                                                                            <  3:1530967/150‑1 (MQ=14)
     aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATAGGTCAGTTTCACCTGATTTACGTAAAAACCCGCGTCTGAGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                               >  3:1466084/1‑147 (MQ=255)
     aCTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGATTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGaa                                                                                                                                               >  2:1150332/1‑147 (MQ=255)
       tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCATCGGCGCCGTTTTGCGTTTGGAGGGTGAGATAGatg                                                                                                                                          >  3:1253643/1‑150 (MQ=37)
       tGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAAACCCCTTATAAGTAAGTTTCACCTTATTTAAGTAAAAAACAGATTCGGCATGTTTTTGATTTTGGAGGGGCAGAAAGatt                                                                                                                                          >  2:1779602/1‑149 (MQ=25)
        gAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                        <  4:594663/151‑1 (MQ=255)
         aaCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATGGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaa                                                                                                                                        <  1:529516/150‑1 (MQ=255)
          aCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTAACCTGATTTACGTAAAAACCCGCTTCGGAGGGTTATTGCTTTTGGAGGGGCAGAAAGatgaatg                                                                                                                                      >  7:839904/1‑151 (MQ=255)
          aCCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCACGCATATCCGTCAGTTTCACCTAATTTAAGTAAAAAACAGCTTCTCACGGTTTTAGCATTTTGAGTAAAATAAAGatgaat                                                                                                                                       >  6:254560/1‑150 (MQ=255)
               ggTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGt                                                                                                                                 <  7:1079649/151‑1 (MQ=255)
                gTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTc                                                                                                                                <  2:332295/151‑1 (MQ=255)
                  aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga                                                                                                                                     <  4:901452/144‑1 (MQ=255)
                  aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatgaatga                                                                                                                                     >  3:899400/1‑144 (MQ=255)
                  aGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                              <  1:212875/151‑1 (MQ=255)
                   ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCac                                                                                                                             <  2:1328143/151‑1 (MQ=255)
                   ggCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCa                                                                                                                              <  2:704109/150‑1 (MQ=255)
                      ggctTGAGCCAGTGAGCGATTGCTGGGCTAGATGAATGACTGTCCACGACAGAACCAGGCTTATCGGTCAGTTTCACCTGATTTACGTAACAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          <  4:1214785/150‑1 (MQ=255)
                       gctTGAGCCTGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACATAAGCCTGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                         <  2:53207/151‑1 (MQ=255)
                       gctTGAGCCAGTGAGCGATTGCTGGGCTAGATGAATGACTTTCCACTACAGAACCAGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgac                                                                                                                          <  8:617372/150‑1 (MQ=255)
                       gctTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacga                                                                                                                           <  7:869881/149‑1 (MQ=255)
                       gctTGAGCCAGTGAGCGATTGCTGCCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGaggaatgccatgcacacacg                                                                                                                         >  1:823822/1‑138 (MQ=21)
                          tGAGCCATTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCt                                                                                                                       <  1:987740/150‑1 (MQ=255)
                           gcgCCAGGGAGCGATGGCGGGCCTAGATCAATGACTGCCCCCGACAGAACCCGGCTTAGGGGTCAGTTTCCCCTGATTTACGTAAAAACCCGCTTCGGCGGTTTTTTGCTTTTGGAGGGGCACAACTATGAATGACTGTCCACGACGCta                                                                                                                      <  2:1212046/148‑1 (MQ=18)
                            aGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCATTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCtata                                                                                                                    <  3:1308921/151‑1 (MQ=255)
                            aGCCAGTGAGCGATGGCTGGCCTCGAGGAATGACGGCCCACGACACAACCCGGCTAATCGCGCAGTTGCCCCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTTGATGGGCAGAAAGATGAATGACTGTCCACGACGCtata                                                                                                                    <  2:991298/151‑1 (MQ=18)
                             gCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGCTTATCGGTCAGTTTCACCTGATTTACCTAAAAACCCCCATCGGAGGGGATGACCTATGAGAGCGGCTCAAAGAACAAAGACTGTGCACCACGCTAGAc                                                                                                                   >  1:1825248/1‑151 (MQ=19)
                               cAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAAACCGGCTTATCCGTCAGTTTCACCTGATTTACGTAAAAACCAGCTTCGGAGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCaagagcggatcccc                                                                                                                 >  7:1594891/1‑141 (MQ=18)
                                  tGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACACGGCTAATAGCTCACTTTCCCCAGATTTAACTAAAAACCAGCTTACGCAGGATTTTGGTTTTGGAGGGGCAGAAAGATGCATCACTGTCAACGCCGCTATACCCaa                                                                                                               >  7:1485446/1‑150 (MQ=255)
                                  ggAGCGATTGCTGGCCTAGATTAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCaaa                                                                                                              <  1:1069445/150‑1 (MQ=255)
                                   ttgcttttGCTGGCCTAGATTAATGACTGTCCCGCACAGAACCCGGCTTATCGGTCTGTTGCACGTGATTTACGTAAAACCCCGGTTCGGCGTGTTTTTGCTTTTGGCGGGGCATAAAGATGATTGAGTGTCCACGACGCTATACCCaaa                                                                                                              <  3:85815/144‑1 (MQ=255)
                                    cgCGATTGCTGGCGTAGCTGAATGACTGGCCCCCACAGAACCCGGCTTAGCGGTCAGTTTCACGTGATTTACGTAAAACCCCGGTTCGGGGGGTTTTTGCTTTTGCAGGGGAACAAAGATTAATGCCTGTCCACGACGCTATACCCAaaag                                                                                                            <  4:1635752/150‑1 (MQ=255)
                                      ggtttGCTGGCATATTTTAATGACTGGCCACGACAGAACCCGGGTTATCGGTCAGTTTCACCTTATTTACGTAAAAACCCGCTTCGGCGGTTTTTTGCTTTTGGTGGGGCAGAAACTTCAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                                          <  7:1685646/148‑1 (MQ=17)
                                      cGATTGCGGGCCTAGATTAATGACTGGCCACCACAGAACCCGGCTTATCGGTCAGTTGCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAaaagaa                                                                                                          <  1:1408986/151‑1 (MQ=255)
                                       gATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGGAGAAAGATGAATGACTGTCCACGACTCTATACCCAaaagaa                                                                                                          >  7:903152/1‑150 (MQ=255)
                                        aTTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                         >  6:1698161/1‑134 (MQ=255)
                                        aTTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCacgacg                                                                                                                         <  5:1697215/134‑1 (MQ=255)
                                         ttgggggcgtattttaatgaCTGCCCACCCCACAACCCGGGTTTTCGGTCTTTTTCACGTTATTTACTTAAAAAGGCGCTGGGGGGGTTTTTTGTTTTTGTTGGGGCACAAAGTTTAATGAGTGTCCACGACGCTATACCCAaaagaaag                                                                                                        <  5:93645/135‑1 (MQ=255)
                                         ttGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCATAAAAATCAATGACTGTCCACGAAGCTATACCCAaaataaatc                                                                                                       >  2:743048/1‑149 (MQ=18)
                                          tGCTGGCCTAGATGAATGACGGGCCACGACACAACCCGGCTTATCGGTCAGTTTCACCTTTTTTACGTAAAACCCCGCTGCGGGGGGTTTTTGCTTTTGGAGGGGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGc                                                                                                       <  5:1296358/150‑1 (MQ=17)
                                             gggCCTAGATGAATCACTTTCCACGACAGAACCCGGCTTAGCGGTCTGTTTCACCTGATTTACGTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGGTGGGGCAGAAAGATTAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                   <  6:1317738/150‑1 (MQ=17)
                                              ggCCTTGATGAATGACTGCCCACGACATAACCCGGCTTATCGGTCAGTTTCACGTTATTTAGTTAAAACCCCGCTTGGGGGGTTTTTTGGTTTTGGAGGGGCAGAAAGATGATTGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCt                                                                                                   <  3:259852/150‑1 (MQ=17)
                                                catatatgaatgaCTGTCCACCACTTAACCGGGCTTCTGGGTCTGTTTCACCTGATTTACTTAAAAACCCACTTGGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTAt                                                                                                <  4:749695/146‑1 (MQ=14)
                                                 agaGATTAATGACTGGCCACCACAGAACCCGGCTTATCGGGAAGTTTCACCTGATTTACTTAAACACCCGCTTCGGCGGGTTTTTGGTTTTGGTGGGGCAGAAAGATGAATCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATc                                                                                               <  3:1671934/149‑1 (MQ=17)
                                                   aGATGAATGACTGTCCACGACACAAAACCCCTTATAAGTCATTTTAACCTGATTTACGTAAAAAACCGCTTCGCACGTTTTTTACTTTTGGAGGGGCAGAAAGATGAATGAATCTCAACGACACTATACCCAaaagaaaagcgcggaagcg                                                                                             >  8:479146/1‑139 (MQ=9)
                                                    catcactgaCTGCCCACGACAGAACCCGGCTTTTCGGTCAGTTTCACCTGATTTACGTAAAACCCCGCTTCGGCGGTTTTTTGCTTTTGTAGGGGCAGAAAGATCAACCACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCgg                                                                                             <  4:1413973/144‑1 (MQ=17)
                                                     attaatgaCTGTCCACTACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTAGTTAAAAACCCGCTTCGGCGGTTTTTTGCTTTTGTAGGGTCATAAAGATTAATTACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTc                                                                                           <  4:17517/148‑1 (MQ=17)
                                                       gaatgaCTGCCCACCACACAACCCGGCTTATCGGTCAGTTTCACCTGATTTACTTAAAAACCCTCTTCGGCGGGTTTTTGCTTTTGGTGGGGCCCAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCa                                                                                          <  2:1637596/150‑1 (MQ=14)
                                                       aaattcgtttaccaaaCATAACCCGTCTTATCGGTAATTTTAACCTGTTTTAGTTAAAAACCCGATTTTGCGGGTTTTTGTTTTTGTTGGGGCTTAAATATGATTGACTTTCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                         <  8:841009/136‑1 (MQ=2)
                                                        aatcaCTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAg                                                                                         <  6:913444/150‑1 (MQ=255)
                                                           gtgtatcgcaaaCAGAACCAGTCTTTTGGTGCAATTTCACGTGCTTAAGGTAAAACGGCTCTGCGGCGTGTTTTTGGTTTTGATGGTGGCCAAATTTGAATGACTTGCCACGACACTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                      <  1:1547683/139‑1 (MQ=2)
                                                           gaggggccaccaCAGAACCCGGCTTAGCGGTAATTTGCACCTGATTAACGTAAAACCCCGCTTCGGCGGGTTTTTGATTTTGGTGGGGCACAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                      <  3:573614/144‑1 (MQ=14)
                                                           gacgggccaccaCAGAACCCGGCTTAGCGGTCAGTTTCACCTGATTTACGTAAAACCGCGCTTCGGCGGGTTTTTGTTTTTGCAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGttt                                                                                      <  1:1278579/144‑1 (MQ=14)
                                                            aCTGTCCACGACTTAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAACCCGGCTTCGGCGGGTTTTTGTTTTTGTAGGGGCTCATATATGATTGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCa                                                                                    <  6:1787363/151‑1 (MQ=14)
                                                              cctcctttaCAGAACCAGTCTTATCTGTCAGTTTCACTGTCTTTAGTTAAACCCGGGCTGCGGCGGGTTTTTGGTTTTGTTGGGTCCGAAATATGATTTACTGCCCACGACGCTATACACAAAAGAAAGCGGCTTATCGGTCAGTTTCaac                                                                                  <  8:1279177/143‑3 (MQ=2)
                                                               tcccaccaCACAACCCGGCTTATCGCTCAGTTTCACGTTATTTACTTAAAAACCCTCTTCGGGGGGTTTTTGTTTTTGGCGGGGCTTAAATATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCAcc                                                                                  <  4:1832011/148‑1 (MQ=14)
                                                               cgccACTACAGAACCCGGCTTAGGGGCCAGTTGCACGTCATTTCGTTAAAACCCGGCTGCGGCGGGTTTTTGCTTTTGCCGGCGCACAAAGATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCt                                                                                 <  2:454960/149‑1 (MQ=14)
                                                                   ccGACACAACCCGGCTTATCGGTCATTTTCCCCTGATTTAGTTAAAAACCGGCTTCGGCGGGTTTTTGTTTTTGGTGGGGCTCAACTATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGtt                                                                             <  2:1490641/150‑1 (MQ=14)
                                                                    caacacAACCCGGCTTATGGGCAATTTGCCCGTCATTTACTTAAAACCGGTCTGCGGCGGGTTTTTGTTTTTGTTGGGGCACAATTATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  2:603763/149‑1 (MQ=9)
                                                                     gACAGAACCCGGATTATGGGTAATTTTCACCTGATTTACTTAAAAACCGCCTGCGGCGGGTTTTTGCTTTTGTAGGTGCCCAAAGATGAATGACTGTCCACGACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTa                                                                           <  1:730214/151‑1 (MQ=14)
                                                                      aCAGAACCCGGCTTATCGGTCATTTTCACGTGATTTAGGTAAAAACCCGCTGCGGCGGGTTTTTGGTTTTGGCGGGGCATAAAGATGAATGACTCTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  7:765919/149‑1 (MQ=14)
                                                                      aCAGAACCCGGCTTATCGGTCATTTGCACCTGATTTAGTTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGGCGGGGCCCAAACATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGttt                                                                            <  4:1663201/149‑1 (MQ=14)
                                                                        aGAACCCGGCTTATCGGTCAGTTTCTCCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGaa                                                                                                                                               <  3:1075384/80‑1 (MQ=255)
                                                                         caaCCCGGCTTATCGGGCAGTTTCACCTGATTTACGTAAAAACCGGCTGCGGCGGGTTTTTGCTTTTGGAGGGGCACAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTa                                                                       <  3:737756/150‑1 (MQ=255)
                                                                          aaCCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTACTTTTGGAGGGGCAGAAAGat                                                                                                                                           >  5:776306/1‑82 (MQ=21)
                                                                             ccGGATTATGGGTCAGTTTCCCGTGATTTACTAAAAAACCCCCTTCGGGGGGTTTTTGGTTTTGCTGGCGCAGAAAGATGAATGACTTTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTaaaaa                                                                   <  7:587333/151‑1 (MQ=11)
                                                                               ggctgatggctaactttCAGCTGATTTCGGTAACAACCGGCTGCGGCGGGTTTTTGCTTTTGTCGGGGCACAACTATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAcc                                                                 <  6:1430403/137‑1 (MQ=11)
                                                                               ggCTTATCGGTCAGTTTCACCTGATTACCGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTTGCAGGCGCAGAAAGATGAAGGACTGGCCACGACGCTATACCCAAACGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAcc                                                                 <  7:1845602/151‑1 (MQ=11)
                                                                                gcgtATGGGACAGTTGCCCGTGATTTACGAAACACCCCGCTGCGGCGGGTTTTTGCTTTTGCCGGGGCCCAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAAccc                                                                <  1:1501531/148‑1 (MQ=11)
                                                                                   tATCGGCAAGTTTCCCCTGATTTAGGTAACACCCCGCTTCGGCGGGTTTTTGCTTTTGTTGGGGCCGAACCATGATGGCCTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGc                                                              <  5:1188361/150‑1 (MQ=11)
                                                                                    atgggaatgTTTCCCCTCATTACCTTAAAAACCCGCTGCGGCGGGTTTTTGCTTTTGCCGGGGCCTAAAGTTGATTGCCTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                                            <  6:666360/143‑1 (MQ=11)
                                                                                     tcgttaattttccgatGATTTAGTTAAAACCCGTCTGCGGCGGGTTTTTGTTTTTGTAGGGGATTAATTATGAATTACTGGCCACGTCTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCtt                                                            <  2:495792/135‑1 (MQ=2)
                                                                                        ggaaGTTGCACGTGATTTCCTTAAACACGGGCTGCGGCGGGTTTTTGGTTTTGTAGGGGCAAAAAGATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTcggc                                                        <  1:338052/148‑1 (MQ=11)
                                                                                              ttCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGCTGTACGTAAAAACCCGGTGCGGCGGGttt                                                  >  7:374853/1‑151 (MQ=255)
                                                                                              tgccccTGATTTACGAAAAAACCCGCTTCGGCGGGTTTTTGGTTTTGCAGGGCCAGAAAGATGAAGCACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTACGTAAAAACCCGCTTCGGCGGGttt                                                  <  3:416740/146‑1 (MQ=11)
                                                                                                cACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGAGGAATGACGGTCCACCACCCTAGACCCAAAAGAAAGCGGCGTATCGGTCAGTATGACCTGGTGTACGTAAAAAACCGATGCGGCGGGtttcg                                                >  5:1128102/1‑149 (MQ=11)
                                                                                                  ccgcaTTTACGTAACCCCCGGCGGCGGCGGGTTTTTGCGTTTGGCGGCGCACACAGATGCATGACTCTCCACGACGCTATACCCAACAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGc                                              <  5:869559/147‑1 (MQ=11)
                                                                                                     gATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCCCGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTAACATGGTGTACGTAATAAACCGCTTCGGCGGGTTTGTGCttg                                           >  7:1138773/1‑150 (MQ=11)
                                                                                                      aTTTCCGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTTGGCGGCGCCGAAACATGAAGGCCTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                          <  5:80562/151‑1 (MQ=11)
                                                                                                      aTTTACGTAAACACCCGCTTCCGCGGGTTTTTGCTTTTGGAGCGGCCGAAAGATGAATGACTCGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCtttt                                          <  1:579939/151‑1 (MQ=11)
                                                                                                       tttACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTGACCTGGTGTACGTAATAACCCGGTCCGGCgggagttgggttg                                           >  1:360543/1‑139 (MQ=11)
                                                                                                         tccgaaaaacCCCGCTGCGGCGGGTTTTTGCTTTTGGCGGGGCCCAAACATGAATGACTGGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTgg                                        <  3:460829/145‑1 (MQ=11)
                                                                                                          ccGTAACACCCCGCTGCGCCGGGTTTTTGGTTTTGCCGGCGCCGAACAATGAAGGCCTCGCCACGACGCTCTACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAg                                      <  1:1355550/150‑1 (MQ=11)
                                                                                                          aCGTAAAACCCCGCTGCGGCGGGTTTTTGCTTTGGTCGGGGCAGAACGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGa                                       <  5:1517228/150‑1 (MQ=11)
                                                                                                             tacacacgcGCTTCGGCGGGTTTTTGCTTTTGTTGCTGCTGAAATATTATTGACTGTCCACTACTCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAggg                                    <  2:256448/142‑1 (MQ=11)
                                                                                                             aaaaacCCCGCTGCGGCGGGTTTTTGCTTTTGCCGGGGCACAAAGATGAAGGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgggg                                   <  5:1813225/150‑1 (MQ=11)
                                                                                                               acaaCCCGCTTGCTCGGTTTTTTGCTTTTGGAGTTTTATTTTTATGATTGTCCGTCCACGTCGCTATACACATACTTAATCTGGTTATCGGTAACTTTCAGCTGTTTTTCTTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGc                                  <  6:580968/148‑1 (MQ=2)
                                                                                                               aaaaCCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGAAGGATCACTGGGCACGACGCGATACCCAAAACAAACGGGCTTATCGGTGAGTTTCCCCGGGTTTACGTAAATCCCGGCTTCGGGGGGTGTTGGCGTGGGGAgggggc                                 >  5:580922/1‑149 (MQ=2)
                                                                                                                 aaCCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGGAGAAAGATGAATGACTGTCCAAAACGCTATACCCAAAAGAAAGAGGATTATCGGTCAGTTTCACATATTTTAAGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAgggggcgg                               >  6:964705/1‑147 (MQ=11)
                                                                                                                   acccgtcttatgcggGTTTTTGTTTTTGTTGGTGCTTAATTTTGATTGACTTGCCAAGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGa                               <  8:1068647/140‑1 (MQ=11)
                                                                                                                       gCTGCGGCGGGTTTTTGCTTTTGTTGGGGCTTAAACTTGATTGACTTGCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAg                            <  8:1403248/149‑1 (MQ=11)
                                                                                                                           tCGGCGGGTTTTTGTTTTTGGAGGGGCATAAATATGAAGGACTGGCCACGACGCTATACCCAAAGGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGatg                         <  7:1585724/149‑1 (MQ=11)
                                                                                                                                 ggTTTTTGCTTTTGCAGGGGCAGAAACATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTg                 <  5:236931/151‑1 (MQ=25)
                                                                                                                                            ttGGAGGGGCAGAAAGATGAATGACTGTCCACGACCCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTTACCTGATGTACGTAATAAACCGTTCCGGCGGGTTTAAGATTGTGGAGGGCCCGTAAGACAAGAGCCAGTCCACGacagc      >  8:1792311/1‑149 (MQ=17)
                                                                                                                                                aGGGGCAGAAAGATGAATGAATGTCCAAGACGCTATAACAAAAAGAAAGCGGATTATCGATCAGTTTAACATGGTTTACGTAAAAACAATATTATAAGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACGGTcggcagcgcgtggg  >  6:466763/1‑138 (MQ=2)
                                                                                                                                                       |                                                                                                                                               
CCCGTACTGAACCCGGGTAGGCTGCTTGAGCCAGTGAGCGATTGCTGGCCTAGATGAATGACTGTCCACGACAGAACCCGGCTTATCGGTCAGTTTCACCTGATTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACGCTATACCCAAAAGAAAGCGGCTTATCGGTCAGTTTCACCTGGTTTACGTAAAAACCCGCTTCGGCGGGTTTTTGCTTTTGGAGGGGCAGAAAGATGAATGACTGTCCACGACACTATAC  >  DH1/613448‑613740

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