Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,398	T→G	L29R (CTT→CGT)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,398	0	T	G	100.0%	89.1 / NA	31	L29R (CTT→CGT)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/17);  total (14/17)

CATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACG  >  NC_000913/3041271‑3041533                                                                                                                                |                                                                                                                                        catctgtcagtcaTATGCCGTTAAGCCGCTCCGGACTTTATTATGGCTTTTACATTTTTTTCTCCACGTCAGCGTTCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGGTTGATGACTGGGCGCGTGCCACCATTACTgg                                                                                                                          <  2:184298/138‑1 (MQ=255)   tcaGTCATTGATACGCCGTTAAGCAACTCAGGACGCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGgcg                                                                                                                        <  1:89882/141‑1 (MQ=255)           tAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGc                                                                                                                <  2:213814/143‑1 (MQ=255)            aaTACGCCGTTAAGCAACCCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCCTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCg                                                                                                               <  2:469749/143‑1 (MQ=255)                                aGGACTCGAGTCTGGCTTTTACCCCTTTTCCTCCCCGTCAGCCTACGGCTGCTGCCCGTTCACCGCTGACGCTGATGACGCTTGATGCCTGGGCGCGTGCCGCCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAg                                                                                           <  1:68008/143‑1 (MQ=255)                                    ctctATTATGGCTTGTACCCCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTGGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGa                                                                                       <  2:157623/143‑1 (MQ=255)                                     tctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTATCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATAGAGCCAGGGTCAGGTGAc                                                                                      >  2:108508/1‑143 (MQ=255)                                           atGGCTTTTACACCTTTGCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGa                                                                                <  1:396716/143‑1 (MQ=255)                                            tGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGAt                                                                               <  1:289203/143‑1 (MQ=255)                                                ttttACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGCCAGCAGTTGTCa                                                                           >  1:115483/1‑143 (MQ=255)                                                        cTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGAAGATGTCAGCCAGATg                                                                   >  2:509743/1‑143 (MQ=255)                                                              ctccCCGTCAGCCTGCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCagaa                                                             <  2:189712/143‑1 (MQ=255)                                                               tccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGCTGGcagcag                                                            >  2:653451/1‑143 (MQ=255)                                                                ccccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAACATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCagaaga                                                           >  1:543084/1‑143 (MQ=255)                                                                 cccGTCAGCCTACGGCTGCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCCCAATTACTCGCGCGGACAGCGGAAAATATACGCAGGGTCAGGGGACAGCAGATGCCAGCCACATGGAagaagac                                                          >  2:239536/1‑142 (MQ=255)                                                                  ccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATc                                                         <  2:253744/143‑1 (MQ=255)                                                                   cGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCa                                                        <  1:128335/143‑1 (MQ=255)                                                                     tctgCCGACCGCGTCTGCCCGTGTCCCGCTGACGCTGGTGCCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGc                                                      <  2:359445/140‑1 (MQ=255)                                                                          cTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTg                                                 >  1:580104/1‑143 (MQ=255)                                                                               gCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTcgc                                            >  1:430941/1‑143 (MQ=255)                                                                                           ttACAGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACAATGACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGa                                <  1:382778/143‑1 (MQ=255)                                                                                            tACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCAATCGCGAc                               >  2:632169/1‑143 (MQ=255)                                                                                              ccGCTGACGCTGATGACGCTTGATGACTGGGTGCGTGCCACCATTACTGGCGCGGACAGCGAAGAATGTATGGCGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTCCTCGCGGCATATTGCGACGc                             >  2:550640/1‑143 (MQ=255)                                                                                                 ctgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCaa                          >  1:104390/1‑143 (MQ=255)                                                                                                 ctgacgctgaTGACGCTTGACGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGCCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCaa                          <  1:12992/143‑1 (MQ=255)                                                                                                   gacgctgaTGGCGCTGGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAGAAATATATGCAGGCTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGTCGCTCTTTGCGACGCCAAAg                        <  2:443264/143‑1 (MQ=255)                                                                                                         gatgaCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGAAAGCGAGAAATACATGCAGGGGAAGGTGATGGCAGCTGTCAGCCAGATGGCAGACGATCAGCAACTGCTCGCCGCCCATAACGACGCAAAAGGCaaaa                  >  2:548072/1‑143 (MQ=255)                                                                                                             aCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtgtg              <  1:112834/143‑1 (MQ=255)                                                                                                               gCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGa            >  1:114699/1‑143 (MQ=255)                                                                                                                        cTGGGCGCGTGCCACCGTGACTGGCGCGGAGAGCGAAAAATATATGCAGGGTCAGGTGACGGCAGATGGCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTAc   <  2:591982/143‑1 (MQ=255)                                                                                                                         tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGGCCATTGCGACGCCAAAAGTAAAAAGAGGAGAAATTCACg  >  2:271683/1‑143 (MQ=255)                                                                                                                                |                                                                                                                                        CATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACG  >  NC_000913/3041271‑3041533

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