Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,894	G→A	Q242* (CAA→TAA)	pdeB ←	c‑di‑GMP phosphodiesterase PdeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,894	0	G	A	100.0%	104.7 / NA	31	Q242* (CAA→TAA)	pdeB	c‑di‑GMP phosphodiesterase PdeB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

CAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCATCAATTTGCCATGAGCTATAGGGAATGACATCAATAA  >  NC_000913/477670‑478136                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   caaccaacTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACAttt                                                                                                                                                                                                                          >  2:158496/1‑251 (MQ=255)                   ctgCGGCCAGACCGCCAGTGCCTCAGCCCCGACAATTTTGCCATTGGCTAAGGAGACAACCGGCTGATAGCGCACGCCAATATCACGATTTTCGATAGCCTCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAACTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACAGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGACCAtgtgt                                                                                                                                                                                                        >  2:177080/1‑250 (MQ=255)                    ggCGGCCAGCGAGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACAGCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                                                                                                                                                                                                       <  2:153552/250‑1 (MQ=255)                      cGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAACCCATGtgatg                                                                                                                                                                                                      >  2:331975/1‑250 (MQ=255)                       ggCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTgatgat                                                                                                                                                                                                    <  1:101259/251‑1 (MQ=255)                           aGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCATCCCGAGCGGTAACCAAATAAAGACTTACCGATGCCa                                                                                                                                                                                                                                           <  1:394145/208‑1 (MQ=255)                           aGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCa                                                                                                                                                                                                                                           >  2:394145/1‑208 (MQ=255)                                gcCAGTGCTTCAGCACCGACAATTATTCCATTGGCTAAGGAGACAATCGGCTGATAGTGCCCGCAAATATCCCGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACTGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCTGATCGGTAACCAAATAAAGCCTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATTATCGAAATAtt                                                                                                                                                                                           <  2:17967/251‑1 (MQ=255)                                          cAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCTCGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATTCGCAGCACAAACATCGCTGCCAGCAGTCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCgg                                                                                                                                                                                  <  1:268779/250‑1 (MQ=255)                                           aGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTa                                                                                                                                                                                <  1:421844/251‑1 (MQ=255)                                                           ccATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGt                                                                                                                                                                 >  2:80834/1‑250 (MQ=255)                                                                tgCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATg                                                                                                                                                            <  2:113851/249‑1 (MQ=255)                                                                       gagaCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCatt                                                                                                                                                    >  1:7783/1‑251 (MQ=255)                                                                        agaCAATCGGCTGATAGTGCAAGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGAGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCatt                                                                                                                                                    >  2:104872/1‑250 (MQ=255)                                                                          aCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCattatt                                                                                                                                                 <  1:391245/251‑1 (MQ=255)                                                                          aCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCattatt                                                                                                                                                 <  2:259035/251‑1 (MQ=255)                                                                                              gCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGgtg                                                                                                                             >  1:337383/1‑251 (MQ=255)                                                                                                       aCGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGt                                                                                                                     >  1:242143/1‑250 (MQ=255)                                                                                                        cGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGt                                                                                                                     >  1:182065/1‑249 (MQ=255)                                                                                                                        ccTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTccc                                                                                                   <  1:80834/251‑1 (MQ=255)                                                                                                                                    gatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCa                                                                                       <  1:169983/251‑1 (MQ=255)                                                                                                                                         ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCgct                                                                                  >  1:51573/1‑251 (MQ=255)                                                                                                                                                                   cagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCAt                                                        >  2:233775/1‑251 (MQ=255)                                                                                                                                                                    aggaccaacctCGCTGCCAGCAGGCCAATCACCAACCCGCGCGGTAAACAAATACAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGATGAAGCCAATGTGATGATCGAAATATTCATCTCCGCTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCAACGGGTGTTTTTTGTAGCATGGTAATAATTCCATGATAAATTTCATCGATGCTGGTTATCACAACGTTATGGGNAtt                                                       >  2:415353/11‑251 (MQ=255)                                                                                                                                                                             aTCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTg                                                                                                                      <  1:351466/179‑1 (MQ=255)                                                                                                                                                                             aTCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTg                                                                                                                      >  2:351466/1‑179 (MQ=255)                                                                                                                                                                                gCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGGTATGACAACGTTATGGGCATTGCCAATAATggc                                           <  1:116820/251‑1 (MQ=255)                                                                                                                                                                                     cagcagGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCAt                                      >  2:178071/1‑251 (MQ=255)                                                                                                                                                                                     cagcagGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCAt                                      >  1:91615/1‑251 (MQ=255)                                                                                                                                                                                                               gTAACCAAATAAAGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCATCAATTTGCCATGAGCTATAGGGAATg            >  1:257993/1‑251 (MQ=255)                                                                                                                                                                                                                         aaaGACTTACCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCATCAATTTGCCATGAGCTATAGGGAATGACATCAATaa  <  2:51573/251‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   CAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCTGCTGGTTATGACAACGTTATGGGCATTGCCAATAATGGCGGCATCAATTTGCCATGAGCTATAGGGAATGACATCAATAA  >  NC_000913/477670‑478136

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