Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,318,237	T→A	F124Y (TTT→TAT)	miaB →	2‑methylthioadenine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,318,237	0	T	A	100.0%	87.2 / NA	28	F124Y (TTT→TAT)	miaB	2‑methylthioadenine synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (20/8);  total (20/8)

ATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAG  >  CP000730/1318111‑1318362                                                                                                                               |                                                                                                                              aTGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTa                                                                                                                 >  1:340692/1‑141 (MQ=255)     cctcagtgTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAAttt                                                                                                             >  1:410663/8‑141 (MQ=255)        tGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGaa                                                                                                          >  2:199688/1‑141 (MQ=255)          aGGTCATTGCTGGTATACTTGAAGCATTAGGCTATAAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGc                                                                                                        >  2:64962/1‑141 (MQ=255)           ggTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCa                                                                                                       <  2:242489/141‑1 (MQ=255)            gTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATACTGAAATAGGTAATTTGAAGCAt                                                                                                      >  1:344824/1‑141 (MQ=255)               aTTGCTTGTATACTTGAAGCATGAGGCTATCAAGCAACGACTTATATTAACACTTCAGATGTTATTTTAATTAATACATCTGCAATTAGAGAAAATGTCGAGAAAAAAGTGTATAGTGAAATAGGTAATGTGAAGCATTTg                                                                                                   >  2:55378/1‑141 (MQ=255)                     ggTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGaaaaaa                                                                                             <  2:283043/141‑1 (MQ=255)                        ataCTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGaa                                                                                          >  2:303409/1‑141 (MQ=255)                                aGCATTAGGCTAAAAAGAAAAGACTAATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAAAACAAAGTGTATAGTGAAATAGGTAATTAAAAAAAAAAAAAAAAAGAACGACCTGa                                                                                  <  2:124414/141‑1 (MQ=255)                                 gCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAAATTGAAGCATTTGAAAAAAGAACGACCTGat                                                                                 >  1:191868/1‑141 (MQ=255)                                     tAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATAttt                                                                             >  1:239182/1‑141 (MQ=255)                                      aGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATAtttt                                                                            <  2:189887/141‑1 (MQ=255)                                          tATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATc                                                                        >  2:378303/1‑141 (MQ=255)                                            tCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCgg                                                                      <  2:331870/141‑1 (MQ=255)                                             cAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGgt                                                                     >  2:240714/1‑141 (MQ=255)                                             cAAGCAACGACTGATATTAAAACTGCAAATGTTATCTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTAAAGTGAAATAGGTAATTTGAAGCATTAGAAAAAATAACGACCTGATATTTTAATCgat                                                                     >  2:394273/1‑139 (MQ=255)                                                  aaCGACTTATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGAATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtg                                                                <  2:259935/141‑1 (MQ=255)                                                     gACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTgg                                                             >  1:399728/1‑141 (MQ=255)                                                       cTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGtt                                                           >  1:191762/1‑141 (MQ=255)                                                                              aTTTTAATTAATACATGTGCGATTAGAGAAAATCCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACATGATATTTTAATCGGTGCTTGTGCTTGTATGTCACAAGACGAGTCagta                                    >  2:210710/1‑141 (MQ=255)                                                                               ttttAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCagtag                                   <  2:308718/141‑1 (MQ=255)                                                                                 ttaatgaaTACATGTGCGATTAGAGCAAATGCCGAGAACAACGTGTATAGTGCAATAGGTCATTTGAAGCATTGGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTg                                 <  2:206588/141‑1 (MQ=255)                                                                                     ttaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATa                             >  2:21803/1‑141 (MQ=255)                                                                                              gtgCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTaa                    >  2:296430/1‑141 (MQ=255)                                                                                              gtgCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTaa                    >  2:40426/1‑141 (MQ=255)                                                                                                             gtGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATg     >  1:316454/2‑141 (MQ=255)                                                                                                                ccGAGAACAAAGTGTATAGTGAAATAGGTAATTTGCGGCATTTGAAAAGAGAGCGACCTGATATGTTAATCGGTGTTCGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGGATAAAAGTTTAAAATCGTAtaggtaaattt  >  2:200303/1‑131 (MQ=255)                                                                                                                               |                                                                                                                              ATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAG  >  CP000730/1318111‑1318362

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