Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,680,898 T→C R546G (AGA→GGA)  bcsC ← cellulose synthase complex outer membrane protein BcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,680,8980TC100.0% 62.0 / NA 18R546G (AGA→GGA) bcsCcellulose synthase complex outer membrane protein BcsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/13);  total (5/13)

CGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATT  >  NZ_CP009273/3680814‑3680983
                                                                                    |                                                                                     
cGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCt                                                                                  <  2:320061/90‑1 (MQ=255)
       cAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATg                                                                           <  1:93485/90‑1 (MQ=255)
            gATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCag                                                                      <  1:208050/90‑1 (MQ=255)
            gATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCag                                                                      <  1:375761/90‑1 (MQ=255)
            gATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCag                                                                      <  1:488195/90‑1 (MQ=255)
              ttAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCagag                                                                    <  2:475557/90‑1 (MQ=255)
                            aaTATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTa                                                      <  1:216687/90‑1 (MQ=255)
                                  gCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTa                                                >  2:10447/1‑90 (MQ=255)
                                     gTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAAc                                             >  2:40773/1‑90 (MQ=255)
                                      ttCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAAcc                                            <  1:355466/90‑1 (MQ=255)
                                      ttCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAAcc                                            <  2:461693/90‑1 (MQ=255)
                                               gcACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAAcc                                            <  2:332781/81‑1 (MQ=255)
                                               gcACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAAcc                                            >  1:332781/1‑81 (MQ=255)
                                                      gCAGGCTATTGATATGCGCCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCg                            <  2:261035/90‑1 (MQ=255)
                                                                     gcgcCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACt                         >  1:353519/1‑78 (MQ=255)
                                                                     gcgcCAGCGCCGCTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACt                         <  2:353519/78‑1 (MQ=255)
                                                                             gccgcTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTgcgc     <  2:302922/90‑1 (MQ=255)
                                                                                gcTCCGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCAtt  >  2:255624/1‑90 (MQ=255)
                                                                                    |                                                                                     
CGCTTTGCAGTCGATTAACCAGCTCCTGAATATTGCTGTTCCACTGCGCACGCGGCAGGCTATTGATATGCGCCAGCGCCGCTCTGTCCTGGTCATGACCAGAGAGGTACAGCCCGTAAGCGTAAACCTGCTCCGGGTCGTTCGACTTCTGCTGCGCCAGATTGCGCATT  >  NZ_CP009273/3680814‑3680983

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