Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3280CA44.7% 101.1 / 51.0 76intergenic (+337/‑76)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (40/2);  new base A (0/34);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.39e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  NC_000913/345180‑345476
                                                                                                                                                      |                                                                                                                                                    
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     >  1:1076556/1‑150 (MQ=255)
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     >  1:879957/1‑150 (MQ=255)
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                                                     >  1:571153/1‑150 (MQ=255)
   aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaaa                                                                                                                                                  >  2:1010698/1‑150 (MQ=255)
         ctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc                                                                                                                                            >  2:16625/1‑150 (MQ=255)
         ctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc                                                                                                                                            >  1:758047/1‑150 (MQ=255)
            aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTCCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                         >  2:93923/1‑150 (MQ=255)
                 gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc                                                                                                                                    >  1:467098/1‑150 (MQ=255)
                 gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtc                                                                                                                                    >  2:248626/1‑150 (MQ=255)
                    tttttttGAATGTGTTGGTTGAAATTAATATAAAAAAAATCATATTCCCATCATTTTTGTGTGGGGTTTTTTCTTATTCTTACGCAGCCCACTCAACCTTAACAAAAAATATGTTAAAAGCCTCATATAAAAAACCGCCATTAGCtcatc                                                                                                                                 <  1:653313/149‑1 (MQ=255)
                     ttttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                 >  2:259552/1‑149 (MQ=255)
                            attgttTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCATTTTTGTGGGGGGTTTTTTTTGTATTTTACGCATCGCACCCAAGCCGGAAAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAaga                                                                                                                          <  2:1253552/147‑1 (MQ=255)
                              tgtgTTGGTTGACATTAATATAAAAAAAACAATAATCCCATAATTTTTGTGGGGGTTTTTTTCTTTATCTTCCCCATCGCACACAAGCGTGACAAAAAATATGCTGTAAGGCTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                       <  1:338928/150‑1 (MQ=255)
                                  ttggttgATATTAATAAAAAAAAAATAATAATCCCATCTTGTTTGTGGGGGTTTTTTTCTTTATCTTCCCCCCCCCACCCAACCCTGAAAAAAAATATGCTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                       <  1:609875/146‑1 (MQ=255)
                                  ttggttcACTTTAATAAAAAAAAAAAAATAATCCCATCATTTTTGTGGGGGTTTTTTTTTTTTTCTTCCCCACCCCCCCAAAGCCGGAAAAAAAATATGCTGAAGGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                        <  1:676596/145‑1 (MQ=255)
                                    ggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                 >  2:437286/1‑150 (MQ=255)
                                      ttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              >  2:1071194/1‑151 (MQ=255)
                                       tgACATTCATATAAAAAAAATCATAATTCCATCATGTTTGGGGGGGTTTTTTTCTGTATCTTACGCACCGCACTCAACCCTAAAAAAAAATATGCTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              <  2:545492/150‑1 (MQ=255)
                                        gACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                             >  2:21710/1‑150 (MQ=255)
                                         actttCATATAAAAAAAATCATAATTCCATCATTTTTGTGTGGGGTCTTTTCTTTTCCTTACGCACCGCACCCAACCCTAAAAAAAAATAGGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa                                                                                                            <  2:873356/147‑1 (MQ=255)
                                          cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                          >  2:658142/1‑151 (MQ=255)
                                            ttCATATAAAAAAAATCATATTCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTACCCACCCCACCAAAGCCTGAAAAAAAATATGTGGAAAGCCTCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              <  2:823711/145‑1 (MQ=255)
                                               ataaaaaaaaaaTCATAATCCCATCTTTTTTGTGGGGGTTTTTTTCTTTTTCTACCGCACCGCCCTCACGCCTAAAAAAAAATATGCTGAAAGCCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                                                        <  1:808714/143‑1 (MQ=255)
                                                   aaaaaaaaTCATAATTCCATCATGTTTGTGGGGGGTTTTTTCTGTATCTTCCGCACCCCACTCAACCCTAACAAAAAATATGCTGTAAGGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                       <  2:660796/144‑1 (MQ=255)
                                                   aaaaaaaaTCATAATCCCATCTTTTTTGGGTGGGTTTTTTTTTGTATCTTCCGCATCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgtt                                                                                                    <  1:616872/147‑1 (MQ=255)
                                                   aaaaaaaaCAAAAATCCCATCTTTTTGGTGTGGGTTTTTTTTTATTCTTACCGCCCCCCACCAAACCCTGAAAAAAAATATGGTGAAGGCAAAATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa                                                                                                 <  1:412360/150‑1 (MQ=255)
                                                     aaaaaaTCATAATCCCATCTTTTTTGTGGGGGTTTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                               <  2:529651/151‑1 (MQ=255)
                                                     aaaaaaTAATAATCCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                               <  2:734548/151‑1 (MQ=255)
                                                        aaaTTATAATTCCATCATTTTGGTGGGGGATTTTTTTTATATCTTACGCACCCCACTCAAGCCTGAAAGAAAATATGCTGAAATTCTCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                            <  2:880429/151‑1 (MQ=255)
                                                        aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                              >  1:647812/1‑149 (MQ=255)
                                                        aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                             >  2:130624/1‑150 (MQ=255)
                                                             aaaactccttattttttGTGTGGGGTCTTTTCTTTACTTTACGCTTCCCACACATGCCGGAAAAAAAATATCCTGTAAGGCTAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                            <  2:370910/132‑1 (MQ=255)
                                                                   ccATCATTTTTGTGGGGGTTTTTTTTTGTATCTTACCCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                     <  1:586497/147‑1 (MQ=255)
                                                                   ccATCATGTTTGTGTGGGGTTTTTTCTGTATCTTCCGCACCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa                                                                                  <  2:207320/150‑1 (MQ=255)
                                                                       catGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGGTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                             >  2:1268897/1‑151 (MQ=255)
                                                                        atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                            >  2:1286685/1‑151 (MQ=255)
                                                                           tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                         >  1:1261741/1‑151 (MQ=255)
                                                                                 tgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                   >  2:800535/1‑151 (MQ=255)
                                                                                  ggggTTTTTTCTGTATCTTCCGCTCCGCACTCAAGCCTGAAAAAAAATAGGCTGTAGGGCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                  <  2:308323/151‑1 (MQ=255)
                                                                                      ttttttCTTTATCTTCCGCACCGCACAAAACCCTGACAAAAAATATGCTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcct                                                                <  1:899993/147‑2 (MQ=255)
                                                                                      tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                >  2:67111/1‑149 (MQ=255)
                                                                                      tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                              >  2:815625/1‑151 (MQ=255)
                                                                                        ttttCTGTATTTTACGCATCCCACTCAAGCCTGAAAGAAAATATGCTGAAGGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                             <  2:249887/150‑1 (MQ=255)
                                                                                         tttttttttCTTACGCACCCCACTCAACCCTGAAAAAAAATATGCTTTAAGGCTAAAAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                            <  2:678854/142‑1 (MQ=255)
                                                                                          ttttttATCTTCCGCACCCCACTCAACCCAAAAAAAAATTTTGCTGTAAGCCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                           <  2:794529/145‑1 (MQ=255)
                                                                                          ttCTGTATCTTCCGCATCGCACTAAAGCCTGAAAAAAAATTTGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                           <  2:155041/150‑1 (MQ=255)
                                                                                              gTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                                       >  2:1249956/1‑150 (MQ=255)
                                                                                               tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                      >  1:127535/1‑150 (MQ=255)
                                                                                                aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                    >  1:119857/1‑151 (MQ=255)
                                                                                                 tCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCACAAct                                                   >  2:763505/1‑151 (MQ=255)
                                                                                                 tCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                                   >  2:400716/1‑151 (MQ=255)
                                                                                                     aCGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                               >  2:596451/1‑151 (MQ=255)
                                                                                                           tCGCACCCAAGCCGAAAAAAAAAGATGCTGAAAGCCTAATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                  <  1:689828/143‑1 (MQ=255)
                                                                                                            cccaccaaagctcgACAAAAAATGTGCTGGAAGCCCCAAATAAAACGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                            <  1:1070163/135‑1 (MQ=255)
                                                                                                             gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                         >  2:921439/1‑150 (MQ=255)
                                                                                                                  ccaGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                   >  1:134515/3‑151 (MQ=255)
                                                                                                                   agcccTGAAAAAAAATTGGCTGTAGGCCCAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                <  1:800476/134‑1 (MQ=255)
                                                                                                                    aGCCTTACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc                                 >  2:1142805/1‑151 (MQ=255)
                                                                                                                      cctaaaaaaaaaTTTGCTGTAAGGCTCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                               <  2:390981/143‑1 (MQ=255)
                                                                                                                       cTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                               >  1:1143444/1‑150 (MQ=255)
                                                                                                                         gACAAAAAATGTGCTATAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                              <  1:830186/149‑1 (MQ=255)
                                                                                                                          aaaaaaaaTATCCTGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                              <  2:1118447/144‑1 (MQ=255)
                                                                                                                          aCAAAAAATAGGCTGAAGGCCAAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                             <  1:737879/149‑1 (MQ=255)
                                                                                                                           cAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                          >  2:827430/1‑151 (MQ=255)
                                                                                                                            aGAAAAGATGCTGTAGGGCTGATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                          <  2:594594/150‑1 (MQ=255)
                                                                                                                              aaatTATCCGGTAAGGCAAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                       <  2:831750/151‑1 (MQ=255)
                                                                                                                               aaaTATGTTGTAAGGCCCATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcatc                      <  1:1210064/151‑1 (MQ=255)
                                                                                                                                    tGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                   <  2:1076625/149‑1 (MQ=255)
                                                                                                                                     gCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                 >  2:1056621/1‑150 (MQ=255)
                                                                                                                                       tGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                >  1:603303/1‑149 (MQ=255)
                                                                                                                                          aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg           >  2:1078230/1‑151 (MQ=255)
                                                                                                                                          aaGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt            <  2:827288/150‑1 (MQ=255)
                                                                                                                                           aGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg          >  1:1196300/1‑151 (MQ=255)
                                                                                                                                                   tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTg   >  1:543164/1‑151 (MQ=255)
                                                                                                                                                   tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTg   >  2:1016789/1‑151 (MQ=255)
                                                                                                                                                    atCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGc  >  2:460284/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  NC_000913/345180‑345476

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