Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,123,110 G→A P77L (CCA→CTA)  ftsN ← cell division protein FtsN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,123,1100GA95.0% 95.2 / NA 40P77L (CCA→CTA) ftsNcell division protein FtsN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/23);  minor base C (2/0);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCA  >  NC_000913/4122970‑4123250
                                                                                                                                            |                                                                                                                                            
gTTCAAGAAGCTCACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGACGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtagt                                                                                                                                            <  1:570371/143‑1 (MQ=255)
  tCAAGGCGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtagtgg                                                                                                                                          <  2:799182/143‑1 (MQ=255)
    aagaagCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTa                                                                                                                                        <  1:218582/143‑1 (MQ=255)
                ttgttCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTCCTGTAGCGCCCGCGTTCTTCTGGTTTTAGGGATAGTCCGTTGccag                                                                                                                            >  1:283035/1‑141 (MQ=255)
                      tGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTCTCCAGCTCTTTAATGTAGCGCCAGCCTTCTTCTGGTCTTCGTGGTAATCCATTTCCgctaacgt                                                                                                                      >  1:256058/1‑136 (MQ=255)
                        gtgtCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTc                                                                                                                    >  1:595298/1‑143 (MQ=255)
                           tCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTgg                                                                                                                 <  1:694870/143‑1 (MQ=255)
                             aGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTCATGTCGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCAGTCACTTTATTGCt                                                                                                               >  1:325018/1‑143 (MQ=255)
                              gTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGCCTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTCGTCCGTTTCCGGTCACTTTCTGGCtt                                                                                                              >  1:272179/1‑143 (MQ=255)
                                 gCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGc                                                                                                           <  1:187911/143‑1 (MQ=255)
                                   tCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAg                                                                                                         >  1:655820/1‑143 (MQ=255)
                                    ccgccgCTTTCGCTTCACCTCCGGCAGAAGGTTCTTTGGGCTCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTCATTTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTAACTGTGTGGCTTTGCAGc                                                                                                        <  2:59673/143‑1 (MQ=255)
                                     cGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCCGCGTGCTTCTGGTGTTCGTGGTGGTCGGTCTCCGGTCGCGTTCTGGCTTTgtcgcg                                                                                                       >  2:705110/1‑138 (MQ=255)
                                     cGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTACTGGTCGTCCGTTTCCGGTCACTTTCTGCCTTTGccccg                                                                                                       >  2:337136/1‑139 (MQ=255)
                                         gTTTTCACTTCACCACCGGCAGAAGGTGCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTGTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGtctc                                                                                                   <  1:337136/143‑1 (MQ=255)
                                         gTTTTCACTTCACCACCGGCAGAAGGCTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGtctc                                                                                                   <  1:465305/143‑1 (MQ=255)
                                          ttttCACTTCACCCCCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTGGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCg                                                                                                  <  1:98540/143‑1 (MQ=255)
                                           tttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTcc                                                                                                                              <  2:291670/114‑1 (MQ=255)
                                           tttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTcc                                                                                                                              >  1:291670/1‑114 (MQ=255)
                                           tttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCgg                                                                                                 >  1:742/1‑143 (MQ=255)
                                                         cGGCAGAAGGTTCTGTGGGCGCTCGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTtg                                                                                   <  1:309572/143‑1 (MQ=255)
                                                                    tCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAAt                                                                        <  1:21592/143‑1 (MQ=255)
                                                                     cTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTATGATGCGTAATg                                                                       >  1:748243/1‑143 (MQ=255)
                                                                          ggCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTa                                                                  <  1:604913/143‑1 (MQ=255)
                                                                          ggCGCACGCACTCCCGGCTGGCGACTTGCCAGGTCGTTAACGTAGCGCCAGCGTTCTTCTGGCCGTAGTGGTAGTCCGTTTCGCGTCACTTTTTGGCTTTGGAGCGGCTCGGCCTCTTCTCTCCTGTGATGCGTCATGACGtg                                                                  >  2:115513/1‑142 (MQ=255)
                                                                           gcgcACGCGCTCCCGGCTGGCGACTTGCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTAc                                                                 <  2:628556/143‑1 (MQ=255)
                                                                                     tCCCGGCTGGCGACTTTCGAGCCCTTTGATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGa                                                       <  2:237496/143‑1 (MQ=255)
                                                                                     tCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCGTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGa                                                       >  2:169966/1‑143 (MQ=255)
                                                                                        cGGCTGGCGATTTTCCAGCTCTTTAATGTAGCGCCATCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATaa                                                    <  2:645031/143‑1 (MQ=255)
                                                                                             ggCGCCCTTGCAGCTCGTTAATGTAGCGCCCGCGTTCTTCTGGTTTTAGTGGTACTCCGTTTCCTGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTTATGCTTAATGAAGTACAGACCACCGATAAAGGTc                                               <  1:455360/143‑1 (MQ=255)
                                                                                              gCGTCTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGCTAGTCCGTGCCTGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTca                                              <  2:693482/143‑1 (MQ=255)
                                                                                                gACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATTCGCAACTAAACACAGACCACCGATAAAGGTcaca                                            <  2:655820/143‑1 (MQ=255)
                                                                                                   tttCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGTTGCGTGATGAAGTACAGACCACCGATCACGGGCACAAGa                                         >  2:53515/1‑143 (MQ=255)
                                                                                                         gCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATACAGGTCACAAGAAcggcg                                   >  1:86585/1‑143 (MQ=255)
                                                                                                          ctctTTAATGTGGCGCCAGCGTTCTTCTGGTTTTAGGGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAAcggcgg                                  <  2:371178/143‑1 (MQ=255)
                                                                                                          ctctTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAGGTACAGCCCACCGATAAAGGTCACAAGAAcggcgg                                  >  1:62773/1‑143 (MQ=255)
                                                                                                             tttAATGTATCGCCAGCGTTCTTCTTGTTTTAGTGGTTGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGgcag                               <  1:123617/143‑1 (MQ=255)
                                                                                                                          cAGCGTTCTTCTGGTTTTAGTGGTAGTCCGGTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCAt                  <  2:595298/143‑1 (MQ=255)
                                                                                                                               ttcttcTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCgg             <  1:169966/143‑1 (MQ=255)
                                                                                                                                 cttccGGTTGTAGTGGTAGTCCGTGTCCGGGCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGga           <  1:311096/143‑1 (MQ=255)
                                                                                                                                    ctGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTCCAGACCACCGATAAAGGTCACACGAACGGCGGCAGCAATAGCGCCCATAGCTGgagca        >  2:557950/1‑141 (MQ=255)
                                                                                                                                     tgtttttAGTGGTAGTCCGTTTCGCGTCACTTTCTGGCTTTGCAGCGTCTCTGACTCTTGTTTCTTGTGCTTCTTAATGAACTACAGACCACCGCTAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGaaa       <  2:283035/140‑1 (MQ=255)
                                                                                                                                          tgggtggtAGTGCGTTGGCGGTGCTTTTGTGGCTTTGCAGCGTCTGGGACTCTGCCTGCTTGTGATGCGTAAGGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCgca  <  2:383959/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGAGAAACCGCA  >  NC_000913/4122970‑4123250

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