Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,400,225 T→C E288G (GAG→GGG)  mreB ← dynamic cytoskeletal protein MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,2250TC97.2% 116.6 / ‑4.2 36E288G (GAG→GGG) mreBdynamic cytoskeletal protein MreB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGG  >  NC_000913/3400148‑3400301
                                                                             |                                                                            
gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCgg                                                                           <  2:644279/81‑1 (MQ=255)
gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCgg                                                                           <  2:644277/81‑1 (MQ=255)
gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCg                                                                            >  2:644333/1‑80 (MQ=255)
  aTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGgag                                                                         <  2:644281/81‑1 (MQ=255)
  aTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGgag                                                                         <  2:644280/81‑1 (MQ=255)
       ggTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTc                                                                    <  1:644282/81‑1 (MQ=255)
         tttCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCgg                                                                  <  2:644283/81‑1 (MQ=255)
            cttcCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAg                                                               >  2:644334/1‑81 (MQ=255)
               ccATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGcc                                                             <  1:644284/80‑1 (MQ=255)
                cATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTTGGAAGCCAg                                                           <  1:644285/81‑1 (MQ=255)
                  ttAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGtt                                                         <  2:644286/81‑1 (MQ=255)
                    aaCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTcc                                                       >  2:644332/1‑81 (MQ=255)
                        aaCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCcggcg                                                   >  2:644348/1‑81 (MQ=255)
                          cGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCggg                                                 <  2:644287/81‑1 (MQ=255)
                            gTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCa                                               >  1:644352/1‑81 (MQ=255)
                               aaGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACt                                             <  2:807660/80‑1 (MQ=255)
                                 ggTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGtt                                          >  1:644335/1‑81 (MQ=255)
                                  gTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc                                         >  1:644337/1‑81 (MQ=255)
                                  gTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc                                         <  1:644288/81‑1 (MQ=255)
                                   ttACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc                                         >  2:644356/1‑80 (MQ=255)
                                        cagcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGc                                   <  1:644289/81‑1 (MQ=255)
                                          gcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCaa                                 <  1:644290/81‑1 (MQ=255)
                                          gcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCaa                                 >  2:644359/1‑81 (MQ=255)
                                             gTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCa                              >  1:644338/1‑81 (MQ=255)
                                               gcgcCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAt                             >  1:644341/1‑80 (MQ=255)
                                                cgcCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAtt                            <  2:644291/80‑1 (MQ=255)
                                                    accaccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACcgc                       <  1:644292/81‑1 (MQ=255)
                                                      caccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCgcgc                     >  1:644355/1‑81 (MQ=255)
                                                        ccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTc                   <  2:644293/81‑1 (MQ=255)
                                                          ggTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTca                  >  2:644343/1‑80 (MQ=255)
                                                           gTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTcac                 >  1:644364/1‑80 (MQ=255)
                                                                   cATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGt        >  1:644336/1‑81 (MQ=255)
                                                                    aTGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTc       <  2:644294/81‑1 (MQ=255)
                                                                       ccGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAg     <  1:644295/80‑1 (MQ=255)
                                                                        cgcgCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCg   >  1:644340/1‑81 (MQ=255)
                                                                         gcgcCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCgg  >  1:644351/1‑81 (MQ=255)
                                                                             |                                                                            
GAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGG  >  NC_000913/3400148‑3400301

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