Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,162,216 G→A Q152* (CAG→TAG)  ftsP ← cell division protein FtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,162,2160GA95.5% 64.5 / NA 22Q152* (CAG→TAG) ftsPcell division protein FtsP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/11);  minor base T (0/1);  total (10/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.96e-01

GCAAGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGACGCGGATAGTCGGCCCCAGGTAACG  >  NC_000913/3161976‑3162453
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
gCAAGCCCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGGATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTg                                                                                                                                                                                                                                     <  2:242497/251‑1 (MQ=255)
              ccACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTAtt                                                                                                                                                                                                                       <  2:144420/251‑1 (MQ=255)
                                                             cGTAACACCAACAAAGGCGCTGCTTCCCGGTTCGTTTTATTCTGGCGTACCAAAGTGATCCAGCCGTGTATCCTGGATAATGACCGGAAAATCATCCACCCCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCTAATGGGCAGTTCgg                                                                                                                                                                        <  2:118550/251‑1 (MQ=255)
                                                                              cgccgcTTCCCGGTTCGGTGTATGCTGGCGGACCAAAGTTATCCAGCCGTTTATCCTGGATAATGTCCGTAAAATCATCCACACCATAATGGTGGGGGATAGGCAGCGCCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGTCCGTTATAGACCTACTGCGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGTCAGTACTGGTGCCCAGTCAGCGttt                                                                                                                                                       <  2:12045/251‑1 (MQ=255)
                                                                                                        ggCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGcc                                                                                                                             <  1:165196/251‑1 (MQ=255)
                                                                                                           gTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCAt                                                                                                                          >  1:286872/1‑251 (MQ=255)
                                                                                                                           gCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGGAcc                                                                                                          >  2:118880/1‑251 (MQ=255)
                                                                                                                               tttATCCTGGATAATGACCGGAAAATCATCCACACCATCATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGt                                                                                                       >  2:277761/1‑250 (MQ=255)
                                                                                                                                       ggATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCgg                                                                                                <  1:222181/249‑1 (MQ=255)
                                                                                                                                         aTAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTTCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGGAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGTACCGCGCATCAGCGTGCCTGGGTCCGGTAGCCCGGCCAc                                                                                            <  2:20252/251‑1 (MQ=255)
                                                                                                                                                  ggAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTga                                                                                   <  1:94970/251‑1 (MQ=255)
                                                                                                                                                    aaaaTCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTgaga                                                                                 <  2:12546/251‑1 (MQ=255)
                                                                                                                                                                             gggATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTa                                                        <  2:286872/251‑1 (MQ=255)
                                                                                                                                                                                aTAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAAt                                                     >  1:115921/1‑251 (MQ=255)
                                                                                                                                                                                     cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCt                                                >  1:155294/1‑251 (MQ=255)
                                                                                                                                                                                      aGCGACTTGCTGACTTCATCTGCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCtt                                               <  2:95877/251‑1 (MQ=255)
                                                                                                                                                                                                            ccaccaGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGAcg                         >  1:281095/1‑251 (MQ=255)
                                                                                                                                                                                                              accaGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGAcgcg                       >  1:29849/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGACGCGGATAGTCGGccc           >  1:116840/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGACGCGGATAGTCGGCCCCAGGTaa    <  1:287647/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTACCATCGTCGCCCTTCCAGACGCGGATAGCCGGCCCCCGGTAACg  >  2:271453/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGCAAATACGCTTAACATCGTCGTCCTTCCTGAcgctggtagttggctcacgggaacg  >  2:3821/1‑228 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GCAAGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATCGTCGCCCTTCCAGACGCGGATAGTCGGCCCCAGGTAACG  >  NC_000913/3161976‑3162453

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