Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,123,110 G→A 100% 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,1100GA99.0% 83.4 / ‑6.8 35P77L (CCA→CTA) ftsNcell division protein FtsN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/17);  minor base C (1/1);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGA  >  NC_000913/4122970‑4123241
                                                                                                                                            |                                                                                                                                   
gTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtagt                                                                                                                                   <  2:846330/143‑1 (MQ=255)
 ttCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGGTCCGGCGTTGTCACTTCACCACCGGCTGAAGGTTCTGTGGGCGCGCGCGCTCCGGGCTGGCGACGTTCCAGCTCTTGAATGTAGCGCCAGCGTTCTTCTGGTTTtagtg                                                                                                                                  >  2:702986/1‑143 (MQ=255)
 ttCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGCCTTTCCAGCTCTTTAATGTAGCGCCAGCGCTCCTCTGGTTTtagtg                                                                                                                                  >  2:238262/1‑143 (MQ=255)
 ttCAAGAAGCTGACGTTGCTCTGGTGTCAGTTGCTCCCGCGTTTTCACTTCCACACCGGCCGAAGGTTCTGTGGGCCCACGCACTCCCGGCTGGCCAGTTTTTATTCCTTTTATGTAGCGCCAGCGTTCTTCTGGTTTtagtg                                                                                                                                  <  2:316377/143‑1 (MQ=255)
    aagaagCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTGCCAGGTCTTTAATGTAGCGCCCGCGTGCATCGGGTGTTAGGGGTa                                                                                                                               >  2:842083/1‑143 (MQ=255)
        agCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTcc                                                                                                                           <  1:395167/143‑1 (MQ=255)
         gCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCg                                                                                                                          >  1:524140/1‑143 (MQ=255)
            gACGTTGTTCTGGTGACAGTTGCTCAGGCACTGGCACGCCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTCTCCAGCTCTTTAACGTCGCGCCAGCGTCCTTCTGGTTTCCGTGGTAGTCCGttt                                                                                                                       >  2:275551/1‑143 (MQ=255)
               gttgttGTGGTGTCCGTTGCTCCGGGGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACACACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGGCCGTTTCCg                                                                                                                    <  2:857131/143‑1 (MQ=255)
                    tCTGGTGTCAGTTGCTCCGGCGTTTTCGCTGCGCCTGCGGCAGAGGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCAc                                                                                                               <  2:331683/143‑1 (MQ=255)
                               ttGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCttt                                                                                                    >  1:808965/1‑143 (MQ=255)
                                  cTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCAGCTGGCGACTATCCAGCTCTTTAGTGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCGGGTCACTTACTTGCTTTGCa                                                                                                 >  1:61051/1‑143 (MQ=255)
                                    ccGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGc                                                                                               >  1:487688/1‑143 (MQ=255)
                                        cGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCCCCTTCTGGCTTCGCACCGcct                                                                                           >  1:171279/1‑143 (MQ=255)
                                          ttttCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCg                                                                                         <  1:458389/143‑1 (MQ=255)
                                            ttcacttcacCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTTCAGCGCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGa                                                                                       >  1:697817/1‑143 (MQ=255)
                                               cctcggccaccGGCAGAAGGGTCTGTGGGCGCACGCGCTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTGGCGCCAGCGTTTTTCTGGTTTGAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGActc                                                                                    <  2:85296/137‑1 (MQ=255)
                                                 ttcacCACCGGCAGAAGGTTCTGTGGGCGGACGCACTCCCGACTGGGGACTTTCCAGCTCTTCAATGTAGCGCCAGCGTTCGTCTGGTTTTAGTGGTAGTCCGTTTCCGGTGACAATTCGGATTTGCAGAGTCTCGGATtctt                                                                                  >  2:218592/1‑143 (MQ=255)
                                                       accGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGCCCTCttctttct                                                                            >  1:583900/1‑143 (MQ=255)
                                                         cGGCAGGAGGTTCGCTGGGCGCGCGCTCGCCCGGCTGGCGACTGTCCAGCGCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTtg                                                                          <  2:752572/143‑1 (MQ=255)
                                                            ctgtacgTTCTGTTGGCGCACGCACTGCCGGTTGTGGACTTTTCAGCTCTTTTTTGTGGCGCTAGCGTTCTTCTGTTTTTCGTCGTACTCCGTTTCCTGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGa                                                                       <  2:125367/137‑1 (MQ=255)
                                                            cAGAAGGTGCTGTGGGCGCACGCACCCCCGGCTGGCGACTTTCCAGCTCTTTCATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGa                                                                       <  1:238262/143‑1 (MQ=255)
                                                                       gtgGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGAGGTGGGCCGGTGGTACTGGTCGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTGCTTTCTTGTGCTGCGTAATGac                                                            >  1:847858/1‑142 (MQ=255)
                                                                       gggggCGCACGGGCGCCCGGCTGGCGGCTTGCCGGCGCTTTGGTGTGGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCGGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGaa                                                            <  2:266446/141‑1 (MQ=255)
                                                                          ggCGCACGCGGTGCCGTCTGGCGACTTTCCGGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTa                                                         <  2:527358/143‑1 (MQ=255)
                                                                             gcacgcacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAg                                                      <  1:348259/143‑1 (MQ=255)
                                                                                cgcacTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTCCAGacc                                                   >  2:358761/1‑143 (MQ=255)
                                                                                         ggCTGGCGACTTTCCAGCTCTTTAATGTGGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTTATGCGTAATGAAGTACAGACCACCGATaaa                                          <  2:601632/143‑1 (MQ=255)
                                                                                             ggCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGGTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGCAGTACAGACCACCGGTAAAGGTc                                      >  1:152059/1‑143 (MQ=255)
                                                                                                    ttCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCAATAAAGGTCAaaaaaa                               >  1:521412/1‑143 (MQ=255)
                                                                                                      ccAGCTCTTTAATGTAGCGCCAGCGTTCTTCAGGTATTAGCGGCACTCCCTTTCCGGTCACTTTCTGGCTTTGCAGCGTCGCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAGAGGGCACAAGAAcg                             >  1:896904/1‑143 (MQ=255)
                                                                                                        aGCTCTTTAATGTCGCGCCAGCGCTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAAcggc                           <  1:388209/143‑1 (MQ=255)
                                                                                                               tAATGTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGGGTCTCGGACTCTTCTTTCTTGTGATGCGTGATGAAGTCGGGACCCCCGATAAAGGTCACCAGCAcggcgctagcc                    >  2:600828/1‑137 (MQ=255)
                                                                                                                  ggTAGCGCCTGCGTTCTTCTTGTTTTAGTGGTAGTGCGTTTCCGTTCACTTTCTGGCTTTGCAGGGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATa                 <  2:75788/142‑1 (MQ=255)
                                                                                                                   gTAGCGCCAGCGTTCTTCTGGTTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAg                <  1:759610/143‑1 (MQ=255)
                                                                                                                      gcgcCAGCGTTCTTCTGGTCTTAGTAGTAGTTCGATGCGGGTCACTTTCTGCCTTTGCAGCTTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGACACAAGAACGGCAGCAGCAATAGCGa             >  2:492585/1‑143 (MQ=255)
                                                                                                                      gcgcCAGCGTTCTTCTGGGTTTAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACGCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGa             <  1:304498/143‑1 (MQ=255)
                                                                                                                                 tttctGGTTTGAGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGAGAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGga  <  2:164856/142‑1 (MQ=255)
                                                                                                                                 cttgtGGTTTTAGTGGTAGTCTGTGTCCGGTCGCTTTCTGGTGTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGga  <  2:362114/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                   
GTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGTTTCCGGTCACTTTCTGGCTTTGCAGCGTCTCGGACTCTTCTTTCTTGTGATGCGTAATGAAGTACAGACCACCGATAAAGGTCACAAGAACGGCGGCAGCAATAGCGACCATAGCGGGA  >  NC_000913/4122970‑4123241

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