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

ACATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATAT  >  NC_000913/3761718‑3761983                                                                                                                                     |                                                                                                                                      aCATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTtata                                                                                                                                 <  1:10685/139‑1 (MQ=255)    tCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAggg                                                                                                                              <  2:467987/139‑1 (MQ=255)       ccGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAg                                                                                                                           <  1:83447/139‑1 (MQ=255)         ggAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGtt                                                                                                                         <  1:401287/139‑1 (MQ=255)             cGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCAtt                                                                                                                     >  1:446202/1‑139 (MQ=255)                cattcattAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTACATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                                                  <  1:67190/139‑1 (MQ=255)                cattcattAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaat                                                                                                                  <  2:402537/139‑1 (MQ=255)                       tAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCg                                                                                                           >  1:251264/1‑139 (MQ=255)                                atatTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATg                                                                                                  >  1:405416/1‑139 (MQ=255)                                    tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaata                                                                                                                 >  2:502469/1‑120 (MQ=255)                                    tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATtaata                                                                                                                 <  1:502469/120‑1 (MQ=255)                                    tCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              >  1:176120/1‑139 (MQ=255)                                      aGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCttt                                                                                            >  1:671146/1‑139 (MQ=255)                                       gCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCtttt                                                                                           >  2:303790/1‑139 (MQ=255)                                         gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  1:608082/1‑139 (MQ=255)                                         gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         <  2:671146/139‑1 (MQ=255)                                         gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         <  2:446202/139‑1 (MQ=255)                                         gatgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCt                                                                                         >  2:406507/1‑139 (MQ=255)                                           tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACTAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                 <  2:507787/129‑1 (MQ=255)                                           tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGc                                                                                                 >  1:507787/1‑129 (MQ=255)                                           tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                                       <  2:251264/139‑1 (MQ=255)                                           tgatCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaa                                                                                       <  2:33620/139‑1 (MQ=255)                                              tCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaa                                                                                    <  1:220714/139‑1 (MQ=255)                                               cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  1:442322/139‑1 (MQ=255)                                               cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  2:405416/139‑1 (MQ=255)                                               cGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTaacaac                                                                                   <  2:176120/139‑1 (MQ=255)                                                ggCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAacaaca                                                                                  >  2:480959/1‑139 (MQ=255)                                                            cAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGtt                                                                      <  2:269746/139‑1 (MQ=255)                                                                    gCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACa                                                              >  1:632405/1‑139 (MQ=255)                                                                     cccTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACaa                                                             <  2:85248/139‑1 (MQ=255)                                                                         ctTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         <  1:544597/139‑1 (MQ=255)                                                                          ttCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                                        <  2:89099/139‑1 (MQ=255)                                                                          ttCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACggg                                                        <  1:85863/139‑1 (MQ=255)                                                                                cACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGgtg                                                  <  2:161348/139‑1 (MQ=255)                                                                                      cACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTa                                            >  2:101603/1‑139 (MQ=255)                                                                                                tttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGtt                                  >  1:282907/1‑139 (MQ=255)                                                                                                tttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGtt                                  >  1:112043/1‑139 (MQ=255)                                                                                                 ttCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGttt                                 >  1:743362/1‑139 (MQ=255)                                                                                                  tCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTc                                <  2:296275/139‑1 (MQ=255)                                                                                                           gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         <  2:326450/105‑1 (MQ=255)                                                                                                           gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACgg                                                         >  1:326450/1‑105 (MQ=255)                                                                                                           gTTAAATTGCGCCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGt                       >  1:259346/1‑139 (MQ=255)                                                                                                                     gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCg             <  1:673844/139‑1 (MQ=255)                                                                                                                     gcCACGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCg             <  1:406507/139‑1 (MQ=255)                                                                                                                         cGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              >  2:393087/1‑54 (MQ=255)                                                                                                                         cGCCGTTGTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCt                                                                                              <  1:393087/54‑1 (MQ=255)                                                                                                                                gTCCGTGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAatat  <  2:259346/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      ACATCGCCGGAGCGACATTCATTAATTCAATATATTCAGCGATGATCGGCGAATCAAACCAGCATTCGCCCTCTTCGGTCACCAGCACCGGCACTTTTCCTAACGGGTTAAATTGCGCCACGCCGTTGTCCGCGTTATAGGGCAGTTCATTAATAAATTCGAAAGTTATGCCCTTTTCTAACAACAGAATAGAAAGTTTGCGTACAAACGGGCTGGTGTAGCTACCGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATAT  >  NC_000913/3761718‑3761983

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