Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	521,821	T→A	L469Q (CTG→CAG)	ybbP →	putative ABC transporter membrane subunit YbbP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	521,821	0	T	A	100.0%	171.6 / NA	49	L469Q (CTG→CAG)	ybbP	putative ABC transporter membrane subunit YbbP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (19/30);  total (19/30)

GATGCTGCTGAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGA  >  NC_000913/521684‑521955                                                                                                                                          |                                                                                                                                       gATGCTGCTGAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagc                                                                                                                                      >  2:1225921/1‑140 (MQ=255) gATGCTGCTGAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagc                                                                                                                                      >  2:326049/1‑140 (MQ=255) gATGCTGCTGAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagc                                                                                                                                      >  1:362685/1‑140 (MQ=255)      tgctgAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagctgctg                                                                                                                                 >  2:554414/1‑140 (MQ=255)      tgctgAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagctgctg                                                                                                                                 >  1:28829/1‑140 (MQ=255)      tgctgAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagctgctg                                                                                                                                 >  2:771349/1‑140 (MQ=255)      tgctgAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCAcagctgctg                                                                                                                                 <  1:326049/140‑1 (MQ=255)         tgAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGgtg                                                                                                                              <  2:405153/140‑1 (MQ=255)             tGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGc                                                                                                                          >  1:1209854/1‑140 (MQ=255)              gTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCg                                                                                                                         <  1:554414/140‑1 (MQ=255)                aCTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTg                                                                                                                       >  2:439434/1‑140 (MQ=255)                     gccgcATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGAt                                                                                                                  >  1:1004503/1‑140 (MQ=255)                     gccgcATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGAt                                                                                                                  >  2:800173/1‑140 (MQ=255)                       cgcATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCt                                                                                                                >  2:1013947/1‑140 (MQ=255)                          aTGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCt                                                                                                             <  2:76352/140‑1 (MQ=255)                            gACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCg                                                                                                           <  1:567299/140‑1 (MQ=255)                             aCGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGa                                                                                                          <  1:433263/140‑1 (MQ=255)                             aCGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGa                                                                                                          <  1:599343/140‑1 (MQ=255)                                   aaaTCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGTCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTg                                                                                                    >  2:839691/1‑140 (MQ=255)                                          tGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGgcagcag                                                                                             <  1:439434/140‑1 (MQ=255)                                            ccTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGgcagcagca                                                                                           <  2:271539/140‑1 (MQ=255)                                                 gcgcCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTAc                                                                                      >  2:999041/1‑140 (MQ=255)                                                        gCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGa                                                                               <  1:44089/140‑1 (MQ=255)                                                                  gccTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAac                                                                     <  2:784665/140‑1 (MQ=255)                                                                  gccTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACACCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAac                                                                     >  2:956404/1‑140 (MQ=255)                                                                      gTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAactact                                                                 <  2:408331/140‑1 (MQ=255)                                                                              aGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTttaatt                                                         <  2:370377/140‑1 (MQ=255)                                                                               gCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTttaatta                                                        >  2:703912/1‑140 (MQ=255)                                                                               gCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTttaatta                                                        >  2:519670/1‑140 (MQ=255)                                                                                ccGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTttaattaa                                                       >  1:888944/1‑140 (MQ=255)                                                                                   tGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACAt                                                    <  2:51389/140‑1 (MQ=255)                                                                                   tGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACAt                                                    <  2:362685/140‑1 (MQ=255)                                                                                        aaCGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCca                                               <  1:484402/140‑1 (MQ=255)                                                                                            ttAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGa                                           <  2:759823/140‑1 (MQ=255)                                                                                               aGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCca                                               <  2:618476/133‑1 (MQ=255)                                                                                               aGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCca                                               >  1:618476/1‑133 (MQ=255)                                                                                                    gCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAg                                   >  1:551360/1‑140 (MQ=255)                                                                                                     cTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGc                                  <  2:28829/140‑1 (MQ=255)                                                                                                          ggCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGc                             <  1:1015908/140‑1 (MQ=255)                                                                                                          ggCATTTTCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGc                             >  2:404579/1‑140 (MQ=255)                                                                                                               tttCGCTCTCCTTTATGCTGCTAGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTaaaa                        >  2:466910/1‑140 (MQ=255)                                                                                                                ttCGCTCTCCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAg                       <  1:469456/140‑1 (MQ=255)                                                                                                                       tcCTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGCAGCACCGCTAAAAGCGTTCCt                >  2:874181/1‑140 (MQ=255)                                                                                                                         cTTTATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcg              <  2:888944/140‑1 (MQ=255)                                                                                                                           ttATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcgcg            <  1:703912/140‑1 (MQ=255)                                                                                                                           ttATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcgcg            <  1:999041/140‑1 (MQ=255)                                                                                                                           ttATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcgcg            <  1:839691/140‑1 (MQ=255)                                                                                                                            tATGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCgg           <  2:551360/140‑1 (MQ=255)                                                                                                                             aTGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGa          <  1:827639/140‑1 (MQ=255)                                                                                                                             aTGCTGCTGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGa          <  2:1004503/140‑1 (MQ=255)                                                                                                                               gctgctGGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAAc        >  1:1083221/1‑140 (MQ=255)                                                                                                                                ctgctgGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACa       <  1:404579/140‑1 (MQ=255)                                                                                                                                 tgctgGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcgcg            >  1:327237/1‑134 (MQ=255)                                                                                                                                 tgctgGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTcgcg            <  2:327237/134‑1 (MQ=255)                                                                                                                                   ctgGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCa    >  2:1010335/1‑140 (MQ=255)                                                                                                                                   ctgGCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCa    >  2:1219487/1‑140 (MQ=255)                                                                                                                                     ggCACAGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGa  <  1:322135/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       GATGCTGCTGAATGTACTTCGCCGCATGACGCTGAAATCGCTGCCTCTGCGCCTGGCGGTTAGCCGCCTGTTACGTCAGCCGTGGTCAACGTTAAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGA  >  NC_000913/521684‑521955

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: ‑