Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,213,4340GA35.0% 61.4 / 53.5 60V130I (GTT→ATT) rpoDRNA polymerase sigma factor RpoD
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (10/11);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACT  >  NC_000913/3213167‑3213717
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                            
atcGAAGACATCATCCAAATGATCAACGACATGGTCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCATGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGaa                                                                                                                                                                                                                                                                          >  1:616457/1‑288 (MQ=255)
                 aaTGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATc                                                                                                                                                                                                                                                           >  1:22742/1‑286 (MQ=255)
                              aTGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGaa                                                                                                                                                                                                                                            <  1:271321/288‑1 (MQ=255)
                                   agggaggggaagggAAGGAGAACGGGATGTCGTTGAACTGAAGAGTGCGTAAAGAACAGCGGATCAGGAAGCGGCCGGAGGCGCTGCACAGGTGCTTTCCAGCTTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGGTGACCCGCGAAGGCGAAATTGACATCGATAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTa                                                                                                                                                                                                                                       <  2:776727/276‑1 (MQ=255)
                                           tGATAGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTa                                                                                                                                                                                                                                                      <  1:119008/265‑1 (MQ=255)
                                                  agaagCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAg                                                                                                                                                                                                                         <  1:649202/287‑1 (MQ=255)
                                                           ggATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGaagaag                                                                                                                                                                                                                >  1:229215/1‑287 (MQ=255)
                                                                           aTGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGGTGACCCGCTAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACCGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTg                                                                                                                                                                                               >  1:752274/1‑288 (MQ=255)
                                                                                 gctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCTAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACg                                                                                                                                                                                                                                     >  1:339337/1‑244 (MQ=255)
                                                                                 gctgAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCACCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCGTGTACATGCGTGAAATGGTCTCCGTTGAACTGTTGACCCGCGAAGGCGAAAGTGACATCGCTAAGCGTCTTTAACCCGGGACCAACCAGGt                                                                                                                                                                                                                                                                                             >  2:192360/1‑188 (MQ=255)
                                                                                           ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGa                                                                                                                                                                                                                                                                           >  1:309201/1‑196 (MQ=255)
                                                                                           ccGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGa                                                                                                                                                                                                                                                                           <  2:309201/196‑1 (MQ=255)
                                                                                            cgcgGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTAACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGAACACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTttgttg                                                                                                                                                                               >  1:515034/1‑287 (MQ=255)
                                                                                               ggACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGAccc                                                                                                                                                                           >  1:841012/1‑288 (MQ=255)
                                                                                                          ctgcCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGa                                                                                                                                                                                                                                                                           <  1:589638/181‑1 (MQ=255)
                                                                                                           tgcCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCACTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCagaaga                                                                                                                                                               <  1:831524/288‑1 (MQ=255)
                                                                                                            gcCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGCACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCGGAACGCagaaga                                                                                                                                                               <  1:528667/287‑1 (MQ=255)
                                                                                                                                      cAGCGTGGAATTTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGt                                                                                                                                    <  1:175434/288‑1 (MQ=255)
                                                                                                                                            ggAATCTGAAATCGGGCCTACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGTTCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGtt                                                                                                                               <  1:326997/287‑1 (MQ=255)
                                                                                                                                              aaTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGAt                                                                                                                                                                                                                                                            >  1:277957/1‑160 (MQ=255)
                                                                                                                                              aaTCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGAt                                                                                                                                                                                                                                                            <  2:277957/160‑1 (MQ=255)
                                                                                                                                                                     gACCCGGTACGCATGTACATGCGTGAAATGGGCACGGTTGAACTGTTTACTCGCGGAGGCGTAATTGACATCGCTAAGCATATCGAAGACGGGATCAGCCAGGTACAATGCTCCGTTGCTGTATATCCGGAAGCGATCACCTGTCTGCTGGAACAGTATGATCGTGTTGAAGCAGAAGACGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGAACCGAACGCAGAAGAAGATCTGGc                                                                                                                                                      >  1:501951/1‑239 (MQ=255)
                                                                                                                                                                      aCCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGaagaag                                                                                                                                                                                                                >  2:208767/1‑180 (MQ=255)
                                                                                                                                                                      aCCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGaagaag                                                                                                                                                                                                                <  1:208767/180‑1 (MQ=255)
                                                                                                                                                                       cccGGTACGCATGTACATGCGTGAAATGGGCACCGTTGATCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCa                                                                                                                                           <  1:128366/248‑1 (MQ=255)
                                                                                                                                                                           gTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGCAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCgg                                                                                                                                                                                                                                                                   >  1:814113/1‑124 (MQ=255)
                                                                                                                                                                           gTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATTCGGAAGCGATCATCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAATAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCTGAACGCagaagtag                                                                                                                                                             >  1:244057/1‑226 (MQ=255)
                                                                                                                                                                           gTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTa                                                                                                                                                                                                                                                      >  1:627003/1‑137 (MQ=255)
                                                                                                                                                                           gTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTa                                                                                                                                                                                                                                                      <  2:627003/137‑1 (MQ=255)
                                                                                                                                                                           gTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCgg                                                                                                                                                                                                                                                                   <  2:814113/124‑1 (MQ=255)
                                                                                                                                                                              agaagagtaaggCGGGAAATGGACCCCGGTTAAACGGTTGACCCCAGAAGGCGAAATTGACTTCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACg                                                                                              <  2:359344/276‑1 (MQ=255)
                                                                                                                                                                                      cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGTCATCGCTAAGCGTATTGAAGACGGGATCGACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTa                                                                                                                                                                                                                                       <  1:338521/141‑1 (MQ=255)
                                                                                                                                                                                      cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTa                                                                                                                                                                                                                                       >  2:338521/1‑141 (MQ=255)
                                                                                                                                                                                      cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGAt                                                                                                                                                                                                                                                            >  1:16501/1‑120 (MQ=255)
                                                                                                                                                                                      cATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGAt                                                                                                                                                                                                                                                            <  2:16501/120‑1 (MQ=255)
                                                                                                                                                                                          cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTgg                                                                                                                                                                                                                                              <  2:335994/130‑1 (MQ=255)
                                                                                                                                                                                          cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTgg                                                                                                                                                                                                                                              >  1:335994/1‑130 (MQ=255)
                                                                                                                                                                                          cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGcgcg                                                                                                                                                                                                            <  1:547876/164‑1 (MQ=255)
                                                                                                                                                                                          cGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTACCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGcgcg                                                                                                                                                                                                            >  2:547876/1‑164 (MQ=255)
                                                                                                                                                                                             gAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAg                                                                                     >  1:359344/1‑280 (MQ=255)
                                                                                                                                                                                              aaaTGGGCACCGTTGAACTGTTGACCCGCGAAGGCGTAATCGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACTTCGGTTCTGAGCTTTCCCAGGAAGACCTGGACGATGACGAATATGAAAACgaagaa                                                                            >  1:564009/1‑287 (MQ=255)
                                                                                                                                                                                                 tGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAgg                                                                                                                >  1:547893/1‑249 (MQ=255)
                                                                                                                                                                                                 tGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGGCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAgg                                                                                                                <  2:547893/249‑1 (MQ=255)
                                                                                                                                                                                                              aaCTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTGGTTGACCCGAACGCAGAAGAAGAACTGGCACCTACCGCCACCCACGTCGGTTCTGAGTCTTGCCAGGAAGTACCGGACGATGACGGACATGAAGTAGGAGGAAGAGGGGGAGgacgag                                                             >  2:317281/1‑287 (MQ=255)
                                                                                                                                                                                                                     gACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCga                                                     >  1:13732/1‑288 (MQ=255)
                                                                                                                                                                                                                         cgcgAAGGCGAAATTGACATCGCTAAGCGGATTGAGAACAGGATAAACCAGATTAAATGCTCCGTTGCAGAATACACGGAAACGATCACCCATCTGCTGAAACAGTACGATCGTGTGGAAGAAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCGACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCgatgat                                                 <  1:330650/288‑1 (MQ=255)
                                                                                                                                                                                                                           cgAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCgtgt                                                                                                                                                                                                                              >  1:659245/1‑114 (MQ=255)
                                                                                                                                                                                                                           cgAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCgtgt                                                                                                                                                                                                                              <  2:659245/114‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGACAACGCTAAGCGTATTGAAGACGGGATCAACCAGATTAAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGTTCGTGTTGAAGAAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCAt                                      <  1:272226/286‑1 (MQ=255)
                                                                                                                                                                                                                                         acacCGCAAAGCGTATTGAAGACGGAATCACCCAGATTCAATGCTACGTTGCTGAAGATCCGGAAGCGATCACCTATCTGCTGAAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATATGGACGAGGACGAAGATGAAGACAAAGAAGATGGCGATGCCGACAGCGCCGATGATGACAACAGCATCGAcc                                 <  1:232133/288‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTttg                                                                                                                                                                                  >  2:493705/1‑139 (MQ=255)
                                                                                                                                                                                                                                              gCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTttg                                                                                                                                                                                  <  1:493705/139‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCg                                                                                                                                                                          >  1:593810/1‑137 (MQ=255)
                                                                                                                                                                                                                                                        ttGAAGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCg                                                                                                                                                                          <  2:593810/137‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGACGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCAt                                      <  1:230583/264‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCGGGATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGAACGTGTTGAAGCAGAAGAAGCACGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGc                                        >  1:135909/1‑260 (MQ=255)
                                                                                                                                                                                                                                                                 ggATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTc                                                                                                                     <  1:215027/181‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ggATCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTc                                                                                                                     >  2:215027/1‑181 (MQ=255)
                                                                                                                                                                                                                                                                   aTCAACCAGATTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTTTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCCCGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATTACGACAGCGCCGAGGATGACAACAGCAt                                      <  1:475393/257‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ccATGTTTAATGCTTTGTTTCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATATGATCACCGGCTTTGTTGACCCGAACGAAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAATAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGAACCGGAACTGGCTCGCGAAAAATTTGCGGAACt  <  1:513730/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                            
ATCGAAGACATCATCCAAATGATCAACGACATGGGCATTCAGGTGATGGAAGAAGCACCGGATGCCGATGATCTGATGCTGGCTGAAAACACCGCGGACGAAGATGCTGCCGAAGCCGCCGCGCAGGTGCTTTCCAGCGTGGAATCTGAAATCGGGCGCACGACTGACCCGGTACGCATGTACATGCGTGAAATGGGCACCGTTGAACTGTTGACCCGCGAAGGCGAAATTGACATCGCTAAGCGTATTGAAGACGGGATCAACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACT  >  NC_000913/3213167‑3213717

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