Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,684,604	G→A	Q339* (CAG→TAG)	yhjK ←	cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,684,604	0	G	A	100.0%	177.4 / NA	45	Q339* (CAG→TAG)	yhjK	cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/26);  total (19/26)

CGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGAATTTCTCGTTGCG  >  NC_000913/3684456‑3684750                                                                                                                                                     |                                                                                                                                                   cGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAc                                                                                                                                                  <  2:544047/151‑1 (MQ=255)     cGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCt                                                                                                                                                <  1:359694/149‑1 (MQ=255)         aaaTGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTg                                                                                                                                          <  2:404479/151‑1 (MQ=255)           aTGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGAt                                                                                                                                        <  2:166006/151‑1 (MQ=255)                   gAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGc                                                                                                                                 >  1:16346/1‑150 (MQ=255)                   gAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGc                                                                                                                                 >  1:347624/1‑150 (MQ=255)                               cTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTa                                                                                                                    >  1:498224/1‑151 (MQ=255)                                tGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                                                                                   <  1:563210/151‑1 (MQ=255)                                  tAAAGCTTTTCGGCGGTGAGAGCGCAGTGGAACATCGCCACGCCAATGCTCCAGTGCGGACGGAGTTGGATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTTCCACGGTTCCTGTgcac                                                                                                                 <  1:744092/151‑1 (MQ=255)                                  tAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTacac                                                                                                                 <  2:307935/151‑1 (MQ=255)                                  tAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTacac                                                                                                                 <  2:48476/151‑1 (MQ=255)                                         gTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTgg                                                                                                           >  1:552230/1‑150 (MQ=255)                                          ttCGGCGGTGAGAGCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCa                                                                                                         <  1:477501/151‑1 (MQ=255)                                                       tCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGt                                                                                            <  1:765100/151‑1 (MQ=255)                                                        cGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGt                                                                                            >  1:156449/1‑150 (MQ=255)                                                                aaCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATGGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAcc                                                                                     <  2:280296/149‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                                                                               >  1:563576/1‑150 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                                                                               >  2:438202/1‑150 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  1:449849/151‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  2:539344/151‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  1:61365/151‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  2:216416/151‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  1:640372/151‑1 (MQ=255)                                                                     cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATc                                                                              <  1:210375/151‑1 (MQ=255)                                                                            ccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGa                                                                       <  1:379468/151‑1 (MQ=255)                                                                            ccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGa                                                                       >  2:598646/1‑151 (MQ=255)                                                                               aTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACCCCGTTGGCAATGCCAGCAAAGTCATAACCGCTAATCTGCGCGAGGa                                                                    <  1:272365/151‑1 (MQ=255)                                                                               aTGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGa                                                                    >  1:361867/1‑151 (MQ=255)                                                                                  cTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCa                                                                 >  1:202177/1‑151 (MQ=255)                                                                                        tGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg                                                             <  2:66051/149‑1 (MQ=255)                                                                                        tGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg                                                             <  2:323530/149‑1 (MQ=255)                                                                                        tGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg                                                             <  2:194332/149‑1 (MQ=255)                                                                                        tGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg                                                             <  1:134675/149‑1 (MQ=255)                                                                                                gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGt                                                   >  1:320015/1‑151 (MQ=255)                                                                                                  gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTa                                                  >  2:464129/1‑150 (MQ=255)                                                                                                  gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTa                                                  >  2:621818/1‑150 (MQ=255)                                                                                                   ttGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTAc                                                 <  1:629348/150‑1 (MQ=255)                                                                                                       aTACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGAtt                                            <  2:139386/151‑1 (MQ=255)                                                                                                                                 cATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTcagcag                   >  1:578725/1‑150 (MQ=255)                                                                                                                                     aTAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGaa              >  2:660083/1‑151 (MQ=255)                                                                                                                                       aGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGAAtt            >  2:479318/1‑151 (MQ=255)                                                                                                                                             cACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCGCCAGCGTCAGCAGCAGAATTtctc        <  1:586877/149‑1 (MQ=255)                                                                                                                                               cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGAATTTCTCGTTg    >  2:723072/1‑151 (MQ=255)                                                                                                                                               cTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCAACCAGCGTCAGCAGCAGAATTTCTCGTTg    >  2:758918/1‑151 (MQ=255)                                                                                                                                                 tGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGAATTTCTCGTTgcg  >  2:495111/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACAGTACCGATTTGAGTTTTTCCACCAGCGTCAGCAGCAGAATTTCTCGTTGCG  >  NC_000913/3684456‑3684750

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