Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,704,816	A→T	D325E (GAT→GAA)	dppB ←	dipeptide ABC transporter membrane subunit DppB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,704,816	0	A	T	100.0%	151.0 / NA	42	D325E (GAT→GAA)	dppB	dipeptide ABC transporter membrane subunit DppB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/21);  total (21/21)

TATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGC  >  NC_000913/3704740‑3704892                                                                             |                                                                             tATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTc                                                                            <  2:655913/79‑1 (MQ=255) tATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCg                                                                           <  2:655912/80‑1 (MQ=255) tATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGa                                                                          >  2:655975/1‑81 (MQ=255)    tttCAGTAAGCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGAc                                                                         >  2:655962/1‑79 (MQ=255)          tAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGt                                                                 >  2:655973/1‑81 (MQ=255)            aCCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGtt                                                                >  1:655959/1‑80 (MQ=255)                   aCATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAgg                                                         <  2:655914/80‑1 (MQ=255)                   aCATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGa                                                        >  2:655971/1‑81 (MQ=255)                    cATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGa                                                        >  1:656002/1‑80 (MQ=255)                      tgatgTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATaa                                                     >  1:655973/1‑81 (MQ=255)                       gatgTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAat                                                    >  2:655970/1‑81 (MQ=255)                        atgTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAatc                                                   <  2:655915/81‑1 (MQ=255)                        atgTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAat                                                    >  2:655963/1‑80 (MQ=255)                           tCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAatcatc                                                <  1:655917/81‑1 (MQ=255)                             cccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCAtcgt                                              <  2:655917/81‑1 (MQ=255)                              ccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCAtcgtc                                             <  2:655918/81‑1 (MQ=255)                              ccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCAtcgtc                                             >  1:655977/1‑81 (MQ=255)                                ttACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGc                                           >  2:655968/1‑81 (MQ=255)                                  aCTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGcc                                          <  2:655919/80‑1 (MQ=255)                                   cttcttATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGcca                                         <  1:655921/80‑1 (MQ=255)                                       ttATGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAg                                     >  1:655975/1‑80 (MQ=255)                                         aTGACGAATACGCGGGTTCACTACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCaa                                  >  2:655977/1‑81 (MQ=255)                                         aTGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGc                                    <  1:655922/79‑1 (MQ=255)                                         aTGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCa                                   <  1:655923/80‑1 (MQ=255)                                          tGACGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAAt                                 <  2:655923/81‑1 (MQ=255)                                            aCGAATACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATAc                               <  1:655925/81‑1 (MQ=255)                                               aaTACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATAcgcc                            <  1:655926/81‑1 (MQ=255)                                                aTACGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATAcgccg                           <  1:655927/81‑1 (MQ=255)                                                  aCGCGGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACgccgcc                         >  2:655986/1‑81 (MQ=255)                                                   cgcgGGTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGccc                        <  2:655927/81‑1 (MQ=255)                                                      gggTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGc                     <  2:655928/81‑1 (MQ=255)                                                        gTTCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCAc                   <  1:655930/81‑1 (MQ=255)                                                         ttCACCACGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACt                  <  2:655930/81‑1 (MQ=255)                                                           caccacGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTAc                >  1:655988/1‑81 (MQ=255)                                                             ccacGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCg              >  2:655984/1‑81 (MQ=255)                                                              cacGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCTTCGTCGCCACCAGCAATACGCCGCCCTGCACTACCg              >  1:655992/1‑80 (MQ=255)                                                              cacGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCg              >  2:655988/1‑80 (MQ=255)                                                              cacGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCg              >  2:655996/1‑80 (MQ=255)                                                                cGCCGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGAt           <  2:655931/81‑1 (MQ=255)                                                                   cGTACAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGt        <  2:655932/81‑1 (MQ=255)                                                                      aCAGCAGTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTcgc     >  2:655997/1‑81 (MQ=255)                                                                         gcagTTCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCgcggc  >  1:655996/1‑81 (MQ=255)                                                                             |                                                                             TATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGC  >  NC_000913/3704740‑3704892

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