Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,038,7410TA68.9% 62.2 / 94.5 91noncoding (1223/2905 nt)rrlA23S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base T (7/21);  new base A (32/30);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 2.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

AGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038591‑4038867
                                                                                                                                                       |                                                                                                                              
agcCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCtgt                                                                                                                                >  1:540190/1‑151 (MQ=39)
    atcatTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAgg                                                                                                                            <  1:1382/151‑1 (MQ=39)
       ctttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGc                                                                                                                       >  2:161482/2‑151 (MQ=2)
           aaagaaagCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTg                                                                                                                     <  1:225074/151‑1 (MQ=39)
             agaaagCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAGCGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAgg                                                                                                                 >  1:356130/1‑151 (MQ=2)
                    gTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCa                                                                                                            <  1:443090/151‑1 (MQ=39)
                      aaTAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAg                                                                                                           <  1:30973/150‑1 (MQ=40)
                         aGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGt                                                                                                       <  1:112661/151‑1 (MQ=39)
                          gCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCg                                                                                                    >  2:363266/1‑151 (MQ=2)
                           cTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑AGGTAGAGGAGCGTTCTGTAAGCCTATGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCg                                                                                                    >  2:146204/1‑150 (MQ=2)
                            tCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCg                                                                                                    >  2:332195/1‑151 (MQ=39)
                                 ggTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGc                                                                                               <  1:89429/151‑1 (MQ=39)
                                 agTCGAGTCGGCCTGCGCGGAAGATGTAGCGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                                                                                             <  2:356088/150‑1 (MQ=2)
                                  gTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                                                                                             >  2:191977/1‑150 (MQ=2)
                                    cGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                                                                           <  1:363309/150‑1 (MQ=2)
                                      aGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                                                                           >  1:620189/1‑150 (MQ=40)
                                         cGGTCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGCGCGAATGCTGACATaa                                                                                       <  2:733605/151‑1 (MQ=255)
                                         cGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGt                                                                                     >  2:674196/1‑151 (MQ=2)
                                            ccTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTa                                                                                    <  2:488565/151‑1 (MQ=39)
                                               gcgcGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGac                                                                               <  1:161499/151‑2 (MQ=2)
                                               gcgcGGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑AGGTAGAGGAGCGTTCTGTAAGCCTATGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 <  1:146220/149‑1 (MQ=2)
                                                cgcgGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGac                                                                               >  2:674088/1‑149 (MQ=2)
                                                 gcgGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAActaca                                                                              >  1:449907/1‑146 (MQ=2)
                                                  cgGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 >  2:353698/1‑146 (MQ=2)
                                                  cgGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 <  1:353740/146‑1 (MQ=2)
                                                  cgGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGAt                                                                               >  2:481149/1‑150 (MQ=40)
                                                  cgGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATa                                                                              <  1:473115/151‑1 (MQ=39)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCt                                                                                              <  2:240344/132‑1 (MQ=2)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCt                                                                                              >  1:240375/1‑132 (MQ=2)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                                                                           >  1:383876/1‑135 (MQ=2)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                                                                                           <  2:383833/135‑1 (MQ=2)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACaatg                                                                           >  2:332836/1‑149 (MQ=2)
                                                   ggAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACaatg                                                                           <  1:332876/151‑3 (MQ=2)
                                                    gAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAActacaatgg                                                                          <  2:449856/151‑9 (MQ=2)
                                                     aaGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGTAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAATCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACaatggg                                                                         <  2:217967/151‑5 (MQ=2)
                                                      aGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGc                                                                          <  1:421917/151‑1 (MQ=39)
                                                       gATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACaatgggtg                                                                       <  1:674275/151‑7 (MQ=2)
                                                        aTGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCgg                                                                        <  2:155426/151‑1 (MQ=39)
                                                            aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                                                     <  1:756852/150‑1 (MQ=40)
                                                            aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                                                     <  2:784492/150‑1 (MQ=40)
                                                               ggggCTAAATCATGCCCCGAAGCTGCGGCAGCGACACTATGGGTTGTTGGGTAGGGGAGCGTTCTGTAAGGCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTATCGATAAAGCGGGTGaaaa                                                                 <  2:333066/151‑1 (MQ=255)
                                                                         ccATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCg                                                        <  2:542539/151‑1 (MQ=39)
                                                                             gCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCgg                                                    <  2:387085/151‑1 (MQ=39)
                                                                                  gAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAc                                               >  2:586312/1‑151 (MQ=39)
                                                                                         cGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGt                                        >  2:268130/1‑151 (MQ=39)
                                                                                            cAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTcc                                     <  1:33674/151‑1 (MQ=39)
                                                                                                     cgtaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgc                                          <  2:451022/125‑2 (MQ=2)
                                                                                                     cgtaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgc                                          >  1:451074/13‑136 (MQ=2)
                                                                                                     cgcaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  1:331104/125‑1 (MQ=2)
                                                                                                     cgcaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          >  2:331064/13‑137 (MQ=2)
                                                                                                     cgcaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACg                            >  1:522501/13‑151 (MQ=2)
                                                                                                        aagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 >  1:297712/10‑95 (MQ=11)
                                                                                                        aagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 <  2:297675/86‑1 (MQ=11)
                                                                                                          gttgttGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTaa                        <  2:31878/150‑1 (MQ=31)
                                                                                                            gacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCg                     >  2:649077/6‑150 (MQ=2)
                                                                                                              ttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                   >  1:694417/1‑151 (MQ=31)
                                                                                                                 ggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACgcgc                                                                              <  1:292966/89‑4 (MQ=11)
                                                                                                                 ggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACgcgc                                                                              >  2:292929/1‑86 (MQ=11)
                                                                                                                 ggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCgacgc                 >  1:565988/1‑146 (MQ=2)
                                                                                                                  gTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 >  2:320888/1‑117 (MQ=2)
                                                                                                                  gTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 <  1:320927/117‑1 (MQ=2)
                                                                                                                   tAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCgacg                  >  2:31275/1‑144 (MQ=2)
                                                                                                                   tAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCgacg                  <  1:31279/147‑4 (MQ=2)
                                                                                                                    aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg             <  1:341675/151‑1 (MQ=31)
                                                                                                                      aggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTgcgc                                >  2:347037/2‑127 (MQ=11)
                                                                                                                      aggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTgcgc                                <  1:347078/129‑4 (MQ=11)
                                                                                                                       ggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGaaaaacacccacg                                                        >  1:244113/1‑96 (MQ=11)
                                                                                                                       ggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGaaaaacacccacg                                                        <  2:244082/105‑10 (MQ=11)
                                                                                                                        gAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAATACCCACGCCGAAAGACCAAg                                           >  1:6396/1‑117 (MQ=2)
                                                                                                                        gAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAATACCCACGCCGAAAGACCAAg                                           <  2:6396/117‑1 (MQ=2)
                                                                                                                        gAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAg                                           <  1:172942/117‑1 (MQ=2)
                                                                                                                        gAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAg                                           >  2:172924/1‑117 (MQ=2)
                                                                                                                         aGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGt        <  1:331710/151‑1 (MQ=31)
                                                                                                                          gCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTc       <  1:666725/151‑1 (MQ=31)
                                                                                                                          gCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCaa                                            >  2:369618/1‑114 (MQ=2)
                                                                                                                          gCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCaa                                            <  1:369661/114‑1 (MQ=2)
                                                                                                                             ttCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAc    >  1:20630/1‑151 (MQ=31)
                                                                                                                               cTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggttc  >  1:368788/1‑147 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  2:192672/109‑1 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          >  1:192692/1‑109 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGttggttt  <  2:565923/149‑7 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggttc  <  2:522440/149‑5 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCAGAAAGACCAAg                                           >  2:337243/1‑108 (MQ=2)
                                                                                                                                 gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCAGAAAGACCAAg                                           <  1:337283/108‑1 (MQ=2)
                                                                                                                                    aGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCaa                                            <  2:103726/104‑1 (MQ=2)
                                                                                                                                    aGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCaa                                            >  1:103738/1‑104 (MQ=2)
                                                                                                                                       cTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  1:578779/103‑1 (MQ=2)
                                                                                                                                       cTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          >  2:578712/1‑103 (MQ=2)
                                                                                                                                         gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 >  2:393364/1‑94 (MQ=2)
                                                                                                                                         gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 <  1:393407/94‑1 (MQ=2)
                                                                                                                                         gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggt    >  1:342162/1‑137 (MQ=2)
                                                                                                                                         gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggt    <  2:342122/139‑3 (MQ=2)
                                                                                                                                                       |                                                                                                                              
AGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038591‑4038867

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