Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,323,424	G→T	T310N (ACC→AAC)	phnD ←	phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,323,424	0	G	T	100.0%	116.7 / NA	34	T310N (ACC→AAC)	phnD	phosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/18);  total (16/18)

GCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCG  >  NC_000913/4323346‑4323493                                                                               |                                                                      gCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGttt                                                                     >  1:800584/1‑81 (MQ=255)  cTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGtttt                                                                    >  2:793431/1‑81 (MQ=255)        tCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGc                                                                >  2:793407/1‑79 (MQ=255)         cACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAg                                                             >  1:793410/1‑81 (MQ=255)           ccGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCt                                                           <  1:793376/81‑1 (MQ=255)           ccGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCt                                                           <  1:793375/81‑1 (MQ=255)            cGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCtt                                                          <  2:793377/81‑1 (MQ=255)             gAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTg                                                         >  2:793378/1‑81 (MQ=255)              aaCTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGt                                                        >  1:793424/1‑81 (MQ=255)                 tCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCt                                                     >  2:793417/1‑81 (MQ=255)                   aTCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCt                                                     <  1:793379/79‑1 (MQ=255)                    tCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCt                                                  <  2:793380/81‑1 (MQ=255)                     cgcgCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCt                                                  >  1:793423/1‑80 (MQ=255)                      gcgcTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCa                                                <  2:793381/81‑1 (MQ=255)                        gcTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCAtt                                              >  2:793419/1‑81 (MQ=255)                             gcgcGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTc                                         >  1:793418/1‑81 (MQ=255)                              cgcgTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTcc                                        <  1:793382/81‑1 (MQ=255)                                 gttgttCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCttt                                     <  2:793383/81‑1 (MQ=255)                                     ttCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTg                                 >  2:793422/1‑81 (MQ=255)                                       cAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCt                               >  2:793416/1‑81 (MQ=255)                                        aGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCtt                              <  2:793384/81‑1 (MQ=255)                                        aGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCtt                              <  1:793385/81‑1 (MQ=255)                                               ccAGGTCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGc                       >  1:793421/1‑81 (MQ=255)                                                    tCATCCAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGc                  <  2:793386/81‑1 (MQ=255)                                                         cAGCTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCAtct              <  1:793387/80‑1 (MQ=255)                                                            cTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCttt          <  2:793388/81‑1 (MQ=255)                                                            cTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCttt          <  2:793389/81‑1 (MQ=255)                                                            cTGCGCCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCttt          <  2:793390/81‑1 (MQ=255)                                                               cgcCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTaaa       >  2:793429/1‑81 (MQ=255)                                                               cgcCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTaaa       >  1:793426/1‑81 (MQ=255)                                                               cgcCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTa         <  2:793391/79‑1 (MQ=255)                                                                gcCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTaaa       >  1:793434/1‑80 (MQ=255)                                                                gcCTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAAc      <  1:793392/81‑1 (MQ=255)                                                                  cTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAg    <  1:793393/81‑1 (MQ=255)                                                                  cTGAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAg    >  1:793433/1‑81 (MQ=255)                                                                    gAATCTCGGTGTTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAgcg  >  1:793428/1‑81 (MQ=255)                                                                               |                                                                      GCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCG  >  NC_000913/4323346‑4323493

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