Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,609,008	T→C	P5P (CCA→CCG)	yneE ←	conserved inner membrane protein YneE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,609,008	0	T	C	100.0%	199.2 / NA	54	P5P (CCA→CCG)	yneE	conserved inner membrane protein YneE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/24);  total (30/24)

GATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTT  >  NC_000913/1608873‑1609144                                                                                                                                        |                                                                                                                                         gATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGAc                                                                                                                                      >  2:11783/1‑140 (MQ=255)   tGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGa                                                                                                                                    <  1:468531/140‑1 (MQ=255)    gCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGaa                                                                                                                                   <  1:320700/140‑1 (MQ=255)          aTGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCa                                                                                                                             <  1:168469/140‑1 (MQ=255)             cGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGATGCGGACGAACAATCATgg                                                                                                                          <  1:11783/140‑1 (MQ=255)                      gggCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAAt                                                                                                                               <  1:841817/126‑1 (MQ=255)                      gggCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAAt                                                                                                                               >  2:841817/1‑126 (MQ=255)                       ggCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAg                                                                                                                >  1:1057897/1‑140 (MQ=255)                       ggCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAg                                                                                                                >  1:1088679/1‑140 (MQ=255)                         cagcaTGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGc                                                                                                              >  1:655515/1‑140 (MQ=255)                            caTGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTg                                                                                                           >  1:868377/1‑140 (MQ=255)                            caTGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTg                                                                                                           >  1:646090/1‑140 (MQ=255)                                      aCAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGc                                                                                                 <  1:620275/140‑1 (MQ=255)                                       cAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGcc                                                                                                >  1:929885/1‑140 (MQ=255)                                          cGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTa                                                                                             >  2:658114/1‑140 (MQ=255)                                            aTAGAAAAGAGAAAATTGACTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTAtt                                                                                           >  2:585347/1‑140 (MQ=255)                                             tAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTa                                                                                          <  2:646090/140‑1 (MQ=255)                                                    gagaAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGa                                                                                   <  2:152944/140‑1 (MQ=255)                                                    gagaAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGa                                                                                   >  1:866246/1‑140 (MQ=255)                                                     agaAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGAt                                                                                  >  2:266667/1‑140 (MQ=255)                                                       aaaaTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGGGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTa                                                                                >  1:911674/1‑140 (MQ=255)                                                       aaaaTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTAtt                                                                                           >  1:775624/1‑129 (MQ=255)                                                       aaaaTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTAtt                                                                                           <  2:775624/129‑1 (MQ=255)                                                       aaaaTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTa                                                                                >  2:85534/1‑140 (MQ=255)                                                          aTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTa                                                                                <  1:219912/137‑1 (MQ=255)                                                          aTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTa                                                                                >  2:219912/1‑137 (MQ=255)                                                          aTTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGa                                                                             >  2:342900/1‑140 (MQ=255)                                                            tGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATc                                                                           <  2:866246/140‑1 (MQ=255)                                                            tGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATc                                                                           <  2:868377/140‑1 (MQ=255)                                                               gtagtaAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATcaca                                                                        >  1:814471/1‑140 (MQ=255)                                                                          gAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACg                                                             <  2:929885/140‑1 (MQ=255)                                                                          gAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACg                                                             >  1:106768/1‑140 (MQ=255)                                                                          gAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACg                                                             <  1:96222/140‑1 (MQ=255)                                                                             gATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTa                                                          >  2:394265/1‑140 (MQ=255)                                                                                   ttGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACa                                                    <  1:476778/140‑1 (MQ=255)                                                                                          ataCTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTAcc                                             <  1:85534/140‑1 (MQ=255)                                                                                           taCTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCt                                            >  1:258166/1‑140 (MQ=255)                                                                                                aGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAg                                       >  2:963921/1‑140 (MQ=255)                                                                                                aGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAg                                       >  2:33221/1‑140 (MQ=255)                                                                                                 gCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGc                                      >  2:93396/1‑140 (MQ=255)                                                                                                     tGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTac                                  >  2:1061275/1‑140 (MQ=255)                                                                                                     tGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTac                                  >  1:623492/1‑140 (MQ=255)                                                                                                     tGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTac                                  >  2:745315/1‑140 (MQ=255)                                                                                                      gCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTaca                                 <  1:311601/140‑1 (MQ=255)                                                                                                              aaaaTACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGc                         <  2:814471/140‑1 (MQ=255)                                                                                                               aaaTACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCa                        >  2:738919/1‑140 (MQ=255)                                                                                                                  tACGGCGCAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGgc                     <  2:258166/140‑1 (MQ=255)                                                                                                                        gcAGCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCa               <  1:562926/140‑1 (MQ=255)                                                                                                                           gCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgt            <  2:655515/140‑1 (MQ=255)                                                                                                                           gCCAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgt            >  2:777637/1‑140 (MQ=255)                                                                                                                            ccAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgtg           <  1:738919/140‑1 (MQ=255)                                                                                                                             cAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTgg          <  1:93396/140‑1 (MQ=255)                                                                                                                             cAGTGTTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTgg          <  1:342900/140‑1 (MQ=255)                                                                                                                               gtgtTGTTGCGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCa        >  2:1077720/1‑140 (MQ=255)                                                                                                                                 gttgttgcGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCAcc      <  2:1057897/140‑1 (MQ=255)                                                                                                                                    gttgcGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGt   >  1:951445/1‑140 (MQ=255)                                                                                                                                     ttgcGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGtt  <  2:623492/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTT  >  NC_000913/1608873‑1609144

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