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%	102.3 / NA	29	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 (10/19);  total (10/19)

CAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGAT  >  NC_000913/4109814‑4109961                                                                         |                                                                            cAAACGCTGATTTCCCCAAGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGtt                                                                     <  1:762221/81‑1 (MQ=255)   aaCGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTAc                                                                   >  1:762265/1‑81 (MQ=255)        tGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTg                                                              >  1:762261/1‑81 (MQ=255)          aTTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGc                                                             <  2:762223/80‑1 (MQ=255)          aTTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGc                                                             <  2:762222/80‑1 (MQ=255)            ttCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTCATTTGTCACGACTTTGCGTTACCCTTGCAtc                                                          >  2:762224/1‑81 (MQ=255)                 cAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCttt                                                     <  2:762225/81‑1 (MQ=255)                   ggATAATTATCAAACCCGGTGGTCTCTCGCGACCGGGTTTCTTATTTGTCACGACTCTGCGTTACCCTTGCATCTCctt                                                     >  2:762281/1‑79 (MQ=255)                              aaaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGa                                         <  1:762226/80‑1 (MQ=255)                               aaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaa                                       <  1:762228/81‑1 (MQ=255)                               aaCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaa                                       <  1:762227/81‑1 (MQ=255)                                aCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaa                                      <  1:762231/81‑1 (MQ=255)                                aCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaa                                      <  1:762229/81‑1 (MQ=255)                                aCCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGaaaa                                      <  2:762230/81‑1 (MQ=255)                                      tggtTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGAt                                >  1:762266/1‑81 (MQ=255)                                       ggtTTCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGAt                                >  2:762267/1‑80 (MQ=255)                                         tttCTCGCGACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGa                              >  2:762280/1‑80 (MQ=255)                                               gcgACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAg                        <  2:762232/80‑1 (MQ=255)                                               gcgACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGc                       <  1:762233/81‑1 (MQ=255)                                                cgACCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCg                      >  1:762284/1‑81 (MQ=255)                                                  aCCGGGTTTTTTATTTGTCACGACTTTGCGTTACCCTGGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCggg                    <  2:762234/81‑1 (MQ=255)                                                       gTTTTTTATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGtcttc               <  1:762235/81‑1 (MQ=255)                                                         tttttATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTcttctt             >  1:762273/1‑81 (MQ=255)                                                           tttATTTGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCttt            <  2:762236/80‑1 (MQ=255)                                                               tttGTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTg       <  1:762237/81‑1 (MQ=255)                                                                  gTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGt     <  1:762238/80‑1 (MQ=255)                                                                  gTCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGt     >  2:762282/1‑80 (MQ=255)                                                                   tCACGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtga   <  2:762239/81‑1 (MQ=255)                                                                     aCGACTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGtgat  <  1:762240/80‑1 (MQ=255)                                                                         |                                                                            CAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCTTGCATCTCTTTGAGGTACAGGGAAAAAAAGATGAAAAAAGCGGGTCTTCTTTTTTTGGTGAT  >  NC_000913/4109814‑4109961

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