Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,109,886	T→C	intergenic (+63/‑44)	sbp → / → cdh	sulfate/thiosulfate ABC transporter periplasmic binding protein Sbp/CDP‑diacylglycerol diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,109,886	0	T	C	100.0%	103.4 / NA	30	intergenic (+63/‑44)	sbp/cdh	sulfate/thiosulfate ABC transporter periplasmic binding protein Sbp/CDP‑diacylglycerol diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/18);  total (12/18)

AGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGAT  >  NC_000913/4109812‑4109961                                                                           |                                                                            agCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCg                                                                       <  2:831604/81‑1 (MQ=255)  gCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGt                                                                      >  1:831640/1‑81 (MQ=255)     aaCGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTAc                                                                   <  1:831605/81‑1 (MQ=255)      aCGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTAcc                                                                  <  1:831606/81‑1 (MQ=255)         cTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCtt                                                               <  1:831607/81‑1 (MQ=255)                 cccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATctct                                                       >  2:831645/1‑81 (MQ=255)                 cccAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATctct                                                       >  2:831641/1‑81 (MQ=255)                   cAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCACCTCttt                                                     <  2:831608/81‑1 (MQ=255)                    aGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCttt                                                     >  1:831649/1‑80 (MQ=255)                             aTCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAgg                                           >  1:831658/1‑81 (MQ=255)                                aaaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaa                                        <  2:831610/81‑1 (MQ=255)                                aaaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGa                                         <  1:831609/80‑1 (MQ=255)                                 aaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaa                                       <  2:831611/81‑1 (MQ=255)                                  aCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaa                                       >  2:831650/1‑80 (MQ=255)                                    ccGGTGGTTTCTCGCGACCGGTTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaaaa                                    <  1:831612/81‑1 (MQ=255)                                           tttCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGaa                             <  2:831613/81‑1 (MQ=255)                                           tttCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGa                              >  2:831659/1‑80 (MQ=255)                                                 gcgACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGaaaaa                          <  2:831614/78‑1 (MQ=255)                                                 gcgACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGc                       >  2:831652/1‑81 (MQ=255)                                                     ccGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGt                   <  1:831615/81‑1 (MQ=255)                                                     ccGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGt                   <  2:831617/81‑1 (MQ=255)                                                     ccGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGt                   <  1:831616/81‑1 (MQ=255)                                                      cGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGtc                  <  2:831618/81‑1 (MQ=255)                                                       gggTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGtct                 <  1:831619/81‑1 (MQ=255)                                                       gggTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAACGGGtct                 >  1:831655/1‑81 (MQ=255)                                                          ttttttATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGtcttct              <  1:831620/81‑1 (MQ=255)                                                              ttATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCttttt          >  2:831654/1‑81 (MQ=255)                                                               tATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCtttttt         >  2:831651/1‑81 (MQ=255)                                                                   tGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGt     >  2:831653/1‑81 (MQ=255)                                                                       aCGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtgat  <  1:831621/80‑1 (MQ=255)                                                                           |                                                                            AGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGAT  >  NC_000913/4109812‑4109961

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