Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	702,903	T→C	M237V (ATG→GTG)	nagB ←	glucosamine‑6‑phosphate deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	702,903	0	T	C	100.0%	119.9 / NA	32	M237V (ATG→GTG)	nagB	glucosamine‑6‑phosphate deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/10);  total (22/10)

ATTGCCCTGAGCAAGGAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCATGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCCAGCAGTGTACCAA  >  NC_000913/702775‑703043                                                                                                                                 |                                                                                                                                             atTGCCCTGAGCAAGGAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTg                                                                                                                         >  2:169045/1‑150 (MQ=255)  ttGCCCTGAGCAAGGAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGc                                                                                                                        >  1:766337/1‑150 (MQ=255)             aaGGAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATg                                                                                                            >  1:78855/1‑151 (MQ=255)                gAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTc                                                                                                         <  2:560976/151‑1 (MQ=255)                     agggcagggATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACat                                                                                                    <  1:703003/151‑1 (MQ=255)                          agggATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGt                                                                                               >  2:141730/1‑151 (MQ=255)                           gggATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTGTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGtt                                                                                              <  2:196324/151‑1 (MQ=255)                              aTAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCAc                                                                                           >  1:523350/1‑151 (MQ=255)                                   aaTTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAAc                                                                                      >  1:627918/1‑151 (MQ=255)                                                tGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACggcggc                                                                         >  1:104032/1‑151 (MQ=255)                                                 gATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACggcggc                                                                         <  1:598786/150‑1 (MQ=255)                                                 gATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGcc                                                                        >  2:270080/1‑151 (MQ=255)                                                   tatTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCt                                                                       <  2:282087/150‑1 (MQ=255)                                                         cTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAgcgc                                                                >  2:247037/1‑151 (MQ=255)                                                                     cATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTg                                                    >  1:293283/1‑151 (MQ=255)                                                                       ttGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGc                                                  >  2:581821/1‑151 (MQ=255)                                                                       ttGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGc                                                  >  1:715354/1‑151 (MQ=255)                                                                          aaaTATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTAc                                               >  1:743355/1‑151 (MQ=255)                                                                           aaTATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTAcc                                              >  1:389239/1‑151 (MQ=255)                                                                                  tAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCAcc                                       >  2:303847/1‑151 (MQ=255)                                                                                      gTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGa                                    >  1:103161/1‑150 (MQ=255)                                                                                           aaCTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAatcatc                              >  1:567763/1‑151 (MQ=255)                                                                                            aCTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAatcatc                              >  1:20629/1‑150 (MQ=255)                                                                                              tttCAGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCAct                           >  1:212537/1‑151 (MQ=255)                                                                                                  aGCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCActtctt                       <  2:235883/151‑1 (MQ=255)                                                                                                   gCTCCATGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACttcttc                      <  1:311772/151‑1 (MQ=255)                                                                                                        aTGGTGGAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCAt                 <  2:627881/151‑1 (MQ=255)                                                                                                         tggtggAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCAAGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATc                >  2:505352/1‑151 (MQ=255)                                                                                                           gtggAAGGTTCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCca              <  2:158588/151‑1 (MQ=255)                                                                                                                   ttCATCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCCAGCAgtgt       >  1:40363/1‑150 (MQ=255)                                                                                                                       tCGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCCAGCAGTGTACCa   >  2:183525/1‑150 (MQ=255)                                                                                                                        cGCACACCACGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCCAGCAGTGTACCaa  <  2:104030/150‑1 (MQ=255)                                                                                                                                 |                                                                                                                                             ATTGCCCTGAGCAAGGAGCCAGGGCAGGGATAACAATTACAGACCTTTGATATTTTCTGCTTCTAATTCATTGAAATATCTTAAAGTCTTAACTTTCAGCTCCATGGTGGAAGGTTCATCGCACACCATGATCGCTTTCGGATGCAGTTGCAGACAGCTGATGGTCCACATATGGTTCACGCAACCTTCAACGGCGGCCTGCAGCGCCAGTGCTTTCTGGCTACCCAGCACCAGAATCATCACTTCTTCGGCATCCAGCAGTGTACCAA  >  NC_000913/702775‑703043

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