Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,255,6620CT36.7% 175.3 / 90.7 128intergenic (+144/+75)dgcZ/yneMdiguanylate cyclase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/73);  new base T (47/0);  total (55/73)
Fisher's exact test for biased strand distribution p-value = 4.73e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

AGTTTAAGCTCCGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATAAA  >  DH1/2255512‑2255811
                                                                                                                                                      |                                                                                                                                                                  
 gTTTAAGCTCCGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCg                                                                                                                                                                   <  1:1726068/151‑1 (MQ=255)
     ctgCTCCGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAACCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGcc                                                                                                                                                                 <  2:36652/147‑1 (MQ=255)
     aaGCTCCGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCg                                                                                                                                                               <  2:1293131/150‑1 (MQ=255)
           cGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCCTaa                                                                                                                                                       >  2:1731867/1‑150 (MQ=255)
            gTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATaaaa                                                                                                                                                     <  2:1347406/151‑1 (MQ=255)
              ttAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATaa                                                                                                                                                       <  1:1404301/147‑1 (MQ=255)
               tAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAt                                                                                                                                                   >  1:1253621/1‑150 (MQ=255)
                 aCATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAg                                                                                                                                               <  2:580951/151‑1 (MQ=255)
                   attcattGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGc                                                                                                                                             <  2:351174/151‑1 (MQ=255)
                     tcattGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGcc                                                                                                                                            <  2:1253623/150‑1 (MQ=255)
                      cattGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGccc                                                                                                                                           <  2:1132135/150‑1 (MQ=255)
                         tGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCtt                                                                                                                                       <  2:2592494/150‑1 (MQ=255)
                         tGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCtt                                                                                                                                       <  2:1762910/150‑1 (MQ=255)
                         tGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCtt                                                                                                                                       <  2:550041/150‑1 (MQ=255)
                          gagaAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGCCATAAAAAAAAAGcccccttt                                                                                                                                      >  2:1446060/1‑147 (MQ=255)
                          gagaAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTg                                                                                                                                      <  2:2508467/150‑1 (MQ=255)
                          gagaAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTg                                                                                                                                      <  2:1404774/150‑1 (MQ=255)
                             aaaaCTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGCCATaaaaaaaa                                                                                                                                                 >  2:1522210/1‑135 (MQ=255)
                             aaaaCTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGTCGGCCTAAAAAAAAAGGCCCTTTTac                                                                                                                                   >  1:1546604/1‑150 (MQ=255)
                             aaaaCTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCa                                                                                                                                                <  1:1522205/138‑1 (MQ=255)
                             aaaaCTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTCTTTTTTTTCCCGCCTAAAAAACACGCCCCCTGtttc                                                                                                                                  >  1:1904966/1‑148 (MQ=255)
                                   gATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCTGGCCTAAAAAACAGGCCCCCTTGACACaaatc                                                                                                                           >  1:1088645/1‑148 (MQ=255)
                                    aTGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATTAAAAAAAGGCCCCCTTTTCACAACTGt                                                                                                                         >  2:972304/1‑151 (MQ=255)
                                    aTGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGt                                                                                                                         <  2:125585/151‑1 (MQ=255)
                                         aCTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAAAAAGCCCCCTTTTCACAAATGTaaaa                                                                                                                    >  1:469513/1‑148 (MQ=255)
                                         aCTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTAAAAAAAAGGCCCCCTTTACAAAAATGTaaaa                                                                                                                    >  2:858911/1‑148 (MQ=255)
                                           tgtgTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAAAAGGGCCCCTGTtca                                                                                                                                  >  2:2109142/1‑137 (MQ=255)
                                           tgtgTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTaca                                                                                                                                  <  1:2109134/137‑1 (MQ=255)
                                              gtCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTAAAAAAAAAGCCCCCTTTACCCAAATTTAAAAAagggg                                                                                                            >  2:2559669/1‑150 (MQ=255)
                                                 aaCGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCgg                                                                                                        <  2:1546609/151‑1 (MQ=255)
                                                     aTTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCCTCATAAAAAAAAGGCCCCTTTTTCACAACTCTAACAAAGGGGCGGtttg                                                                                                    >  1:251728/1‑149 (MQ=255)
                                                        cGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCCTAAAAAACAAGGCCCCTTTACACAACTGTAACAAAGGGGCGGGTTgggg                                                                                                  >  1:895111/1‑150 (MQ=255)
                                                          ttATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTTCGTCCTTAAAAAAAAGGCCCCTTTTTACAAAACTTAAAAAGGGGCGGGttggtggg                                                                                                >  1:2149000/1‑145 (MQ=255)
                                                          ttATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACCGGCGGTTTTTTTTTTCCGTTCTAAAAAAAAAGGCCCTTTTACTCACCAGTAAAAAAGGGGCGGGTTgggg                                                                                                  >  1:2631197/1‑148 (MQ=255)
                                                            aTCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGCCCTCATAAAAAAAAAGCCCCCTGTTCACCAATGGTAAAAAGGGCCCGGTTGGTGgggg                                                                                             >  2:1394580/1‑151 (MQ=255)
                                                            aTCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAgg                                                                                              <  1:972303/150‑1 (MQ=255)
                                                                gTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTTCGGCATAAAAAAAAAGACCCCTGTTCACCAACTCAAAAAAGGGGCGGGTAGGGGAGGGGttt                                                                                         >  1:2365472/1‑150 (MQ=255)
                                                                  gCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGGTTTTTTTTTCCGTCATAAAAAATAAGACCCTTTTTCACAACTGTAAAAAAGGGGCGGGTTGGTGAGGGAGTAtt                                                                                       >  1:2335119/1‑150 (MQ=255)
                                                                    aGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCCTCCTAAAAAAAAAGCCCCCTTTTCACAAAAGTTAAAAAGGGGCGGGTTGGGGAGGGGAtttttt                                                                                      >  1:1986721/1‑147 (MQ=255)
                                                                     ggAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTcc                                                                                    <  2:626247/151‑1 (MQ=255)
                                                                      gAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCTTCATAAAAAACAGGCCCCTTGTACGCCACTTTAACCACGGGGCGGGTATGTGAggggatgtgttc                                                                                    >  1:2405260/1‑141 (MQ=255)
                                                                       aaaaTGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGt                                                                                  <  2:1904971/151‑1 (MQ=255)
                                                                        aaaTGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGtt                                                                                 <  2:2149008/151‑1 (MQ=255)
                                                                          aTGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCa                                                                               <  2:2631213/151‑1 (MQ=255)
                                                                            gCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCAt                                                                              <  2:2405271/150‑1 (MQ=255)
                                                                             ccTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTTAAAAACAAGGCCCCTTTACACCAACGTAACAAAGGGGCGGGTAGGGGGGGGGGTATTTCccttccttt                                                                            >  2:2289632/1‑150 (MQ=255)
                                                                             ccTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTa                                                                            <  1:858910/151‑1 (MQ=255)
                                                                               tGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCg                                                                                                         <  1:575536/120‑1 (MQ=255)
                                                                                gTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCCTCATAAAAATAAAGCCCCCTTTTTACAAAAGTAAAAAAGGGGCGGGTTGGTGAGGGGGTATTTCCGtttctttttt                                                                         >  2:477160/1‑143 (MQ=255)
                                                                                gTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGt                                                                          <  2:837065/150‑1 (MQ=255)
                                                                                     cGTAAAAGCAAAACGCAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGCCCTAAAAAAAAAGCCCCCTTTAAACAAATGTAAAAAAGGGCCGGGTTAGGGGGGGGGTTTTTTCCtttattttttattcc                                                                     >  1:634434/1‑139 (MQ=255)
                                                                                     cGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTAAAATATAAGCCCCTTGTACACCACTGGAACAACAGGGCGGGTTTGGGAGGGGTTATCTTCGCTCTTTTATCATccc                                                                   >  1:2602652/1‑151 (MQ=255)
                                                                                     cGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCCTCATTAAAAACAAGCCCCTTTTTCACAAATGTAAAAAAGGGGCGGGTAGGGGAGGGGTTATTTCCCTTCCTTATTTATccc                                                                   >  2:1087090/1‑151 (MQ=255)
                                                                                     cGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATccc                                                                   <  1:2289623/151‑1 (MQ=255)
                                                                                      gTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAgg                                                                                                   <  2:163614/119‑1 (MQ=255)
                                                                                             cAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTa                                                                            <  1:2231098/135‑1 (MQ=255)
                                                                                                 acacAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGGCGGCCTAAAAAACAAGCCCCCTCTACACAACCTTAAAAAAGGGACGGGTAGGTGGGGGGATATTTTCCTTCCTTTGTTCTCCCAttttttggtgg                                                       >  2:915977/1‑143 (MQ=255)
                                                                                                 acacAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGa                                                       <  1:573903/151‑1 (MQ=255)
                                                                                                  cacaAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGaa                                                      <  2:469513/151‑1 (MQ=255)
                                                                                                   acaAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTAAAAAAAAAGACCCCTCTTCCCCAACCTAACAAGGGGGCGGGTTGGGGAGGGGTTTTTTCCCTTTATTTGTCATCCCCCTTGTgggttta                                                      >  2:16565/1‑145 (MQ=255)
                                                                                                    caAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGCCGTCATAAAAAACAAGCCCCCTTTACACCAATGTAAAAAAGGGGCGGGTTGGTGAGGGGTTATTTCCGTTCATTTGTCATTCCCTTTTTTGCTGaaa                                                     >  2:2662672/1‑150 (MQ=255)
                                                                                                     aaaTCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGt                                                                          <  1:2075067/129‑1 (MQ=255)
                                                                                                      aaTCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATAc                                                  <  2:634435/151‑1 (MQ=255)
                                                                                                       aTCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGGTTTTTTTTTCCCTCCTCAAAAAAAAGACCCCTTTACACAAATGTAACAAAGGGGCGGGTTGGGGGGGGGTTAACTCCTTTCCTTTGTCCTCCCATTTTTTGGTGaaaaaa                                                  >  1:466762/1‑147 (MQ=255)
                                                                                                         cTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCCTAAAAAAAAAGCCCCCTTTTCACCCCAGTAAAAAAAGGGCCGGTTGGGGGGGGGTTATCTCCCTTCATTAGTCATCCCCATTTTTGGTGTAAAAcgcg                                              >  2:2641830/1‑151 (MQ=255)
                                                                                                          tATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATAc                                                  <  1:1435423/147‑1 (MQ=255)
                                                                                                           aTCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCCTAAAAAACAGGGCCCCTTTATACAAAAGTAAAAAAGGGGCGGGTTGGGGGGGGGGTTTCTCCCTTTCTTTTTTATCCCCTTTTTGGGTGGAATAcgaga                                              >  1:1828048/1‑147 (MQ=255)
                                                                                                           aTCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGTCCTCCTAAAAAAAAAGCCCCCTGTTCACAACCGTAACAAAGGGGCCGGTGGGTGGGGGGGTATTTCCCCTTATTAATTATTCCATTTTTGGGTGtaaaaaacgcc                                             >  2:2267282/1‑143 (MQ=255)
                                                                                                           aTCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATAcgcg                                              <  2:2510954/149‑1 (MQ=255)
                                                                                                           aTCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCgg                                             <  2:1986729/150‑1 (MQ=255)
                                                                                                           aTCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGc                                            <  2:1781423/151‑1 (MQ=255)
                                                                                                            tCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTGCCGGCATAAAAAAAAGGCCCCCTTTACACAACTGTAACAAAGGGGCGGGTAGGTGGGGGGTTATTTCCGCTCATTAGTCATCCCCTTTTTTGCTGAAAAAAGCGGcc                                           >  2:2342522/1‑151 (MQ=255)
                                                                                                              cATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAAAAAGCCCCCTTTTCACAAATGTAAAAAAGGGGCGGGTTGGTGAGGGGTTATTTCCGTTCATTTTTTATCCCCTTTTTGGGTGGAAAAAGCGGcccg                                         >  2:1445889/1‑149 (MQ=255)
                                                                                                               aTGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGGCATAAAAAACAAGCCCCCTTTTCACCAACGTAACAAAGGGGCGGGTTGGGGGGGGGTTATTTCCCTTCATTTATTATCCCCTTTTTTGGTGGAATaagagggc                                            >  1:1665071/1‑141 (MQ=255)
                                                                                                                 gCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaaa                                      <  2:1054711/151‑1 (MQ=255)
                                                                                                                 gCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaa                                       <  1:2559654/150‑1 (MQ=255)
                                                                                                                  cAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAACAAGCCCCCTGTTCACAAATGTAAAAAAGGGGCGGGTAGGGGAGGGGTTTTTTCCGTTCCTTAGTCATTCCATTTTTGGGTGAAATAAGCGGcccgaaa                                      >  2:204695/1‑150 (MQ=255)
                                                                                                                  cAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaaa                                      <  1:204696/150‑1 (MQ=255)
                                                                                                                   aaGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATAcgc                                                <  1:1101830/140‑1 (MQ=255)
                                                                                                                       aTTCACCGCCGGTTTACTGGCGGTTTTTTTTTGCCGGCCTAAAAAAAAGGCCCCCTTTTCACAAATGTAAAAAAGGGGCGGGTAGGGGGGGGGGTTTCTCCGTTTTTTTGTCATCCCCTTTTTGGGTGGAAAAAGCGGcccgatcggg                                   >  2:865324/1‑143 (MQ=255)
                                                                                                                         tCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaaaaccagaa                               <  2:2517357/150‑1 (MQ=255)
                                                                                                                          cACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaaaaccagaaaa                             <  2:2335129/151‑1 (MQ=255)
                                                                                                                            ccgccgGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTtgt                                                             <  1:292604/118‑1 (MQ=255)
                                                                                                                            ccgccgGTTTACTGGCGGTTTTTTTTCGCCCTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccagaaaacca                                  <  2:1012550/144‑1 (MQ=255)
                                                                                                                                 ggTTTACTGGCGGTTTTTTTTTTCCGTCATAAAAAACAAGCCCCTTTTACACAACTGTTACAAAGGGGCGGGTTGGTGAGGGGATATTTTCCTTCATTTGTTATCCCCTTTTTTGGTGGAAAAAGCGGGCCGAAAAACCAAaaaaaaaaat                       >  1:864175/1‑151 (MQ=255)
                                                                                                                                      aCTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTcc                 <  1:1445884/151‑1 (MQ=255)
                                                                                                                                        tGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGccaga                                        <  2:396866/126‑1 (MQ=255)
                                                                                                                                        tGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGAATATCTCCGTTCATTAGTCATCCCATTTGTGGATGAAATACGAGGtcaga                                        >  1:396867/1‑126 (MQ=255)
                                                                                                                                         ggCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAatt                      >  2:2272000/1‑143 (MQ=255)
                                                                                                                                         ggCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAatt                      <  1:2271991/143‑1 (MQ=255)
                                                                                                                                         ggCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCa               <  1:2267273/150‑1 (MQ=255)
                                                                                                                                         ggCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAg              <  1:1907131/151‑1 (MQ=255)
                                                                                                                                               tttttttCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGc        >  2:2495865/1‑151 (MQ=255)
                                                                                                                                                 tttttCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGcc       >  2:149937/1‑150 (MQ=255)
                                                                                                                                                 tttttCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGcc       <  2:1840025/150‑1 (MQ=255)
                                                                                                                                                     tCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATaa   >  2:2419551/1‑150 (MQ=255)
                                                                                                                                                      cGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATaaa  <  1:1459944/150‑1 (MQ=255)
                                                                                                                                                      cGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATaaa  <  1:915976/150‑1 (MQ=255)
                                                                                                                                                      cGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATaaa  >  1:1642276/1‑150 (MQ=255)
                                                                                                                                                      cGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATaaa  <  2:1854622/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                                  
AGTTTAAGCTCCGTTTAACATTCATTGAGAAAACTGATGCTACTGTGTCAACGATTCGTTATCAGTGCAGGAAAATGCCTGTTAGCGTAAAAGCAAAACACAAATCTATCCATGCAAGCATTCACCGCCGGTTTACTGGCGGTTTTTTTTCGCCGTCATAAAAATCAGGCCCCTTGTACACAACTGTAACAAGGGGCCGGTTAGGTGAGGGATTATCTCCGTTCATTAGTCATCCCATTTGTGGCTGAAATACGCGGCCAGAAAACCAGAAAATAAAATTATTCCCAGTACGGCCATAAA  >  DH1/2255512‑2255811

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