Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,152	A→G	I5T (ATC→ACC)	phnC ←	phosphonate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,152	0	A	G	100.0%	203.8 / NA	57	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/28);  total (29/28)

ACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGC  >  NC_000913/4325022‑4325286                                                                                                                                   |                                                                                                                                        aCGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGtt                                                                                                                                 >  1:112295/1‑139 (MQ=255)  cGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGttt                                                                                                                                >  1:121976/1‑139 (MQ=255)  cGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGttt                                                                                                                                >  1:429786/1‑139 (MQ=255)   gTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTg                                                                                                                               <  1:1610548/139‑1 (MQ=255)            tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                                  >  2:1351050/1‑127 (MQ=255)            tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                                  <  1:1351050/127‑1 (MQ=255)              gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTg                                                                                                                               >  2:744981/1‑128 (MQ=255)              gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTg                                                                                                                               <  1:744981/128‑1 (MQ=255)              gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    <  1:1028313/139‑1 (MQ=255)              gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    <  2:112295/139‑1 (MQ=255)              gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    >  2:1028313/1‑139 (MQ=255)                    ccGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCt                                                                                                              >  1:6331/1‑139 (MQ=255)                           ccGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAg                                                                                                      <  2:6331/139‑1 (MQ=255)                             gACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTcc                                                                                                                >  2:914928/1‑128 (MQ=255)                             gACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTcc                                                                                                                <  1:914928/128‑1 (MQ=255)                             gACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTg                                                                                                    >  2:1689890/1‑139 (MQ=255)                                            gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCgg                                                                                     >  1:164300/1‑139 (MQ=255)                                                                aaTGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    <  2:580587/89‑1 (MQ=255)                                                                aaTGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    >  1:580587/1‑89 (MQ=255)                                                                    ttCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAATGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                              <  1:760677/75‑1 (MQ=255)                                                                    ttCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                              >  2:760677/1‑75 (MQ=255)                                                                     tCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaatta                                                            >  1:555865/1‑139 (MQ=255)                                                                       aGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGa                                                                        >  2:781020/1‑125 (MQ=255)                                                                       aGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGa                                                                        <  1:781020/125‑1 (MQ=255)                                                                        gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTaaaa                                                         >  1:24811/1‑139 (MQ=255)                                                                         aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    >  2:1605225/1‑80 (MQ=255)                                                                         aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    <  1:1605225/80‑1 (MQ=255)                                                                            aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                                   <  1:1647498/125‑1 (MQ=255)                                                                            aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                                   >  2:1647498/1‑125 (MQ=255)                                                                              ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                                   <  2:555865/139‑1 (MQ=255)                                                                              ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                                   <  2:24811/139‑1 (MQ=255)                                                                                gCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                                 <  1:1689890/139‑1 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    <  1:1395611/70‑1 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                                    >  2:1395611/1‑70 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTg                                                                                         <  1:1304689/96‑1 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTg                                                                                         >  2:1304689/1‑96 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                                 <  2:445986/136‑1 (MQ=255)                                                                                   tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                                 >  1:445986/1‑136 (MQ=255)                                                                                        cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaa                                                               >  1:1201634/1‑117 (MQ=255)                                                                                        cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaa                                                               <  2:1201634/117‑1 (MQ=255)                                                                                        cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                                         <  1:556448/139‑1 (MQ=255)                                                                                          ccTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGAGACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaat                                                              <  2:1018988/116‑1 (MQ=255)                                                                                          ccTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAATGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaat                                                              >  1:1018988/1‑116 (MQ=255)                                                                                             gATGCTGATTGAAGGTTTTGGCGAGCTTGTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                                   >  2:786423/1‑124 (MQ=255)                                                                                             gATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                                   <  1:786423/124‑1 (MQ=255)                                                                                                   gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg                              <  1:1532878/139‑1 (MQ=255)                                                                                                     ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGaa                            >  2:1376020/1‑139 (MQ=255)                                                                                                             tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa                    <  2:164300/139‑1 (MQ=255)                                                                                                              ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCat                                                                                                                            >  2:1246659/1‑35 (MQ=255)                                                                                                              ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCat                                                                                                                            <  1:1246659/35‑1 (MQ=255)                                                                                                              ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                                   >  2:133732/1‑92 (MQ=255)                                                                                                              ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                                   <  1:133732/92‑1 (MQ=255)                                                                                                              ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGAc                   >  1:576137/1‑139 (MQ=255)                                                                                                               tGGCGAGCTTCTCGACACGGGTAATCGTTTGTATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCg                                                                                      <  2:1296158/71‑1 (MQ=255)                                                                                                               tGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCg                                                                                      >  1:1296158/1‑71 (MQ=255)                                                                                                                      cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCggg           >  1:348179/1‑139 (MQ=255)                                                                                                                               acGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAgggcgggc  <  1:1376020/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        ACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGC  >  NC_000913/4325022‑4325286

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