Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 92,197 C→T A262V (GCC→GTC)  ftsI → peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091392,1970CT96.8% 97.1 / NA 31A262V (GCC→GTC) ftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACG  >  NC_000913/92075‑92316
                                                                                                                          |                                                                                                                       
ttcttcTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCtggtg                                                                                                               >  1:354668/1‑133 (MQ=255)
ttcttcTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCtggtg                                                                                                               >  1:439134/1‑133 (MQ=255)
      tACTGACAGGGAGGCTGCGCACAACCGGCCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTc                                                                                                         <  2:286948/133‑1 (MQ=255)
      tACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTc                                                                                                         >  2:22992/1‑133 (MQ=255)
      tACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTc                                                                                                         <  1:425263/133‑1 (MQ=255)
        cTGACAGCCAGGCAGCGCACAATCTGGCGCAGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGGGAACTGAACAACGAGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCaa                                                                                                       <  1:206313/133‑1 (MQ=255)
           aCAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAacac                                                                                                    >  1:554592/1‑133 (MQ=255)
           aCAGCCAGCCAGCGCCCACCCTTGCGCTGACTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCTAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAacac                                                                                                    <  1:301047/133‑1 (MQ=255)
                cAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTCGCGTCGTGCTGATGAATGTCAACACCGGTg                                                                                               >  1:298247/1‑133 (MQ=255)
                  ggCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGTGTCGTGCGGGTGGAGGTCAACTCCGGTGaa                                                                                             >  1:166341/1‑133 (MQ=255)
                       cgcACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCt                                                                                        >  1:521942/1‑133 (MQ=255)
                           caACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCg                                                                                    <  2:521942/133‑1 (MQ=255)
                           caACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCg                                                                                    <  1:515685/133‑1 (MQ=255)
                                 ggCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCt                                                                              >  2:62692/1‑133 (MQ=255)
                                        aGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGcc                                                                       >  1:138836/1‑133 (MQ=255)
                                        aGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGcc                                                                       <  2:439134/133‑1 (MQ=255)
                                              gatgaACGCACGCAGTCGCTGGTTGAGCGCGATCTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCAt                                                                 <  2:298247/133‑1 (MQ=255)
                                               atgaACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATa                                                                <  1:4162/133‑1 (MQ=255)
                                                tgaACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATAc                                                               <  2:138836/133‑1 (MQ=255)
                                                  aaCGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACaa                                                             <  1:22992/133‑1 (MQ=255)
                                                    cGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAAcc                                                           <  1:128093/133‑1 (MQ=255)
                                                                             aaCTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACATTCTGAGCGGGACGCCGaa                                  >  1:302307/1‑133 (MQ=255)
                                                                             aaCTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAAAAGCCCGTCCTACCACCATAACAATCTGAGCGCCACGCCGaa                                  >  2:269932/1‑133 (MQ=255)
                                                                                  aacaacGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGTTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAgg                             >  1:45185/1‑133 (MQ=255)
                                                                                  aacaacGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAgg                             >  2:424017/1‑133 (MQ=255)
                                                                                           gTGGCCTTTAGCAAGGCTGAATCTGGTACCGCCGTGCTGGTCGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTa                    <  2:201480/133‑1 (MQ=255)
                                                                                            tGGCCTTTAACAAGGCTGAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTaa                   >  2:575936/1‑133 (MQ=255)
                                                                                                         ggCTGAATCTGGTAGCGTCGTGCTGGTGGATGTTAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGGGCGGCACGACGAAAGAGGCGATGCGGAACCGTACCATCAcc      >  2:171798/1‑133 (MQ=255)
                                                                                                            tGAATCTGGTAGCGTCGTGCTGGTAGATATCAACACCGGTGCAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGAc   <  2:268653/133‑1 (MQ=255)
                                                                                                             gAATCTGGTATCGTCTTGCTGGTGAATGCCAACACCGGTGACGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACg  <  2:166341/133‑1 (MQ=255)
                                                                                                             gAATCTGGTAGCGTCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACg  >  1:599445/1‑133 (MQ=255)
                                                                                                                          |                                                                                                                       
TTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACG  >  NC_000913/92075‑92316

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑