Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,850	C→T	A36T (GCG→ACG)	yibF ←	glutathione S‑transferase homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,850	0	C	T	100.0%	133.7 / NA	42	A36T (GCG→ACG)	yibF	glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/19);  total (23/19)

CGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGC  >  NC_000913/3761722‑3761971                                                                                                                                 |                                                                                                                          cGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                                 >  1:20965/1‑139 (MQ=255)                    aaTTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                              <  2:1506227/139‑1 (MQ=255)                           tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                              >  2:477213/1‑132 (MQ=255)                           tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                              <  1:477213/132‑1 (MQ=255)                           tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTAt                                                                                       >  1:60362/1‑139 (MQ=255)                                tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                  >  1:1274012/1‑139 (MQ=255)                                     gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAGCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                             <  1:311007/139‑1 (MQ=255)                                       tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                           >  2:310290/1‑139 (MQ=255)                                         atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                  <  2:334577/98‑1 (MQ=255)                                         atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                  >  1:334577/1‑98 (MQ=255)                                         atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                         >  1:952784/1‑139 (MQ=255)                                                           cATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGc                                                       >  2:371229/1‑139 (MQ=255)                                                              tCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTa                                                    >  1:327500/1‑139 (MQ=255)                                                               cGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                  <  1:1158278/108‑1 (MQ=255)                                                               cGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                  >  2:1158278/1‑108 (MQ=255)                                                                  ccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaaa                                                <  1:371229/139‑1 (MQ=255)                                                                  ccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaaa                                                <  2:60362/139‑1 (MQ=255)                                                                    tctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACg                                              <  2:1274012/139‑1 (MQ=255)                                                                      ttCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                            <  1:310290/139‑1 (MQ=255)                                                                                aGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGgtg                                      <  1:729098/135‑1 (MQ=255)                                                                                aGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGgtg                                      >  2:729098/1‑135 (MQ=255)                                                                                  cACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTa                                >  1:1663152/1‑139 (MQ=255)                                                                                       gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCa                                                                                                                <  1:91399/54‑1 (MQ=255)                                                                                       gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCa                                                                                                                >  2:91399/1‑54 (MQ=255)                                                                                         aCTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgacga                         >  2:910962/1‑139 (MQ=255)                                                                                                    cGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCg              <  1:910962/139‑1 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                              >  1:425294/1‑56 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                              <  2:425294/56‑1 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                            >  2:1502884/1‑106 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                            <  1:1502884/106‑1 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt           >  2:502747/1‑139 (MQ=255)                                                                                                       gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt           >  1:403371/1‑139 (MQ=255)                                                                                                        ttAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAg                                                                     <  1:1515941/80‑1 (MQ=255)                                                                                                        ttAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAg                                                                     >  2:1515941/1‑80 (MQ=255)                                                                                                         tAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGa                                                               >  1:1333606/1‑85 (MQ=255)                                                                                                         tAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGa                                                               <  2:1333606/85‑1 (MQ=255)                                                                                                         tAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGt                                                           >  1:1088569/1‑89 (MQ=255)                                                                                                         tAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGt                                                           <  2:1088569/89‑1 (MQ=255)                                                                                                          aaaTTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCt        <  2:327500/139‑1 (MQ=255)                                                                                                                cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGc  >  2:228133/1‑139 (MQ=255)                                                                                                                     cGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                            <  2:260044/92‑1 (MQ=255)                                                                                                                     cGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                            >  1:260044/1‑92 (MQ=255)                                                                                                                                 |                                                                                                                          CGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGC  >  NC_000913/3761722‑3761971

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