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%	261.5 / NA	65	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 (32/33);  total (32/33)

AATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGT  >  NC_000913/4354868‑4354961                                                |                                               aaTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaa                                              <  1:4370244/50‑1 (MQ=255)  aTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaaa                                             <  1:935865/50‑1 (MQ=255)    ggTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:3067991/1‑50 (MQ=255)    ggTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:2856467/1‑50 (MQ=255)    ggTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:627677/1‑50 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                          <  1:2676188/50‑1 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:2794485/1‑49 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           <  1:2657439/49‑1 (MQ=255)      tttAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                          >  1:3750237/1‑49 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAttt                                     <  1:1567470/50‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAttt                                     <  1:476116/50‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                      <  1:3989407/49‑1 (MQ=255)           aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                    <  1:2079715/50‑1 (MQ=255)            tGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAt                                   <  1:362476/50‑1 (MQ=255)             gACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  >  1:536411/1‑50 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTc                                 <  1:689492/50‑1 (MQ=255)               cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                                >  1:481085/1‑50 (MQ=255)                  attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                             >  1:2945317/1‑50 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                             <  1:2102019/49‑1 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                             <  1:1634386/49‑1 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                            >  1:3292924/1‑50 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                            >  1:3874045/1‑50 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                            >  1:3837207/1‑50 (MQ=255)                    taATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                            >  1:1945985/1‑49 (MQ=255)                     aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                           >  1:1514305/1‑49 (MQ=255)                      aTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCaa                         <  1:843436/50‑1 (MQ=255)                      aTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCa                          <  1:3663531/49‑1 (MQ=255)                      aTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCa                          <  1:3802608/49‑1 (MQ=255)                       tctcGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAt                        >  1:4215191/1‑50 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAtt                       <  1:3533826/47‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAt                        <  1:3356663/46‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTg                      <  1:834455/48‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                    >  1:280886/1‑50 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                    <  1:3732552/50‑1 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   >  1:4267108/1‑50 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   <  1:2691032/50‑1 (MQ=255)                             tttGTTAGTGTTATAAGCTAAAAGTTAATTTATTCATAACCAATTGtttt                  <  1:3795682/50‑1 (MQ=255)                             tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                  >  1:2437151/1‑50 (MQ=255)                             tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                  >  1:2234247/1‑50 (MQ=255)                             tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                  >  1:395940/1‑50 (MQ=255)                               tGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttttt                 <  1:3014366/49‑1 (MQ=255)                               tGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa                >  1:118349/1‑50 (MQ=255)                                gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa                <  1:2459563/49‑1 (MQ=255)                                 ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc               <  1:1763477/49‑1 (MQ=255)                                 ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc               <  1:1852243/49‑1 (MQ=255)                                 ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACg              <  1:4342056/50‑1 (MQ=255)                                  tAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:1419840/1‑50 (MQ=255)                                   aGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:2059756/1‑49 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:1001593/1‑48 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:2535323/1‑48 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:2181432/1‑48 (MQ=255)                                     tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg          <  1:2913600/50‑1 (MQ=255)                                     tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg          >  1:1924182/1‑50 (MQ=255)                                      gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         >  1:3233909/1‑50 (MQ=255)                                      gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         >  1:4203947/1‑50 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        <  1:2352245/50‑1 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        <  1:2956184/50‑1 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         <  1:4328082/49‑1 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         >  1:1999401/1‑49 (MQ=255)                                        tataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        >  1:4092058/1‑49 (MQ=255)                                        tataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        >  1:18368/1‑49 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     >  1:1279195/1‑50 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     <  1:3551018/50‑1 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc    <  1:213727/49‑1 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCg   <  1:2120853/50‑1 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCg   <  1:3945901/50‑1 (MQ=255)                                             aCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGt  <  1:2693525/50‑1 (MQ=255)                                                |                                               AATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGT  >  NC_000913/4354868‑4354961

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