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.1% 42.8 / NA 18R81S (CGT→AGT) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (5/11);  total (6/12)
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.36e-01

GGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTC  >  NC_000913/2867065‑2867534
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
gggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCa                                                                                                                                                                                                                             <  1:196262/251‑1 (MQ=255)
              aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAActtcttc                                                                                                                                                                                                               <  1:292697/251‑1 (MQ=255)
                 gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCCAACGCAAATTACTCTCGCTCATCCGGCGGCGAGAGGCGACAACTCCACGCAGGGGGGGGCACCCAAAAAAAActtcttcttt                                                                                                                                                                                                            >  1:197717/1‑250 (MQ=255)
                 ccccGTTCAACCGTCGTGCGAAGCCCCCAGGTTGCGAATGTGTAGACGTGGCAACCACGCTCCGGGGCAAACATCTCTGCCGCGCGGATCAGCCCCAGGTTGCCCTCTTCTATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGAAAGTTACTCGCGATCATCCGGCGGCGAGAGGCGACATCGTCACGCAGTGCGCGACTCGCAAAATAAACttcttcttc                                                                                                                                                                                                            <  2:143969/250‑1 (MQ=255)
                                     aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACGCAAGATACTCTCGCTCATCCCGCGGCGAGAGGCGACATCCCCACCCAGGGCGCGGCTCGCCAAATAAActtcttctttatcggaaaagcacacatccg                                                                                                                                                                                        >  1:286841/1‑230 (MQ=255)
                                     aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCAGGGTCAAACTTCTCTACCGCGCGGACCAGCCCCAGGTTGCACCCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAActtcttct                                                                                                                                                                                                              <  2:286841/229‑1 (MQ=255)
                                         caccaGGTTGTGTATGTTTAGCAGCGGAAACCACGTTCCGGGTCAAACCTCTCTACCGCGAGGATCAGCCCCAGGTTGCCCTCTTCGAAAAGTTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAGTTTTACCACCAGACTCAAGTTACTCTCGATCATCCGGCGGCGAGATGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTTAATAAcc                                                                                                                                                                                    <  2:123915/251‑1 (MQ=255)
                                                      aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCCGAAGCAAGTTACTTTCGATCATCCGGCGGCGCGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAAATTCTTCTTCCGCCGTTTAAAGTGGTGAATAACCAATCTCACCAAgg                                                                                                                                                                       >  1:59691/1‑251 (MQ=255)
                                                                gcgtaaaccacgtaacgtgataagcTTCTCTCCAGCGCGGATCACCCCCAGATTGCCCTCTTCGATAAGGTCAAGCAACGCCAGACCACGATTGTGATAGCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGTCGGCGAGAGGTGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTGCTTCTGCCGTTAACAGTGGTGAATAACTAATCTCACCAAGGTAAAGCTGAg                                                                                                                                                             <  2:175881/227‑1 (MQ=255)
                                                                           gttcccggcacaaCTTCCCTACCGCGCGGATCAGCCCCTGGTTGCCCGCTTCGATAAGGTCTAGCAACGCCAGACCTCGATTGCCATAGCGGCGGGCAATTTTTACCACCATACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGTCGCGCAAAAGAAACTTCTTCTGCGGCTGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTATAGCTGAGTCGCGACCAac                                                                                                                                                  <  2:266291/240‑1 (MQ=255)
                                                                                               aCCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCt                                                                                                                               <  1:168604/250‑1 (MQ=255)
                                                                                               aCCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg                                                                                                                              <  1:138050/251‑1 (MQ=255)
                                                                                                                        cccctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTGCGGCCGTTAACAGTGTTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaaa                                                                                                     <  1:182022/248‑1 (MQ=255)
                                                                                                                                                       ccGATTGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTGCTTCTTCGGCCTTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTTAGTCGCGTCCAACACACGCTGTGTGGCTTCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCTCTGGGTTCCTGTTCTACTAAgg                                                                      <  2:288829/250‑1 (MQ=255)
                                                                                                                                                                          gCAATTTTTACCACCAGACGCCAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACAAGCTGGGTGGCTCCCTGCGATACCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTTCTAAGGCCTTTTCGTCATAAACCt                                                    >  2:20943/1‑250 (MQ=255)
                                                                                                                                                                                            cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaaa                                 <  1:18360/251‑1 (MQ=255)
                                                                                                                                                                                                           atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTAATAAGGCCTTTTCGTCAAAAACCTAAACTCCGTTTTCATCAAATTCCGCATCTTCAtt                  >  2:93478/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggCGACATCTCCACGCAGTGCGCGACTCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACGTTc  >  1:259785/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
GGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTC  >  NC_000913/2867065‑2867534

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