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_000913	3,683,630	0	C	T	100.0%	169.1 / NA	54	*663* (TAG→TAA)	yhjK	cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/27);  total (27/27)

GCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGC  >  NC_000913/3683497‑3683761                                                                                                                                      |                                                                                                                                    gCGGCCATAAAATAACACCCAGACAACAGCACAGTAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTTACTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACtt                                                                                                                                <  2:144408/139‑1 (MQ=255)        tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAGCAAAAAGGTTAAAAAGTGTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTcc                                                                                                                         <  2:252228/139‑1 (MQ=255)               acacCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  >  2:408681/1‑139 (MQ=255)                       acaacaTCACAGAAATGTACCTGGATCATAATTGTAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                          <  1:168872/139‑1 (MQ=255)                           caTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      >  1:1029271/1‑139 (MQ=255)                                      tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAggg                                                                                           >  1:637587/1‑139 (MQ=255)                                        tACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  >  2:1052644/1‑114 (MQ=255)                                        tACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  <  1:1052644/114‑1 (MQ=255)                                             ggATCATAATTGAAACACAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                           >  1:808607/1‑100 (MQ=255)                                             ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                           <  2:808607/100‑1 (MQ=255)                                              gATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCg                                                                                   <  2:1029271/139‑1 (MQ=255)                                               aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                           <  1:1254205/98‑1 (MQ=255)                                               aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                           >  2:1254205/1‑98 (MQ=255)                                               aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                          <  1:1170487/115‑1 (MQ=255)                                               aTCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                          >  2:1170487/1‑115 (MQ=255)                                                        gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                         <  2:1064645/139‑1 (MQ=255)                                                        gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCCCACAAACCAACACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                         <  1:466650/139‑1 (MQ=255)                                                         aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      <  1:1032333/109‑1 (MQ=255)                                                         aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      >  2:1032333/1‑109 (MQ=255)                                                         aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGa                                                                        >  1:313353/1‑139 (MQ=255)                                                             aaaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGcc                                                                    <  1:536003/139‑1 (MQ=255)                                                              aaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGccc                                                                   <  1:408681/139‑1 (MQ=255)                                                                              tGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTTTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                <  2:539548/94‑1 (MQ=255)                                                                              tGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTTTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                >  1:539548/1‑94 (MQ=255)                                                                                          cacaAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      <  2:169952/76‑1 (MQ=255)                                                                                          cacaAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      >  1:169952/1‑76 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  >  1:396071/1‑59 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  <  2:396071/59‑1 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      <  2:248837/71‑1 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                      >  1:248837/1‑71 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAg                                                                                         <  2:548845/84‑1 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAg                                                                                         >  1:548845/1‑84 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag                                  >  2:362783/1‑139 (MQ=255)                                                                                               aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag                                  >  1:265783/1‑139 (MQ=255)                                                                                                   gCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  <  2:620621/55‑1 (MQ=255)                                                                                                   gCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                  >  1:620621/1‑55 (MQ=255)                                                                                                     aCTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGt                            >  1:846035/1‑139 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                          <  2:566940/41‑1 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttcttc                                                                                                                          >  1:566940/1‑41 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                    >  2:992704/1‑63 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                    <  1:992704/63‑1 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                >  1:705240/1‑67 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                <  2:705240/67‑1 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCaa                                                            <  1:676726/103‑1 (MQ=255)                                                                                                         tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCaa                                                            >  2:676726/1‑103 (MQ=255)                                                                                                                   tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAATATTTCAATAGGGAGTGGGCGATCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGtctc              <  2:286464/139‑1 (MQ=255)                                                                                                                   tGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCACGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGtctc              >  1:249746/1‑139 (MQ=255)                                                                                                                              tGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCa                                                                               >  2:469840/1‑63 (MQ=255)                                                                                                                              tGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCa                                                                               <  1:469840/63‑1 (MQ=255)                                                                                                                               ggggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTcc             <  2:1030130/128‑1 (MQ=255)                                                                                                                               ggggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTcc             >  1:1030130/1‑128 (MQ=255)                                                                                                                               ggggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGc  >  1:702055/1‑139 (MQ=255)                                                                                                                                 ggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                     >  2:279659/1‑102 (MQ=255)                                                                                                                                 ggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                     <  1:279659/102‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    GCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGC  >  NC_000913/3683497‑3683761

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