Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,430 G→A A41V (GCC→GTC)  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,4300GA100.0% 55.7 / NA 21A41V (GCC→GTC) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/8);  total (13/8)

CACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCC  >  NC_000913/2867215‑2867656
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  
naCGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACCCAAGTTACTCACGATCATCCGCCGGCGAGCGGCGCCACCTCCACGCAGTCCGCGCCGCCCAAAATACACTTCTTCTTCCGCCGTTTACACTGGTGAATAACCAATCTCCCCAAGGGAAAGCTGCGTCGCCCCCAACACCCGCTGTGTGGCTCCCCGCGATAAAAGTTACTCTTAGACCAAATCGTTAACACTGGGTTCCTGTTttacta                                                                                                                                                                                                   >  1:92164/2‑249 (MQ=255)
         cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCTATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCg                                                                                                                                                                                        <  1:65839/251‑1 (MQ=255)
                                   aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGAGGAGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGAGTTCCTGTTCTACTAAGGCCTTTTCGTCAAACACCTCAACTCCGTTCtcatc                                                                                                                                                              >  1:453530/1‑251 (MQ=255)
                                                                    gagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGt                                                                                                                                                                                       <  1:475235/193‑1 (MQ=255)
                                                                                          cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCt                                                                                                                                                                                             >  1:469347/1‑165 (MQ=255)
                                                                                                ncgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaa                                                                                                                                                                                    >  1:63365/2‑168 (MQ=255)
                                                                                                gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCGGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACGAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtttn                                                                                                                            <  1:60230/224‑2 (MQ=255)
                                                                                                                     cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACg                                                                             <  1:270147/250‑1 (MQ=255)
                                                                                                                      ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACAAGCTGTGTGGCACCCTGCGATAACAGATCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                        >  1:286624/1‑238 (MQ=255)
                                                                                                                         cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTg                                                                                                                                                                                                                   >  1:14648/1‑112 (MQ=255)
                                                                                                                                 gtggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCAGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                 >  1:341484/1‑250 (MQ=255)
                                                                                                                                      gNNTAACCANNCTCACCAAGGGAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTg                                                                                                                                                                                                                   <  2:14648/99‑1 (MQ=255)
                                                                                                                                                          gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTACCTGCGATAACAGTTCCTCTTCGACCAAAACGTTATCCCTGGGTTAATGTTCTACTAACGCCTTTTCGTCAAAAAACTCAACTCCGTTCTCATCAAATTCCGCATCTTCCTTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCCAGCAAAAGCCTGGTTCCGCCGATTTATCg                                       >  1:495249/1‑251 (MQ=255)
                                                                                                                                                                    gTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGGTCCTCTTCGAACAAATCGTTATCAATGGGGTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGGTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTa                                                                                                                      >  1:265784/1‑162 (MQ=255)
                                                                                                                                                                        cgTCCAACACACGCTGTGTGGATCCCTGCGATAACAGCTACTCTTCGACCACATCCTTATCACTGGGTTCCTGTTCTACTATGGACTTTTCGTCAAACACCTCAACTCCGTTCTCATAAACTTCCGCATCTTCATTTAAACCATCAACTTTCAGCGAATTCTGACTCATAAGCTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGATCCGCCGATTTATCGCTGCGGCAAATAAc                         >  1:336534/1‑251 (MQ=255)
                                                                                                                                                                            cAACACACGCTGTGTGGCTCCCTGCGATAATAGTTCCTCTTCGACCAAACGGTTATCCCTGGGTCACTGTTCTACTACGGCCTTTTCGTCAAAACCCt                                                                                                                                                                              >  1:461967/1‑98 (MQ=255)
                                                                                                                                                                                acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATNATGNNNNNTcnnnnnnnnnn                                                                                                                  >  2:306365/1‑144 (MQ=255)
                                                                                                                                                                                       gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCg                                                                                          <  1:438654/171‑1 (MQ=255)
                                                                                                                                                                                       gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAcn           <  1:306365/250‑2 (MQ=255)
                                                                                                                                                                                               cccTGCGATAACAGTTCCTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGTCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCTGGTTTACGGATTTcccn  <  1:150645/251‑2 (MQ=255)
                                                                                                                                                                                                               ccTCTTCGACCAAATCGTTATCACTGGGTTCCTGTTCTACTCAGGCCTTTTCGTCAAAAACCTCAACTCCGTTATCCTCAAATTCCGCATCTTCa                                                                                                                                              >  1:253046/1‑95 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  
CACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCC  >  NC_000913/2867215‑2867656

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