Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,038,7470AT70.0% 55.2 / 86.9 90noncoding (1229/2905 nt)rrlA23S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base A (7/20);  new base T (33/30);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038598‑4038867
                                                                                                                                                      |                                                                                                                        
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                                                                                 <  1:353740/146‑1 (MQ=2)
                                           cgGAAGATGTAACGGGGCTCAAACCACACACCGAAGCTGCGGGTGTCACGCAAGTGACG‑‑CGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACg                                                                                 >  2:353698/1‑146 (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                                                                     <  2:784492/150‑1 (MQ=40)
                                                     aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTg                                                                     <  1:756852/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                                          >  2:331064/13‑137 (MQ=2)
                                                                                              cgcaagtgacgcggTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  1:331104/125‑1 (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                                                 <  1:320927/117‑1 (MQ=2)
                                                                                                           gTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 >  2:320888/1‑117 (MQ=2)
                                                                                                            tAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCgacg                  <  1:31279/147‑4 (MQ=2)
                                                                                                            tAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCgacg                  >  2:31275/1‑144 (MQ=2)
                                                                                                             aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg             <  1:341675/151‑1 (MQ=31)
                                                                                                               aggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTgcgc                                <  1:347078/129‑4 (MQ=11)
                                                                                                               aggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTgcgc                                >  2:347037/2‑127 (MQ=11)
                                                                                                                ggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGaaaaacacccacg                                                        <  2:244082/105‑10 (MQ=11)
                                                                                                                ggAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGaaaaacacccacg                                                        >  1:244113/1‑96 (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                                          >  1:192692/1‑109 (MQ=2)
                                                                                                                          gTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  2:192672/109‑1 (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                                            >  1:103738/1‑104 (MQ=2)
                                                                                                                             aGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCaa                                            <  2:103726/104‑1 (MQ=2)
                                                                                                                                cTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          >  2:578712/1‑103 (MQ=2)
                                                                                                                                cTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAgg                                          <  1:578779/103‑1 (MQ=2)
                                                                                                                                  gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 <  1:393407/94‑1 (MQ=2)
                                                                                                                                  gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAg                                                 >  2:393364/1‑94 (MQ=2)
                                                                                                                                  gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggt    >  1:342162/1‑137 (MQ=2)
                                                                                                                                  gtgAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATGGGTGTGAAAAACACCCACGCCGAAAGACCAAGGGTTCCTGCGCAACGTTAATCGACGCAGGGTTAGTCggt    <  2:342122/139‑3 (MQ=2)
                                                                                                                                                      |                                                                                                                        
ATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCC  >  NC_000913/4038598‑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: