Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,311 G→T R81S (CGT→AGT)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,3110GT94.9% 151.0 / ‑3.6 63R81S (CGT→AGT) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (21/35);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTG  >  NC_000913/2867066‑2867557
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
ggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCCCCCAGCGCGCGACtcgcc                                                                                                                                                                                                                                                     >  1:641209/1‑249 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATaa                                                                                                                                                                                                                                               <  1:2549/251‑1 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATa                                                                                                                                                                                                                                                <  1:226887/250‑1 (MQ=255)
         cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAActt                                                                                                                                                                                                                                           <  1:780957/251‑1 (MQ=255)
             aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAActtcttc                                                                                                                                                                                                                                       <  1:393566/251‑1 (MQ=255)
                 cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCTTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCg                                                                                                                                                                                                                                  <  2:605619/251‑1 (MQ=255)
                     ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAATCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACTCGCAAAAAAAACTTCTTCTTCGGCCg                                                                                                                                                                                                                              >  1:11731/1‑251 (MQ=255)
                         aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTACACCAGACGCAAGTTACTCTCGATCCTCCGGCGGCGGGAGGGGACATCTCCAAGCCGTTCGCGGCTCGCAAAAAAAACATCTTCTTCGGACGTTa                                                                                                                                                                                                                           >  1:340150/1‑250 (MQ=255)
                                gtcGAATCCACCAGGTTGCGGATGTTGAGAAGCGGAAACCACTTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACTCCATACCACGATTTCCATAACGGCGGGCAATTTTTACCTCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACCgtggt                                                                                                                                                                                                                   <  1:477941/249‑1 (MQ=255)
                                 gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtggtg                                                                                                                                                                                                                  <  1:107141/251‑1 (MQ=255)
                                     aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATa                                                                                                                                                                                                              <  2:430009/251‑1 (MQ=255)
                                      aacaccaAGTTGCGTATGTAGAGAAGCGGAAACCACGTTCCAGGACAAAAATCTCAACCGCGCGGAGCAGCCCCAGGTTGCCCTCTTCGAAAAGGTCCAGCAACGCCAGACCACGATAGCCATAACGGCGGGCAATTTTTACCACCAGACGAAAGTTACTCTCGATCAACCGGCGGCGAGAGGCGAAAAATCCACGAAGTGCGCGAAACACAAAAAAAACTTCTTCAGCGGCAATTAACAGTGGTGAATaa                                                                                                                                                                                                             <  2:400216/249‑1 (MQ=255)
                                           caggtgtcGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCACTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACTGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAt                                                                                                                                                                                                        <  2:198233/244‑1 (MQ=255)
                                           caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCGCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAt                                                                                                                                                                                                        <  1:375936/251‑1 (MQ=255)
                                            aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCCATTTTTTACACCCGACGCAAGTTAAACTCCGTCATCCCGCGGCGAGAGGCGACCTCTTCCCGCCATGCGGGCACCGCAAAATAAAATTCTTTCTCGGCCGCTAACAGAGGAGCATAACCaaac                                                                                                                                                                                                       >  1:145040/1‑249 (MQ=255)
                                                tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGCATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCCGCGGCGGGAGGCGACATCCCCCCGCAGGGCGCGACTCGCCAAATAAAATTCTTCTTCTGCCGGTTACAGTTGTGAATAACCAAACTCAc                                                                                                                                                                                                   >  1:9691/1‑251 (MQ=255)
                                                                cGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCCAGTTACTCTTGATTATCCGCCGGCGGGAGGCGGCCTCTCCCCCCCGGGCGAGACTCGCAAAATAAACATTTTCTTTGGCCGTTTACAGGGGTGAATAACCAATCTCACCACGGGAAAGCTGAgg                                                                                                                                                                                   >  1:334002/1‑250 (MQ=255)
                                                                 ggAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCCTCCGGCGGCGAGAGGCGACATCTCCCCCCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCTGCCGTTAACAGTGGTGAATAACCACTCTCACCAAGGTAAAACTGAGTc                                                                                                                                                                                  >  1:468990/1‑251 (MQ=255)
                                                                 ggAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGAACTTTTTACCACCAGACGCAAGGTACTCTCGCTCATCCGGCGGCGAGGGGCGACAACTCCCACCAGTGAGCGACTCGCAAAATAAACTTCTTCTTCGGCCGAGAACAGTGGGGAATAAACAATCTCAcccaagaaaaactgaggc                                                                                                                                                                                  >  1:506131/1‑235 (MQ=255)
                                                                  gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAAGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCCGCCGTTAACAGTGGGGAATAACCAATCTCACCCAGGGAAAGCTGAGTcg                                                                                                                                                                                 >  1:721723/1‑251 (MQ=255)
                                                                   aaaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAACCGCCAGTTACTCTCTATCATCCCGCGGCGGGAGGCGCCCACTCCCCACACTGCGCGCCACGCAAAAAAAACTTCTTCTTTGCGCGTTAAACGTGGGGAATAACCAAtactcccccaggaaaacggtgtgcc                                                                                                                                                                                >  2:807675/1‑227 (MQ=255)
                                                                       cACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTTCCCTCTTCGATAAGGTCCAGCAACGCCAGTCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTcc                                                                                                                                                                            <  2:515034/251‑1 (MQ=255)
                                                                         cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCGATTTTTTCCCCCCGACACCAGTTACTTTCGGTCATCCGGCGGCGGGAGGCGGCATCTCCCCCCAGTGCTCGACTCCAAAAATAAACTTCTTCTTCTGCCGTTAACACTGGTGAGTAACCAATCTCACCaacagaaaactgagtagtgaccc                                                                                                                                                                           >  2:402053/1‑229 (MQ=255)
                                                                             ccGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGGCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacaca                                                                                                                                                                      <  2:281113/251‑1 (MQ=255)
                                                                                   cAAACTACTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAGTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCCATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                                                                                                <  2:1854/251‑1 (MQ=255)
                                                                                               ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCCATTTTTATCACCCAACGCAAGTTACACTCGATCATCCGGCGGCGAGAGGCGCCATCTCCACACCCGGCCCGCCCCCCCAAATAAACATCTTTTTCGGCCGTTTACAGTGGGGAATGAACAATCTCACCCAGGCTAAGCTGTGTTGCGTCCacaccaccccatgtggctcccctc                                                                                                                                                    >  2:668101/1‑228 (MQ=255)
                                                                                                    cgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAATTTTAACACCCGAAACCAATTACTCTCGGACATCCGGCGGCGGGAGGGGACAGCTCCCCGGCCATCGCGACTCGCAAAAAAAACATCCTCTTCGGACGCTAAAAGGGGGGGAAAACAAATCTCAACAAGGGAAAGCTGAGGCGCGTCCAaaacacaccggggggccacccgcgata                                                                                                                                                >  2:715506/1‑229 (MQ=255)
                                                                                                        tCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAACCCCCGACGCAAGATACTCTCGATCATCCGGCGGCGAGAGGCGAACTCTCCACCCCATGCGCGACTCGCAAAAAAAACTTCTTCTTCTGCCCGTTACAGTGGGGGATAAACAATCTCACCAAGGTAAAGCTTAGTCGCGTCCCACCCCCCCTGTGTTGCTCCCTCCGATAAccgt                                                                                                                                           >  1:728041/1‑248 (MQ=255)
                                                                                                                       ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaaa                                                                                                                            <  1:48213/251‑1 (MQ=255)
                                                                                                                          ctTCTATACGGTCCAGCCACTCCAGACCACGATTGCCACAACGGCCGTCAATTTTTACCAACAGACGCAAGTTACTCTCCATCATCCGGCGGCGAGAGGCGTCATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAATCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACTCGCTGTGTGGCACCCTGCGATAACAGTTCCTCTTCGTCCAAAt                                                                                                                           <  1:680874/249‑1 (MQ=255)
                                                                                                                          ctTCGATAAGGTCCAGCAACGCCAGGCCACGATTGCCATAACGGCTGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCACAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGATTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAAt                                                                                                                           <  1:12284/249‑1 (MQ=255)
                                                                                                                            tCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                                                                                        <  1:703654/250‑1 (MQ=255)
                                                                                                                                        atcAACGCCAGACCACGATTGCCATAACGGCGGTCAATTTTTACCACCAGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAAGCAGTGCGCGACTCGCAAAATAACCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTTTGGCGCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTc                                                                                                           <  2:88510/249‑1 (MQ=255)
                                                                                                                                             cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGtt                                                                                                      <  1:223841/251‑1 (MQ=255)
                                                                                                                                              gCCAGACCACGATGGCCACAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAGCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                                                                                                     <  2:721723/251‑1 (MQ=255)
                                                                                                                                                      ccGATTGCCATAACGGCGGGCAATTTTTACCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACTCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                                                                             <  2:83862/250‑1 (MQ=255)
                                                                                                                                                           tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTCCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCtt                                                                                         <  1:426667/250‑1 (MQ=255)
                                                                                                                                                            gcaagaaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTACGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTGATCACTGGGTTCCTGTTCTACTAAGGCCtttt                                                                                       <  2:778820/246‑1 (MQ=255)
                                                                                                                                                                         gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                          <  1:428022/251‑1 (MQ=255)
                                                                                                                                                                                ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCg                                                                   >  1:317622/1‑251 (MQ=255)
                                                                                                                                                                                  accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                 >  1:314060/1‑251 (MQ=255)
                                                                                                                                                                                  accaccAGACGCAAGTTACCCGCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTTCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                 <  2:340150/251‑1 (MQ=255)
                                                                                                                                                                                          aCGCAAGTTACGCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                         <  1:310066/251‑1 (MQ=255)
                                                                                                                                                                                            gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                                                       >  1:290862/1‑251 (MQ=255)
                                                                                                                                                                                             caatttACTCTCGACCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                                                      <  2:668840/247‑1 (MQ=255)
                                                                                                                                                                                               aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc                                                    <  1:789590/251‑1 (MQ=255)
                                                                                                                                                                                                gtgaCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg                                                   <  1:777972/248‑1 (MQ=255)
                                                                                                                                                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTCTCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                >  2:715206/1‑251 (MQ=255)
                                                                                                                                                                                                   aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                <  1:487086/251‑1 (MQ=255)
                                                                                                                                                                                                    cgtccGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc                                               <  2:468990/247‑1 (MQ=255)
                                                                                                                                                                                                     tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                              >  1:462619/1‑251 (MQ=255)
                                                                                                                                                                                                        cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGGCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt                                             <  1:402053/249‑1 (MQ=255)
                                                                                                                                                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                                                                                                 >  1:63354/1‑148 (MQ=255)
                                                                                                                                                                                                         gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTCCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt                                                                                                                                                 <  2:63354/148‑1 (MQ=255)
                                                                                                                                                                                                                cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCCTTTAAATc                                   >  2:312252/1‑251 (MQ=255)
                                                                                                                                                                                                                   cggcgAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                >  1:548595/1‑251 (MQ=255)
                                                                                                                                                                                                                          agGCGACATCTCCACGCCGTGCGCGACTTGCAAAATAAACTTCTGCTACGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAGACCGCATCTTCATTTAAATCAGGAACTTTc                         <  2:145040/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGt                                                                                                             >  1:638502/1‑164 (MQ=255)
                                                                                                                                                                                                                             cGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGt                                                                                                             <  2:638502/164‑1 (MQ=255)
                                                                                                                                                                                                                                     ccACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAt                  <  1:174595/247‑1 (MQ=255)
                                                                                                                                                                                                                                     ccACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAGCTTCATTTAAATCATGAACTTTCAGCGTAt                  <  1:631694/247‑1 (MQ=255)
                                                                                                                                                                                                                                            cAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                                                                                      >  2:545248/1‑78 (MQ=255)
                                                                                                                                                                                                                                            cAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTg                                                                                                                                                                                      <  1:545248/78‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg  >  1:789454/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTG  >  NC_000913/2867066‑2867557

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