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%	175.2 / NA	54	A36T (GCG→ACG)	yibF	glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/24);  total (30/24)

AATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATA  >  NC_000913/3761741‑3761982                                                                                                              |                                                                                                                                     aaTTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                         <  2:908091/139‑1 (MQ=255)   ttCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGaaa                                                                                                       >  2:608378/1‑139 (MQ=255)        tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                         >  1:961496/1‑132 (MQ=255)        tatatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                         <  2:961496/132‑1 (MQ=255)                gCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCtttt                                                                                          >  2:153182/1‑139 (MQ=255)                 cGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTc                                                                                         >  2:1293549/1‑139 (MQ=255)                  gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                        >  1:241312/1‑139 (MQ=255)                  gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                        <  2:117530/139‑1 (MQ=255)                  gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                        >  2:258728/1‑139 (MQ=255)                    tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCATTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                                      >  1:1055115/1‑139 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                                            <  2:734967/99‑1 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGc                                                                                                                            >  1:734967/1‑99 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                <  1:1207088/127‑1 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                >  2:1207088/1‑127 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                        >  2:1323204/1‑135 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                        <  1:1323204/135‑1 (MQ=255)                      atCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaaca                                                                                    >  1:623854/1‑139 (MQ=255)                       tCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaa                                                                                   <  1:326933/139‑1 (MQ=255)                           cGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGa                                                                               >  2:149667/1‑139 (MQ=255)                           cGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGa                                                                               >  2:194507/1‑139 (MQ=255)                           cGAATCAAACCAGCATTCGCCCTCTTCAGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                                                 <  1:398298/105‑1 (MQ=255)                           cGAATCAAACCAGCATTCGCCCTCTTCAGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                                                 >  2:398298/1‑105 (MQ=255)                                  aaCCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGaaa                                                                        >  1:1094505/1‑139 (MQ=255)                                         aTTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCg                                                                 >  2:358072/1‑139 (MQ=255)                                         aTTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCg                                                                 <  1:153182/139‑1 (MQ=255)                                                   ttCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                                       <  1:194507/139‑1 (MQ=255)                                                     cGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCt                                                     <  1:742805/139‑1 (MQ=255)                                                     cGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCt                                                     <  2:1094505/139‑1 (MQ=255)                                                           ccAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                             >  1:756504/1‑61 (MQ=255)                                                           ccAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                             <  2:756504/61‑1 (MQ=255)                                                               cACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTa                                           >  2:972091/1‑139 (MQ=255)                                                                  cGGCACTTTTCCTAACGGGGTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcg                                        <  1:1343120/139‑1 (MQ=255)                                                                    gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                >  1:661101/1‑81 (MQ=255)                                                                    gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                <  2:661101/81‑1 (MQ=255)                                                                    gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                                       <  1:16873/122‑1 (MQ=255)                                                                    gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                                       >  2:16873/1‑122 (MQ=255)                                                                    gCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACcgac                                      >  1:406021/1‑139 (MQ=255)                                                                         tttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAGCGGGCTGGTGTAGCTACCGACGAGtt                                 >  2:1421854/1‑139 (MQ=255)                                                                         tttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGtt                                 >  2:919434/1‑139 (MQ=255)                                                                               aaCGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGc                           <  2:623854/139‑1 (MQ=255)                                                                                  gggTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                             <  2:1287710/70‑1 (MQ=255)                                                                                  gggTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                             >  1:1287710/1‑70 (MQ=255)                                                                                  gggTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGa                        <  1:608378/139‑1 (MQ=255)                                                                                    gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                         <  1:1245039/56‑1 (MQ=255)                                                                                    gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGa                                                                                                         >  2:1245039/1‑56 (MQ=255)                                                                                    gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt                      >  1:1073524/1‑139 (MQ=255)                                                                                    gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt                      >  2:1401531/1‑139 (MQ=255)                                                                                            gcgcCACGCCATTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTtgtg              <  2:1055115/139‑1 (MQ=255)                                                                                             cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                             >  2:1160146/1‑59 (MQ=255)                                                                                             cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                             <  1:1160146/59‑1 (MQ=255)                                                                                             cgcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGc             >  2:1192167/1‑139 (MQ=255)                                                                                              gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCg            >  2:961031/1‑139 (MQ=255)                                                                                                        tGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAata  <  1:149667/139‑1 (MQ=255)                                                                                                        tGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAata  <  1:1401531/139‑1 (MQ=255)                                                                                                              |                                                                                                                                     AATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATA  >  NC_000913/3761741‑3761982

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