Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,346,417	C→T	G109E (GGG→GAG)	pdeR ←	cyclic di‑GMP phosphodiesterase PdeR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,346,417	0	C	T	100.0%	73.1 / NA	22	G109E (GGG→GAG)	pdeR	cyclic di‑GMP phosphodiesterase PdeR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/2);  total (20/2)

TTGCCCAATGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCAT  >  NC_000913/1346281‑1346550                                                                                                                                          |                                                                                                                                      ttGCCCAATGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCgt                                                                                                                                    <  2:729182/140‑1 (MQ=255)        aTGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATaa                                                                                                                             >  1:289510/1‑140 (MQ=255)        aTGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATaa                                                                                                                             >  1:226653/1‑140 (MQ=255)             gtcgtGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCAc                                                                                                                        <  2:264191/140‑1 (MQ=255)                 tgtTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGc                                                                                                                    >  2:347673/1‑140 (MQ=255)                    tttttCACCCCTGTGTAATCTTCACATAACCGATTGAGGCGTGGGATATCCCCCCGGCTGTCGAGAATCACTTTTGCGGAATGTGCTTCAGAAA‑CTCCCTGCTCAGAAAAGGAAAGCTCGTGTGAAAAATCACGAGAAAc                                                                                                                 <  2:138189/138‑1 (MQ=255)                           aCCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGgt                                                                                                          >  2:939769/1‑140 (MQ=255)                              ttgtgTAATCTTCACAAACCCGATTGAGGCGTTGGATTTCCCCCCGGCTGTGAAAAACCATTTTTGGGGAATGTGCTTCAAAAC‑CTCCCTCCCCAGAAAAGAAAAGCTGGTGTAAAAAACCACGAGCAACCGCCGgtgtg                                                                                                       <  2:307372/139‑1 (MQ=255)                              cTGTGTAATCTTCACATACCAGATTGAAGCGTTGGATATCCCCCCGGCTGTCGAAAATCTCTTTTGGGGAATGTGCTTCAAAAA‑CTCCCTGCTCAGAAAAGAAAAGCTGGTGTAAAAAATCACGAGCAACCGCCGgtgtg                                                                                                       <  2:289510/140‑1 (MQ=255)                              cTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATTTCCCCCCGGCTGTCGAGAATCATTATTGCGGAATGTGCTTCAAAAA‑CTCCCTGCTCAGAAAAGGAAAGCTCGTGTAATAAATCACGAGCAACCGCCGgtgtg                                                                                                       <  2:744766/140‑1 (MQ=255)                               tgtgTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGgtgtgt                                                                                                      >  1:441562/1‑140 (MQ=255)                                gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTc                                                                                                     >  1:1158079/1‑140 (MQ=255)                                   aaaTCTCAACAAAGCCAATTGAGGGTTTGAATCTCCCCCCGGCTAGCGAGATTCCTTTTTGGGGAATCTGTTTCAAAAG‑CCCCTTCCTAAGAAAGGAAAAGCTCGTGTGAAAAACCCCGAGCAACCGCCGGTGTGTcatg                                                                                                  <  2:619660/139‑1 (MQ=255)                                       cTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTcatgccat                                                                                              >  1:956182/1‑140 (MQ=255)                                           acaTAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAg                                                                                          >  1:1173602/1‑140 (MQ=255)                                                    aTTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGc                                                                                 >  1:1081853/1‑140 (MQ=255)                                                    aTTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGc                                                                                 >  1:969165/1‑140 (MQ=255)                                                    aTTGAAGCGTTGGATATCCCCCCGGCTGTCGAAAATCACTTTTGCGGAATGTGCTTCAGAAA‑CTCCCTGCTCAGCAAAGGAAAGCTCGTGTGAAAAACCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGc                                                                                 <  2:1081853/140‑1 (MQ=255)                                                       gAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAg                                                                              >  2:1125208/1‑140 (MQ=255)                                                         aGCTTTGAATATCCCCCCGGCTGCGGAGAATCATTTTTGCGGAATGTGCTTCAAAAA‑CTCCCGGCTCAGAAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGcc                                                                            <  2:1173602/140‑1 (MQ=255)                                                                                 cGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGt                                                    >  2:1156134/1‑140 (MQ=255)                                                                                            tttGGGGATTGTG‑TTCAAAAACCTCCCTGCACAGAAAAGAAAAGCTGGTGTGAAAAACCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATg                                         <  2:956182/139‑1 (MQ=255)                                                                                               gCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGcac                                      >  2:642965/1‑140 (MQ=255)                                                                                                  gAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGa                                   >  2:905991/1‑140 (MQ=255)                                                                                                     tgtgCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAg                                >  2:765811/1‑140 (MQ=255)                                                                                                               aaa‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCa                      >  1:865637/1‑140 (MQ=255)                                                                                                                 a‑CTACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATc                    >  1:1173584/1‑140 (MQ=255)                                                                                                                    tACCTGCTCAGCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCACTTTCATCGa                  >  1:654929/1‑140 (MQ=255)                                                                                                                              gCAAAGGAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGCGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACt        >  2:658259/1‑140 (MQ=255)                                                                                                                                   gaaaaGTCCGTGTAAAAAATCCCGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatca   <  1:939769/138‑1 (MQ=255)                                                                                                                                    gAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat  <  1:905991/140‑1 (MQ=255)                                                                                                                                    gAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat  >  1:80340/1‑140 (MQ=255)                                                                                                                                    gAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat  <  1:642965/140‑1 (MQ=255)                                                                                                                                    gAAAGCTCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat  <  1:1125208/140‑1 (MQ=255)                                                                                                                                    gAAAGCTCGTGTAATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat  <  1:892793/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TTGCCCAATGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAA‑CTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCAT  >  NC_000913/1346281‑1346550

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