Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,693,326 C→G R637R (CGG→CGC bcsA ← cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,693,3260CG93.3% 46.2 / ‑5.6 15R637R (CGG→CGCbcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (6/8);  minor base T (1/0);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCA  >  NC_000913/3693217‑3693467
                                                                                                             |                                                                                                                                             
acacCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGt                                                                                                        >  3:33661/1‑149 (MQ=255)
        cTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCtt                                                                                                 <  2:4258/148‑1 (MQ=255)
           acgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTg                                                                                               <  4:33661/147‑1 (MQ=255)
           acgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGc                                                                                             >  2:158388/1‑149 (MQ=255)
              acCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGt                                                                                          <  1:158388/149‑1 (MQ=255)
              acCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGt                                                                                          <  2:769/149‑1 (MQ=255)
                                        gCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTcc                                                                                                      <  7:20281/111‑1 (MQ=255)
                                        gCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTcc                                                                                                      >  8:20281/1‑111 (MQ=255)
                                                                tACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCttt                                                                        <  3:42710/117‑1 (MQ=255)
                                                                tACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCttt                                                                        >  4:42710/1‑117 (MQ=255)
                                                                           ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGTCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCAACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAc                             >  5:32938/1‑149 (MQ=255)
                                                                             accaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTg                                                                             <  5:25864/99‑1 (MQ=255)
                                                                             accaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTg                                                                             >  6:25864/1‑99 (MQ=255)
                                                                                              aaGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCg           >  5:169497/1‑148 (MQ=255)
                                                                                                      tGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCa  <  6:169497/149‑1 (MQ=255)
                                                                                                             |                                                                                                                                             
ACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCA  >  NC_000913/3693217‑3693467

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