Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,626,871 A→G C690R (TGC→CGC)  pdeF ← cyclic di‑GMP phosphodiesterase PdeF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,626,8710AG100.0% 66.9 / NA 20C690R (TGC→CGC) pdeFcyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/9);  total (11/9)

AGAGCACCGCCTCGCGGATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAG  >  NC_000913/2626793‑2626946
                                                                              |                                                                           
agagCACCGCCTCGCGGATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCgaa                                                                           <  2:418841/81‑1 (MQ=255)
   gCACCGCCTCGCGGATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCg                                                                        <  2:418842/81‑1 (MQ=255)
             gcgGATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTAt                                                              >  2:418872/1‑81 (MQ=255)
                   ctctTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGAt                                                         >  1:418873/1‑80 (MQ=255)
                      ttCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATa                                                        <  2:418843/78‑1 (MQ=255)
                         gTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATc                                                     >  1:418868/1‑78 (MQ=255)
                          ttttCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGa                                                 <  2:418844/81‑1 (MQ=255)
                          ttttCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGa                                                 >  1:418869/1‑81 (MQ=255)
                               aaCGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTAtt                                            >  2:418875/1‑81 (MQ=255)
                                  gTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGa                                         >  1:418874/1‑81 (MQ=255)
                                    aCTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGaca                                       >  2:418882/1‑81 (MQ=255)
                                       ctGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACaa                                    >  2:418877/1‑81 (MQ=255)
                                          ccACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAatat                                 >  1:418876/1‑81 (MQ=255)
                                             cTACCAGCATTTTCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGc                              <  1:418845/81‑1 (MQ=255)
                                                      ttttCTTCATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAAc                     <  1:418846/81‑1 (MQ=255)
                                                         tcttcATTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTg                   <  1:418847/80‑1 (MQ=255)
                                                              aTTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTc              >  1:418870/1‑80 (MQ=255)
                                                              aTTCGCGCCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCa             <  1:418848/81‑1 (MQ=255)
                                                                    gcCAGGTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAAtttt        >  1:418871/1‑80 (MQ=255)
                                                                         gTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAg  <  1:418850/81‑1 (MQ=255)
                                                                         gTGGCGAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAg  <  1:418849/81‑1 (MQ=255)
                                                                              |                                                                           
AGAGCACCGCCTCGCGGATCTCTTCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAG  >  NC_000913/2626793‑2626946

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