Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,081	A→C	W83G (TGG→GGG)	yieP ←	putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,081	0	A	C	100.0%	186.9 / NA	51	W83G (TGG→GGG)	yieP	putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/25);  total (26/25)

GCTTTCTGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCGGGCAA  >  NC_000913/3940933‑3941219                                                                                                                                                     |                                                                                                                                           gcTTTCTGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGt                                                                                                                                          >  2:632429/1‑151 (MQ=255) gcTTTCTGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGt                                                                                                                                          >  1:554804/1‑151 (MQ=255)     tCTGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCg                                                                                                                                       >  1:672639/1‑150 (MQ=255)       tGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGAtt                                                                                                                                    >  1:299334/1‑151 (MQ=255)        gTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGAt                                                                                                                                     <  1:210464/149‑1 (MQ=255)         ttCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATtg                                                                                                                                   <  1:839664/150‑1 (MQ=255)           ccGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTgg                                                                                                                                >  1:125009/1‑151 (MQ=255)           ccGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTgg                                                                                                                                <  2:32478/151‑1 (MQ=255)            cgcgGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGc                                                                                                                               <  2:658515/151‑1 (MQ=255)              cgGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGATCGATTGTGGCAt                                                                                                                             <  1:182027/151‑1 (MQ=255)                gTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGa                                                                                                                           <  2:552517/151‑1 (MQ=255)                     aaCCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCg                                                                                                                       >  2:579153/1‑150 (MQ=255)                       ccGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGt                                                                                                                    <  1:626503/151‑1 (MQ=255)                       ccGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGt                                                                                                                    <  1:675690/151‑1 (MQ=255)                        cGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGt                                                                                                                    >  1:776376/1‑150 (MQ=255)                               gcCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAAtt                                                                                                             <  1:287874/150‑1 (MQ=255)                               gcCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATtc                                                                                                            <  2:83094/151‑1 (MQ=255)                                     aggcaggcTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATtcgcggtc                                                                                                      >  1:702353/1‑151 (MQ=255)                                           gcTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTa                                                                                                <  2:428403/151‑1 (MQ=255)                                           gcTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTa                                                                                                <  2:163193/151‑1 (MQ=255)                                           gcTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTa                                                                                                <  2:719053/151‑1 (MQ=255)                                             ttGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTaa                                                                                               >  1:379230/1‑150 (MQ=255)                                                   cTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATc                                                                                        >  1:522987/1‑151 (MQ=255)                                                     ccAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATcc                                                                                       <  2:125758/150‑1 (MQ=255)                                                           agaTACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGc                                                                                <  1:760718/151‑1 (MQ=255)                                                             aTACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCg                                                                              >  2:895101/1‑151 (MQ=255)                                                                    tAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGt                                                                       >  1:754639/1‑151 (MQ=255)                                                                      aCCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGttt                                                                     >  2:145368/1‑151 (MQ=255)                                                                                gATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTc                                                           >  2:901336/1‑151 (MQ=255)                                                                                   cGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTGTTGCCGTTAACGTTTTGACCGCTTcgcg                                                        >  1:726266/1‑151 (MQ=255)                                                                                      tGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                                                     >  1:294429/1‑151 (MQ=255)                                                                                       gACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                                                     <  2:74081/150‑1 (MQ=255)                                                                                          cTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCt                                                 >  1:815318/1‑151 (MQ=255)                                                                                             atgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTa                                              <  2:921870/151‑1 (MQ=255)                                                                                               gaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTAc                                             <  1:697374/150‑1 (MQ=255)                                                                                                aaaGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACg                                            >  1:11083/1‑150 (MQ=255)                                                                                                aaaGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGa                                           <  2:319844/151‑1 (MQ=255)                                                                                                          tctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCa                                  >  2:510666/1‑150 (MQ=255)                                                                                                             gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTCCGACTCACTCCaa                                 <  1:169484/148‑1 (MQ=255)                                                                                                             gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaa                                 <  1:839169/148‑1 (MQ=255)                                                                                                             gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAAtt                              >  2:790243/1‑151 (MQ=255)                                                                                                               cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCATACCGCTGTACGACTCACTCCAAATTGc                            <  2:725788/151‑1 (MQ=255)                                                                                                                           aaGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCa                <  2:156668/151‑1 (MQ=255)                                                                                                                             gCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAAt              <  1:91058/151‑1 (MQ=255)                                                                                                                                  tCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGGACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCAc         >  2:573808/1‑151 (MQ=255)                                                                                                                                    cTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCg       >  2:65968/1‑151 (MQ=255)                                                                                                                                      gATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCggg     >  2:781351/1‑151 (MQ=255)                                                                                                                                      gATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCggg     >  1:614828/1‑151 (MQ=255)                                                                                                                                      gATCAAGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCgg      <  2:672658/150‑1 (MQ=255)                                                                                                                                          aaGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCGGGCaa  >  1:925694/1‑150 (MQ=255)                                                                                                                                          aaGAAAATTCCCGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCGGGCaa  >  2:452161/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                           GCTTTCTGTTCCGCGGTGCCAACCGTTGCTGCCAGTAGGCAGGCTTGCGGCTCCAGACAGATACGCATAACCAGAAAGTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAATTGCTCGCCCAGCTCAATCTCACCGGGCAA  >  NC_000913/3940933‑3941219

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