Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,949,783	(T)7→6	coding (187/1521 nt)	yifB ←	putative magnesium chelatase YifB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,949,777	0	T	.	100.0%	141.9 / NA	35	coding (193/1521 nt)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/20);  total (15/20)

TCTCCGACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACA  >  NC_000913/3949644‑3949910                                                                                                                                      |                                                                                                                                      tctcCGACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttc                                                                                                                                 <  1:370517/140‑2 (MQ=255)  ctcCGACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttcg                                                                                                                                <  1:422009/140‑3 (MQ=255)   tcCGACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGAttttttcgc                                                                                                                               <  2:458724/140‑4 (MQ=255)       aCTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGa                                                                                                                          <  2:8248/140‑1 (MQ=255)                tattcatCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCatat                                                                                                                 <  2:193212/140‑1 (MQ=255)                     atCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCt                                                                                                            >  1:60440/1‑140 (MQ=255)                       cTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTAt                                                                                                          >  2:99744/1‑140 (MQ=255)                       cTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTAt                                                                                                          <  2:368932/140‑1 (MQ=255)                                      tAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTgc                                                                                           >  2:60179/1‑140 (MQ=255)                                                  aGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATcgc                                                                               <  2:428777/140‑1 (MQ=255)                                                     ccgccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAg                                                                            <  1:60179/140‑1 (MQ=255)                                                     ccgccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAg                                                                            <  1:99744/140‑1 (MQ=255)                                                       gccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGct                                                                          >  2:392800/1‑140 (MQ=255)                                                                 cAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCg                                                                <  1:416648/140‑1 (MQ=255)                                                                 cAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCg                                                                <  2:60440/140‑1 (MQ=255)                                                                                aaTCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGc                                                 >  1:288835/1‑140 (MQ=255)                                                                                      atCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATGATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccacca                                           >  2:314073/1‑140 (MQ=255)                                                                                       tCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCAt                                          >  1:281887/1‑140 (MQ=255)                                                                                        cGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATc                                         >  2:449660/1‑140 (MQ=255)                                                                                          taccccTTCTTTGGGCAGATCAGCTGGATCTAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGGGATCGCTAGGTTCTTTTACCGATGTTTCTGGTAAGACCACCATCGt                                       >  2:122602/3‑140 (MQ=255)                                                                                          tCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGt                                       >  2:176116/1‑140 (MQ=255)                                                                                            ccccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccac                                             <  1:176351/132‑1 (MQ=255)                                                                                            ccccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccac                                             >  2:176351/1‑132 (MQ=255)                                                                                                   ttttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCggg                              <  2:288835/140‑1 (MQ=255)                                                                                                     ttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccac                                             <  1:17095/123‑1 (MQ=255)                                                                                                     ttGGCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccac                                             >  2:17095/1‑123 (MQ=255)                                                                                                        gCAGATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATGATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                         <  1:314073/140‑1 (MQ=255)                                                                                                           gATCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCtt                      <  1:387951/140‑1 (MQ=255)                                                                                                             tCAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTa                    <  2:329136/140‑1 (MQ=255)                                                                                                              cAGCTGGAGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTAc                   <  1:460913/140‑1 (MQ=255)                                                                                                                     aGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATg            <  1:176116/140‑1 (MQ=255)                                                                                                                     aGCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATg            <  2:281887/140‑1 (MQ=255)                                                                                                                      gCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATg            <  2:146197/139‑1 (MQ=255)                                                                                                                      gCCAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATg            >  1:146197/1‑139 (MQ=255)                                                                                                                       ccAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTa          >  2:90479/1‑140 (MQ=255)                                                                                                                        cAGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTAc         <  2:446538/140‑1 (MQ=255)                                                                                                                         aGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTg                 >  1:288593/1‑131 (MQ=255)                                                                                                                         aGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTg                 <  2:288593/131‑1 (MQ=255)                                                                                                                         aGGTTGATGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTAcc        >  2:279274/1‑140 (MQ=255)                                                                                                                               aTGGTGA‑TTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGTTTCTTTTACAGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACa  >  2:174192/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      TCTCCGACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACA  >  NC_000913/3949644‑3949910

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