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

AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  NC_000913/345180‑345476
                                                                                                                                                      |                                                                                                                                                    
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCa                                                                                                                                                     >  1:1076556/1‑150 (MQ=255)
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCa                                                                                                                                                     >  1:879957/1‑150 (MQ=255)
aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCa                                                                                                                                                     >  1:571153/1‑150 (MQ=255)
   aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCaaaa                                                                                                                                                  >  2:1010698/1‑150 (MQ=255)
         ctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGcc                                                                                                                                            >  2:16625/1‑150 (MQ=255)
         ctcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGcc                                                                                                                                            >  1:758047/1‑150 (MQ=255)
            aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTCCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCAtt                                                                                                                                         >  2:93923/1‑150 (MQ=255)
                 gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCtc                                                                                                                                    >  1:467098/1‑150 (MQ=255)
                 gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCtc                                                                                                                                    >  2:248626/1‑150 (MQ=255)
                    tttttttGAATGTGTTGGTTGAAATTAATATAAAAAAAATCATATTCCCATCATTTTTGTGTGGGGTTTTTTCTTATTCTTA‑CGCAGCCCACTCAACCTTAACAAAAAATATGTTAAAAGCC‑TCATATAAAAAACCGCCATTAGCtcatc                                                                                                                                 <  1:653313/149‑1 (MQ=255)
                     ttttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                 >  2:259552/1‑149 (MQ=255)
                            attgttTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCATTTTTGTGGGGGGTTTTTTTTGTATTTTA‑CGCATCGCACCCAAGCCGGAAAAAAAATATGCTGTAAGGC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAaga                                                                                                                          <  2:1253552/147‑1 (MQ=255)
                              tgtgTTGGTTGACATTAATATAAAAAAAACAATAATCCCATAATTTTTGTGGGGGTTTTTTTCTTTATCTTC‑CCCATCGCACACAAGCGTGACAAAAAATATGCTGTAAGGC‑TAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                       <  1:338928/150‑1 (MQ=255)
                                  ttggttgATATTAATAAAAAAAAAATAATAATCCCATCTTGTTTGTGGGGGTTTTTTTCTTTATCTTC‑CCCCCCCCACCCAACCCTGAAAAAAAATATGCTGAAAGCC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                       <  1:609875/146‑1 (MQ=255)
                                  ttggttcACTTTAATAAAAAAAAAAAAATAATCCCATCATTTTTGTGGGGGTTTTTTTTTTTTTCTTC‑CCCACCCCCCCAAAGCCGGAAAAAAAATATGCTGAAGGCC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                        <  1:676596/145‑1 (MQ=255)
                                    ggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                 >  2:437286/1‑150 (MQ=255)
                                      ttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              >  2:1071194/1‑151 (MQ=255)
                                       tgACATTCATATAAAAAAAATCATAATTCCATCATGTTTGGGGGGGTTTTTTTCTGTATCTTA‑CGCACCGCACTCAACCCTAAAAAAAAATATGCTGTAAGGC‑TAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              <  2:545492/150‑1 (MQ=255)
                                        gACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                             >  2:21710/1‑150 (MQ=255)
                                         actttCATATAAAAAAAATCATAATTCCATCATTTTTGTGTGGGGTCTTTTCTTTTCCTTA‑CGCACCGCACCCAACCCTAAAAAAAAATAGGCTGAAAGGC‑TCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa                                                                                                            <  2:873356/147‑1 (MQ=255)
                                          cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                          >  2:658142/1‑151 (MQ=255)
                                            ttCATATAAAAAAAATCATATTCCCATCATTTTTGTGTGGGGTTTTTTCTGTATCTTA‑CCCACCCCACCAAAGCCTGAAAAAAAATATGTGGAAAGCC‑TCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                              <  2:823711/145‑1 (MQ=255)
                                               ataaaaaaaaaaTCATAATCCCATCTTTTTTGTGGGGGTTTTTTTCTTTTTCTAC‑CGCACCGCCCTCACGCCTAAAAAAAAATATGCTGAAAGCC‑TCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                                                        <  1:808714/143‑1 (MQ=255)
                                                   aaaaaaaaTCATAATTCCATCATGTTTGTGGGGGGTTTTTTCTGTATCTTC‑CGCACCCCACTCAACCCTAACAAAAAATATGCTGTAAGGC‑TCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                                                       <  2:660796/144‑1 (MQ=255)
                                                   aaaaaaaaTCATAATCCCATCTTTTTTGGGTGGGTTTTTTTTTGTATCTTC‑CGCATCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGCC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgtt                                                                                                    <  1:616872/147‑1 (MQ=255)
                                                   aaaaaaaaCAAAAATCCCATCTTTTTGGTGTGGGTTTTTTTTTAT‑TCTTACCGCCCCCCACCAAACCCTGAAAAAAAATATGGTG‑AAGGCAAAATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa                                                                                                 <  1:412360/150‑1 (MQ=255)
                                                     aaaaaaTCATAATCCCATCTTTTTTGTGGGGGTTTTTTTCTGTATCTTA‑CGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                               <  2:529651/151‑1 (MQ=255)
                                                     aaaaaaTAATAATCCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGC‑TAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg                                                                                               <  2:734548/151‑1 (MQ=255)
                                                        aaaTTATAATTCCATCATTTTGGTGGGGGATTTTTTTTATATCTTA‑CGCACCCCACTCAAGCCTGAAAGAAAATATGCTGAAATTC‑TCATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                            <  2:880429/151‑1 (MQ=255)
                                                        aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                              >  1:647812/1‑149 (MQ=255)
                                                        aaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                             >  2:130624/1‑150 (MQ=255)
                                                             aaaactccttattttttGTGTGGGGTCTTTTCTTTACTTTA‑CGCTTCCCACACATGCCGGAAAAAAAATATCCTGTAAGGC‑TAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                            <  2:370910/132‑1 (MQ=255)
                                                                   ccATCATTTTTGTGGGGGTTTTTTTTTGTATCTTA‑CCCACCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                     <  1:586497/147‑1 (MQ=255)
                                                                   ccATCATGTTTGTGTGGGGTTTTTTCTGTATCTTC‑CGCACCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGGC‑TCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa                                                                                  <  2:207320/150‑1 (MQ=255)
                                                                       catGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGGTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                             >  2:1268897/1‑151 (MQ=255)
                                                                        atGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                            >  2:1286685/1‑151 (MQ=255)
                                                                           tttGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                         >  1:1261741/1‑151 (MQ=255)
                                                                                 tgGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                   >  2:800535/1‑151 (MQ=255)
                                                                                  ggggTTTTTTCTGTATCTTC‑CGCTCCGCACTCAAGCCTGAAAAAAAATAGGCTGTAGGGC‑TAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                  <  2:308323/151‑1 (MQ=255)
                                                                                      ttttttCTTTATCTTC‑CGCACCGCACAAAACCCTGACAAAAAATATGCTGTAAGGC‑TAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcct                                                                <  1:899993/147‑2 (MQ=255)
                                                                                      tCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                                >  2:67111/1‑149 (MQ=255)
                                                                                      tCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                              >  2:815625/1‑151 (MQ=255)
                                                                                        ttttCTGTATTTTA‑CGCATCCCACTCAAGCCTGAAAGAAAATATGCTGAAGGGC‑TCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                             <  2:249887/150‑1 (MQ=255)
                                                                                         tttttttttCTTA‑CGCACCCCACTCAACCCTGAAAAAAAATATGCTTTAAGGC‑TAAAAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                            <  2:678854/142‑1 (MQ=255)
                                                                                          ttttttATCTTC‑CGCACCCCACTCAACCCAAAAAAAAATTTTGCTGTAAGCC‑TAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                           <  2:794529/145‑1 (MQ=255)
                                                                                          ttCTGTATCTTC‑CGCATCGCACTAAAGCCTGAAAAAAAATTTGCTGTAAGGC‑TCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                           <  2:155041/150‑1 (MQ=255)
                                                                                              gTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                                       >  2:1249956/1‑150 (MQ=255)
                                                                                               tATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                      >  1:127535/1‑150 (MQ=255)
                                                                                                aTCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                    >  1:119857/1‑151 (MQ=255)
                                                                                                 tCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCACAAct                                                   >  2:763505/1‑151 (MQ=255)
                                                                                                 tCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                                   >  2:400716/1‑151 (MQ=255)
                                                                                                     a‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                               >  2:596451/1‑151 (MQ=255)
                                                                                                           tCGCACCCAAGCCGAAAAAAAAAGATGCTGAAAGCC‑TAATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                                                  <  1:689828/143‑1 (MQ=255)
                                                                                                            cccaccaaagctcgACAAAAAATGTGCTGGAAGCC‑CCAAATAAAACGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                            <  1:1070163/135‑1 (MQ=255)
                                                                                                             gCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                         >  2:921439/1‑150 (MQ=255)
                                                                                                                  ccaGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                   >  1:134515/3‑151 (MQ=255)
                                                                                                                   agcccTGAAAAAAAATTGGCTGTAGGCC‑CAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                                                <  1:800476/134‑1 (MQ=255)
                                                                                                                    aGCCTTACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc                                 >  2:1142805/1‑151 (MQ=255)
                                                                                                                      cctaaaaaaaaaTTTGCTGTAAGGC‑TCAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                               <  2:390981/143‑1 (MQ=255)
                                                                                                                       cTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                               >  1:1143444/1‑150 (MQ=255)
                                                                                                                         gACAAAAAATGTGCTATAAGCC‑TCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                              <  1:830186/149‑1 (MQ=255)
                                                                                                                          aaaaaaaaTATCCTGTAAGCC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtta                              <  2:1118447/144‑1 (MQ=255)
                                                                                                                          aCAAAAAATAGGCTGAAGGCC‑AAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                             <  1:737879/149‑1 (MQ=255)
                                                                                                                           cAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                          >  2:827430/1‑151 (MQ=255)
                                                                                                                            aGAAAAGATGCTGTAGGGC‑TGATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                          <  2:594594/150‑1 (MQ=255)
                                                                                                                              aaatTATCCGGTAAGGC‑AAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                       <  2:831750/151‑1 (MQ=255)
                                                                                                                               aaaTATGTTGTAAGGC‑CCATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcatc                      <  1:1210064/151‑1 (MQ=255)
                                                                                                                                    tGCTGAAAGGC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                   <  2:1076625/149‑1 (MQ=255)
                                                                                                                                     gCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                 >  2:1056621/1‑150 (MQ=255)
                                                                                                                                       tGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa                >  1:603303/1‑149 (MQ=255)
                                                                                                                                          aaGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg           >  2:1078230/1‑151 (MQ=255)
                                                                                                                                          aaGGC‑TAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt            <  2:827288/150‑1 (MQ=255)
                                                                                                                                           aGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg          >  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)
                                                                                                                                                      |                                                                                                                                                    
AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTA‑CGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGC‑TCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGC  >  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: ‑