Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	875,327	G→A	*783* (TAG→TAA)	yliE →	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	875,327	0	G	A	100.0%	131.7 / NA	41	*783* (TAG→TAA)	yliE	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

TGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCC  >  NC_000913/875192‑875457                                                                                                                                        |                                                                                                                                   tGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAAttt                                                                                                                                 >  1:1055425/1‑139 (MQ=255)           gTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTcc                                                                                                                       >  2:1150953/1‑139 (MQ=255)                      aGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATaa                                                                                                            >  2:942978/1‑139 (MQ=255)                           ttGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTTCAGAATCAATAAGTTCg                                                                                                       <  1:942978/139‑1 (MQ=255)                           ttGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                       <  2:1055425/139‑1 (MQ=255)                                  cgAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAc                                                                                                >  1:1033186/1‑139 (MQ=255)                                     aTTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                             >  1:1144151/1‑139 (MQ=255)                                     aTTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt                                                                                             <  1:1150953/139‑1 (MQ=255)                                         aTCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagtcagt                                                                                         <  2:1033186/139‑1 (MQ=255)                                                         ttCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCaa                                                                                                               >  2:164111/1‑101 (MQ=255)                                                         ttCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCaa                                                                                                               <  1:164111/101‑1 (MQ=255)                                                                    gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATATCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagca                                                              >  2:1247725/1‑139 (MQ=255)                                                                        aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                             <  2:1160760/72‑1 (MQ=255)                                                                        aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGt                                                                                                                             >  1:1160760/1‑72 (MQ=255)                                                                         tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGc                                                         <  1:290870/139‑1 (MQ=255)                                                                           tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATATCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCt                                                       <  1:1247725/139‑1 (MQ=255)                                                                            gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGcc                                                        >  1:253899/1‑137 (MQ=255)                                                                            gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGcc                                                        <  2:253899/137‑1 (MQ=255)                                                                            gTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAG‑GGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                     >  1:1180879/1‑139 (MQ=255)                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                    <  1:691084/107‑1 (MQ=255)                                                                              aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTctgct                                                                                    >  2:691084/1‑107 (MQ=255)                                                                                   aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                               >  2:1204217/1‑139 (MQ=255)                                                                                   aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                               >  1:1353290/1‑139 (MQ=255)                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                   <  2:330004/131‑1 (MQ=255)                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                   >  1:837382/1‑131 (MQ=255)                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                   >  1:330004/1‑131 (MQ=255)                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                   <  2:837382/131‑1 (MQ=255)                                                                                       ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                           <  2:1353290/139‑1 (MQ=255)                                                                                           gCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCa                                       <  2:802975/139‑1 (MQ=255)                                                                                               ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATg                                   >  2:946750/1‑139 (MQ=255)                                                                                                    cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTaaa                              >  1:1248259/1‑139 (MQ=255)                                                                                                     gCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATc                                                                          <  2:301761/94‑1 (MQ=255)                                                                                                     gCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATc                                                                          >  1:301761/1‑94 (MQ=255)                                                                                                           gATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACgg                       >  1:902333/1‑139 (MQ=255)                                                                                                            aTAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                       >  2:1382603/1‑58 (MQ=255)                                                                                                            aTAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg                                                                                                       <  1:1382603/58‑1 (MQ=255)                                                                                                                aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                          >  1:886774/1‑131 (MQ=255)                                                                                                                aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                          <  2:886774/131‑1 (MQ=255)                                                                                                                aaGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta                  <  1:946750/139‑1 (MQ=255)                                                                                                                                aaaGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc  <  2:1144151/139‑1 (MQ=255)                                                                                                                                aaaGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACCTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc  <  2:902333/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   TGCGGAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCC  >  NC_000913/875192‑875457

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