Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,138	(A)8→7	intergenic (‑123/‑122)	yibN ← / → gpmM	putative sulfurtransferase YibN/2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,785,131	0	A	.	100.0%	142.6 / NA	34	intergenic (‑116/‑129)	yibN/gpmM	putative sulfurtransferase YibN/2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (15/19);  total (15/19)

AAAGGAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGAAAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATT  >  NC_000913/3785057‑3785207                                                                           |                                                                             aaaGGAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGaaaaa                                                                          >  2:676740/1‑79 (MQ=255)   aGGAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGaaaaaaagc                                                                      <  1:676698/81‑3 (MQ=255)   aGGAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGaaaaaaagc                                                                      <  2:676698/81‑3 (MQ=255)    ggAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGaaaaaaagcg                                                                     >  1:676738/1‑78 (MQ=255)       gTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTaa                                                                 <  2:676699/81‑1 (MQ=255)         ataCCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAAt                                                               >  2:676741/1‑81 (MQ=255)             cTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTgg                                                           >  2:676744/1‑81 (MQ=255)              tGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTgg                                                           >  2:676747/1‑80 (MQ=255)               gcgcGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCa                                                         >  1:676740/1‑81 (MQ=255)                cgcgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGc                                                          <  1:676701/79‑1 (MQ=255)                cgcgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCAt                                                        <  2:676703/81‑1 (MQ=255)                cgcgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCAt                                                        <  2:676702/81‑1 (MQ=255)                cgcgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCAt                                                        <  2:676701/81‑1 (MQ=255)                 gcgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATc                                                       >  1:676755/1‑81 (MQ=255)                  cgGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCg                                                      <  2:676704/81‑1 (MQ=255)                   ggTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCgg                                                     <  1:676706/81‑1 (MQ=255)                     tGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCt                                                   <  1:676708/81‑1 (MQ=255)                     tGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCt                                                   <  2:676706/81‑1 (MQ=255)                       cGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCtt                                                  <  2:676708/80‑1 (MQ=255)                         gCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGcc                                               <  1:676710/81‑1 (MQ=255)                             aTACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCcgc                                             >  1:676743/1‑79 (MQ=255)                              tACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAg                                          <  2:676710/81‑1 (MQ=255)                                  gccagcGCGTTAACTGAAATGCAATTTAGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGt                                      <  2:676711/81‑1 (MQ=255)                                    cagcGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCg                                    <  1:676713/81‑1 (MQ=255)                                        gcgTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGtata                                >  2:676764/1‑81 (MQ=255)                                          gTTAACTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATAtt                              >  1:676756/1‑81 (MQ=255)                                             aaCTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATAtttt                            <  1:676714/80‑1 (MQ=255)                                             aaCTGGAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATAtttt                            >  2:676746/1‑80 (MQ=255)                                                  gAATGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCt                       <  1:676715/80‑1 (MQ=255)                                                    aTGCAATTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCtt                      >  2:676757/1‑79 (MQ=255)                                                        aaTTTTGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAgg                >  2:676758/1‑81 (MQ=255)                                                             tGCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTAtt           >  2:676760/1‑81 (MQ=255)                                                              gCGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTAttt          >  2:676750/1‑81 (MQ=255)                                                               cGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGGGCGGTCGTATATTTTGATCTTCAGAGGCTAtttt         <  2:676716/81‑1 (MQ=255)                                                               cGGGGCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTAtttt         <  1:676716/81‑1 (MQ=255)                                                                ggggCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTAttttt        <  2:676717/81‑2 (MQ=255)                                                                ggggCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTa        >  1:676767/1‑81 (MQ=255)                                                                 gggCGCGACG‑AAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTAt       >  2:676759/1‑81 (MQ=255)                                                                      cgacg‑aAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGAtt  >  2:676753/1‑81 (MQ=255)                                                                           |                                                                             AAAGGAGTATACCTGCGCGGTGCGGCAAATACAGCCAGCGCGTTAACTGGAATGCAATTTTGCGGGGCGCGACGAAAAAAAAGCGTAAATCTGGCATCGGCTTGCCGCGAGCGGTCGTATATTTTGATCTTCAGAGGCTATTTTATCGATT  >  NC_000913/3785057‑3785207

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