Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,683,630 C→T *663* (TAG→TAA)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,6300CT100.0% 137.0 / NA 43*663* (TAG→TAA) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/20);  total (23/20)

GCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTC  >  NC_000913/3683500‑3683758
                                                                                                                                  |                                                                                                                                
gCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttc                                                                                                                          >  2:446405/1‑139 (MQ=255)
    tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTcc                                                                                                                      >  1:1423758/1‑139 (MQ=255)
                  gACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttc                                                                                                                          >  1:1308989/1‑121 (MQ=255)
                  gACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttc                                                                                                                          <  2:1308989/121‑1 (MQ=255)
                         tCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                 <  2:1423758/139‑1 (MQ=255)
                                  tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                       <  2:577506/108‑1 (MQ=255)
                                  tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                       >  1:577506/1‑108 (MQ=255)
                                  tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        >  2:307990/1‑107 (MQ=255)
                                  tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        <  1:307990/107‑1 (MQ=255)
                                    tACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAg                                                                                      >  2:112272/1‑139 (MQ=255)
                                                     aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGa                                                                     >  2:917304/1‑139 (MQ=255)
                                                       aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                             <  2:1578954/81‑1 (MQ=255)
                                                       aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                             >  1:1578954/1‑81 (MQ=255)
                                                       aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                             >  1:1336405/1‑113 (MQ=255)
                                                       aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                             <  2:1336405/113‑1 (MQ=255)
                                                          aaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGccc                                                                <  1:917304/139‑1 (MQ=255)
                                                                         ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                 <  1:186564/91‑1 (MQ=255)
                                                                         ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                 >  2:186564/1‑91 (MQ=255)
                                                                         ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                             >  1:1216182/1‑95 (MQ=255)
                                                                         ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                             <  2:1216182/95‑1 (MQ=255)
                                                                          tGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        >  2:510016/1‑67 (MQ=255)
                                                                          tGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        <  1:510016/67‑1 (MQ=255)
                                                                                        caAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                  <  1:112272/139‑1 (MQ=255)
                                                                                           aCCAGCTCTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTGCTTCCAGGTAACTCTCTTCGAAGATTTc                                                                                               <  2:1036795/75‑1 (MQ=255)
                                                                                           aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTCGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTc                                                                                               >  1:1036795/1‑75 (MQ=255)
                                                                                                 aCTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                      >  1:166617/1‑94 (MQ=255)
                                                                                                 aCTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                      <  2:166617/94‑1 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        >  2:220897/1‑40 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                        <  1:220897/40‑1 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                     >  2:1536433/1‑59 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                     >  2:341986/1‑59 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                     <  1:1536433/59‑1 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                     <  1:341986/59‑1 (MQ=255)
                                                                                                     tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtgtg                     >  2:1230058/1‑139 (MQ=255)
                                                                                                             tttGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                >  1:880790/1‑136 (MQ=255)
                                                                                                             tttGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                <  2:880790/136‑1 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                 <  2:1473647/53‑1 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                 >  1:1473647/1‑53 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTcgcg                          <  1:811980/124‑1 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTcgcg                          >  2:811980/1‑124 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGtctc           >  2:291952/1‑139 (MQ=255)
                                                                                                               tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGtctc           <  1:1230058/139‑1 (MQ=255)
                                                                                                                        tttGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTc  >  2:849981/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                
GCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTC  >  NC_000913/3683500‑3683758

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