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%	391.7 / NA	95	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 (52/43);  total (52/43)

AATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGT  >  NC_000913/4354868‑4354961                                                |                                               aaTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaa                                              <  1:2777276/50‑1 (MQ=255) aaTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaa                                              <  1:1029946/50‑1 (MQ=255)  aTGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTaaaa                                             <  1:6251125/50‑1 (MQ=255)   tGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAg                                            <  1:5156648/50‑1 (MQ=255)    ggTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:2300653/1‑50 (MQ=255)    ggTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:7010613/1‑50 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAATTAAAAGt                                           >  1:1530448/1‑49 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                          <  1:4928921/50‑1 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                          <  1:1962454/50‑1 (MQ=255)     gTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGt                                           >  1:4158201/1‑49 (MQ=255)      tttAAATGACTTATTAATCTCTATTTGTTAGTGTTATAAACTAAAAGtt                                          >  1:5007094/1‑49 (MQ=255)      tttAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGtt                                          >  1:3889247/1‑49 (MQ=255)        tAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTaa                                        <  1:5074384/49‑1 (MQ=255)        tAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTaa                                        <  1:5302363/49‑1 (MQ=255)        tAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                       >  1:4280220/1‑50 (MQ=255)         aaaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                       <  1:3217372/49‑1 (MQ=255)         aaaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAt                                       <  1:2415271/49‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAttt                                     <  1:4780594/50‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAttt                                     <  1:123067/50‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                      >  1:6117119/1‑49 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                      <  1:3198969/49‑1 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                      >  1:369927/1‑49 (MQ=255)          aaTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAAtt                                      >  1:507590/1‑49 (MQ=255)           aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                    >  1:4382737/1‑50 (MQ=255)           aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                    >  1:5817093/1‑50 (MQ=255)           aTGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTa                                    <  1:410110/50‑1 (MQ=255)            tGACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAt                                   >  1:348996/1‑50 (MQ=255)             gACTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  >  1:308095/1‑50 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  >  1:5381798/1‑49 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  <  1:5029152/49‑1 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  <  1:2910488/49‑1 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  <  1:4841780/49‑1 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTAtt                                  >  1:6429889/1‑49 (MQ=255)              aCTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTc                                 <  1:4422959/50‑1 (MQ=255)               cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                                <  1:394940/50‑1 (MQ=255)               cTTATTAATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                                <  1:176075/50‑1 (MQ=255)                ttattaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCa                                >  1:1637528/1‑49 (MQ=255)                 tattaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCAt                               >  1:4942494/1‑49 (MQ=255)                  attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCAt                               <  1:689492/48‑1 (MQ=255)                  attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                             >  1:2103038/1‑50 (MQ=255)                  attaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATaa                             >  1:1458605/1‑50 (MQ=255)                   ttaATCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAc                            >  1:6259577/1‑50 (MQ=255)                     aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                           <  1:642332/49‑1 (MQ=255)                     aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                           >  1:1525022/1‑49 (MQ=255)                     aaTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAAcc                           <  1:7254940/49‑1 (MQ=255)                      aTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCaa                         >  1:3569442/1‑50 (MQ=255)                      aTCTCGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCa                          <  1:4037473/49‑1 (MQ=255)                        ctcGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAAtt                       >  1:2487430/1‑50 (MQ=255)                         tcGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTg                      >  1:4672439/1‑50 (MQ=255)                          cGATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGt                     >  1:63582/1‑50 (MQ=255)                           gATTTGTTAGTGTTATAAGCTAAAAGTTAATTTATTCATAACCAAtt                       <  1:5919424/47‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTg                      <  1:4407287/48‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                    <  1:6875723/50‑1 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                    >  1:6908711/1‑50 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtt                    >  1:4876236/1‑50 (MQ=255)                           gATTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGt                     <  1:7047708/49‑1 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   <  1:4493652/50‑1 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   >  1:4252346/1‑50 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   >  1:5700941/1‑50 (MQ=255)                            aTTTGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGttt                   >  1:1095218/1‑50 (MQ=255)                             tttGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGtttt                  >  1:686213/1‑50 (MQ=255)                               tGTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa                >  1:4386365/1‑50 (MQ=255)                                gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa                <  1:1334711/49‑1 (MQ=255)                                gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTa                <  1:487285/49‑1 (MQ=255)                                gTTAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTAc               >  1:1638300/1‑50 (MQ=255)                                 ttAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACg              <  1:4547400/50‑1 (MQ=255)                                  tAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:5002297/1‑50 (MQ=255)                                  tAGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:3854548/1‑50 (MQ=255)                                   aGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             <  1:2115928/49‑1 (MQ=255)                                   aGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:6228760/1‑49 (MQ=255)                                   aGTGTTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:1603676/1‑49 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGa             >  1:3406838/1‑48 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGAc            <  1:3202137/49‑1 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACa           <  1:4197461/50‑1 (MQ=255)                                    gtgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACa           <  1:3558563/50‑1 (MQ=255)                                     tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg          <  1:6633136/50‑1 (MQ=255)                                     tgtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAg          <  1:3259807/50‑1 (MQ=255)                                      gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         >  1:5011576/1‑50 (MQ=255)                                      gtTATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGt         >  1:3358829/1‑50 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        <  1:2551907/50‑1 (MQ=255)                                       ttATAAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        <  1:4575594/50‑1 (MQ=255)                                        tataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGtt        >  1:1360631/1‑49 (MQ=255)                                         ataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTaa      >  1:1578965/1‑50 (MQ=255)                                         ataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTaa      >  1:2542642/1‑50 (MQ=255)                                         ataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTaa      >  1:1643823/1‑50 (MQ=255)                                         ataAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTa       <  1:5850014/49‑1 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTaa      <  1:3538459/49‑1 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTaa      <  1:5471144/49‑1 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     >  1:7482824/1‑50 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     >  1:1642152/1‑50 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     <  1:6465971/50‑1 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     >  1:5587346/1‑50 (MQ=255)                                          taAACTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAAt     >  1:4878754/1‑50 (MQ=255)                                           aaaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc    >  1:5325255/1‑50 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc    <  1:5485811/49‑1 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATc    <  1:1318117/49‑1 (MQ=255)                                            aaCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCg   >  1:4888507/1‑50 (MQ=255)                                             aCTAAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGt  >  1:3854652/1‑50 (MQ=255)                                                |                                               AATGGTTTAAATGACTTATTAATCTCGATTTGTTAGTGTTATAAACTGAAAGTTAATTTATTCATAACCAATTGTTTTTACGACAGTTAATCGT  >  NC_000913/4354868‑4354961

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