Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,354,915	G→A	intergenic (‑198/+39)	lysU ← / ← dtpC	lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase/dipeptide/tripeptide:H(+) symporter DtpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,354,915	0	G	A	100.0%	295.6 / NA	72	intergenic (‑198/+39)	lysU/dtpC	lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase/dipeptide/tripeptide:H(+) symporter DtpC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/35);  total (37/35)

ATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCG  >  NC_000913/4354869‑4354960                                               |                                              aTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaaa                                            <  1:595349/50‑1 (MQ=255)  tGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAg                                           <  1:1729666/50‑1 (MQ=255)    gTTTAGATGACTTATTAATCTCGATTTGTTACTGTTATAAACTAAAAGt                                          <  1:3570129/49‑1 (MQ=38)    gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                         <  1:4212150/50‑1 (MQ=255)    gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                         <  1:2609168/50‑1 (MQ=255)     tttAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                         >  1:75711/1‑49 (MQ=255)      ttAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTaa                                       >  1:4789458/1‑50 (MQ=255)       tAAATTACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                      >  1:3095206/1‑50 (MQ=255)       tAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                      >  1:3237239/1‑50 (MQ=255)        aaaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                     >  1:3707972/1‑50 (MQ=255)        aaaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                      <  1:1395327/49‑1 (MQ=255)         aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                     <  1:2328960/49‑1 (MQ=255)         aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                     >  1:4255263/1‑49 (MQ=255)         aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                     >  1:4591951/1‑49 (MQ=255)          aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                   <  1:3395569/50‑1 (MQ=255)          aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                   <  1:5407536/50‑1 (MQ=255)          aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                   >  1:4248365/1‑50 (MQ=255)          aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                   <  1:4186609/50‑1 (MQ=255)           tGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAt                                  <  1:3129762/50‑1 (MQ=255)            gACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                 >  1:5580010/1‑50 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                 >  1:4365891/1‑49 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                 >  1:4941404/1‑49 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                 <  1:750302/49‑1 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                 >  1:1308382/1‑49 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTc                                <  1:3648326/50‑1 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTc                                <  1:583/50‑1 (MQ=255)             aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTc                                <  1:6181902/50‑1 (MQ=255)              cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                               >  1:4972072/1‑50 (MQ=255)              cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                               >  1:1544677/1‑50 (MQ=255)              cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                               >  1:1529567/1‑50 (MQ=255)                tattaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCAt                              >  1:5744327/1‑49 (MQ=255)                 attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                            >  1:2199786/1‑50 (MQ=255)                 attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATa                             <  1:1119574/49‑1 (MQ=255)                 attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATa                             <  1:911553/49‑1 (MQ=255)                  ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                           >  1:5126865/1‑50 (MQ=255)                  ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                           <  1:4300120/50‑1 (MQ=255)                    aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                          <  1:2999134/49‑1 (MQ=255)                    aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                          >  1:372125/1‑49 (MQ=255)                       ctcGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAtt                      >  1:2207138/1‑50 (MQ=255)                        tcGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTg                     >  1:3774426/1‑50 (MQ=255)                        tcGATTTGTTAGCGTTATAAACTAAAAGCTAATTTATTCATAACCAATTg                     >  1:4656943/1‑50 (MQ=255)                          gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAt                       <  1:1069123/46‑1 (MQ=255)                          gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTg                     <  1:1688348/48‑1 (MQ=255)                          gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                   <  1:5272662/50‑1 (MQ=255)                          gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                   >  1:681859/1‑50 (MQ=255)                           aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                  >  1:6037495/1‑50 (MQ=255)                           aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                  <  1:838066/50‑1 (MQ=255)                           aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                  >  1:2109367/1‑50 (MQ=255)                            tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                 >  1:382565/1‑50 (MQ=255)                            tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                 >  1:5835873/1‑50 (MQ=255)                            tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                 >  1:157930/1‑50 (MQ=255)                            tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                 >  1:680696/1‑50 (MQ=255)                            tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                 >  1:879607/1‑50 (MQ=255)                               gTTAGTGTTATAAGCTAAAAGTTAATTTATTCATAACCAATTGTTTTTa               <  1:30007/49‑1 (MQ=255)                               gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa               <  1:3234298/49‑1 (MQ=255)                               gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa               <  1:1050371/49‑1 (MQ=255)                               gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc              >  1:2793873/1‑50 (MQ=255)                                ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc              >  1:1353776/1‑49 (MQ=255)                                ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc              <  1:2691992/49‑1 (MQ=255)                                ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc              <  1:4654278/49‑1 (MQ=255)                                ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc              <  1:1227018/49‑1 (MQ=255)                                ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACg             <  1:2166335/50‑1 (MQ=255)                                  aGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa            <  1:3359421/49‑1 (MQ=255)                                   gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa            >  1:3326296/1‑48 (MQ=255)                                   gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGAc           <  1:2015470/49‑1 (MQ=255)                                    tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg         <  1:4583825/50‑1 (MQ=255)                                    tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg         <  1:935290/50‑1 (MQ=255)                                     gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt        >  1:2940898/1‑50 (MQ=255)                                     gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt        >  1:3725974/1‑50 (MQ=255)                                      ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt        >  1:3824504/1‑49 (MQ=255)                                          aaaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc   >  1:4838150/1‑50 (MQ=255)                                           aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc   <  1:728118/49‑1 (MQ=255)                                           aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCg  <  1:4255199/50‑1 (MQ=255)                                               |                                              ATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCG  >  NC_000913/4354869‑4354960

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