Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	521,921	A→G	V502V (GTA→GTG)	ybbP →	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	521,921	0	A	G	100.0%	157.9 / NA	41	V502V (GTA→GTG)	ybbP	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/27);  total (14/27)

AAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGG  >  NC_000913/521777‑522068                                                                                                                                                 |                                                                                                                                                    aaGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCg                                                                                                                                               <  2:659012/151‑1 (MQ=255)   gTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCt                                                                                                                                             <  1:437422/151‑1 (MQ=255)      aGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTaa                                                                                                                                           >  2:212395/1‑150 (MQ=255)        cTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTaaa                                                                                                                                          >  1:160813/1‑149 (MQ=255)                   ttCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTcgcg                                                                                                                             >  1:401110/1‑151 (MQ=255)                     cGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCgg                                                                                                                            >  2:551528/1‑150 (MQ=255)                     cGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGa                                                                                                                           <  1:494537/151‑1 (MQ=255)                              ttATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATc                                                                                                                      >  2:596220/1‑147 (MQ=255)                              ttATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAg                                                                                                                    <  1:102705/149‑1 (MQ=255)                                aTGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGAt                                                                                                                  >  1:670028/1‑149 (MQ=255)                                  gctgctGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTc                                                                                                              <  1:486487/151‑1 (MQ=255)                                     gctgGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCg                                                                                                           >  2:65227/1‑151 (MQ=255)                                                 gctgGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTa                                                                                                <  2:567912/150‑1 (MQ=255)                                                                      gCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTg                                                                          <  2:551045/151‑1 (MQ=255)                                                                          gACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACgg                                                                      <  2:367653/151‑1 (MQ=255)                                                                           aCCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGc                                                                     <  1:17163/151‑1 (MQ=255)                                                                            ccGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCg                                                                    >  2:61589/1‑151 (MQ=255)                                                                                 ggCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGAtt                                                                 <  1:216249/149‑1 (MQ=255)                                                                                           gcTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCg                                                     >  2:640639/1‑151 (MQ=255)                                                                                               ccTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGAc                                                   >  2:279810/1‑149 (MQ=255)                                                                                                  ccAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAg                                              <  1:142281/151‑1 (MQ=255)                                                                                                          cccGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGaa                                      <  1:389719/151‑1 (MQ=255)                                                                                                            cGAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGa                                    <  1:476760/151‑1 (MQ=255)                                                                                                             gAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGAt                                   <  1:71207/151‑1 (MQ=255)                                                                                                             gAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGAt                                   <  1:82981/151‑1 (MQ=255)                                                                                                             gAACTACTTTTTAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGAt                                   <  1:281620/151‑1 (MQ=255)                                                                                                                     ttttAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGgcgc                            <  2:416877/150‑1 (MQ=255)                                                                                                                     ttttAATTAACATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCt                           >  1:669801/1‑151 (MQ=255)                                                                                                                        taattaaCATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTa                         <  1:568822/150‑1 (MQ=255)                                                                                                                         aattaaCATCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTaa                        <  2:321925/150‑1 (MQ=255)                                                                                                                                aTCGCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAAc                <  1:212394/151‑1 (MQ=255)                                                                                                                                   gCCACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCa             <  1:61590/151‑1 (MQ=255)                                                                                                                                    ccACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa            <  2:165517/151‑1 (MQ=255)                                                                                                                                    ccACAGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCaa            <  1:493697/151‑1 (MQ=255)                                                                                                                                     cacaGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAAt           >  1:560757/1‑151 (MQ=255)                                                                                                                                      acaGAACAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATc          <  2:261640/151‑1 (MQ=255)                                                                                                                                           aCAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTAcc     <  1:65227/151‑1 (MQ=255)                                                                                                                                            cAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCt    <  1:559409/151‑1 (MQ=255)                                                                                                                                            cAGGTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCt    <  2:11352/151‑1 (MQ=255)                                                                                                                                              ggTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTgg  >  1:622543/1‑151 (MQ=255)                                                                                                                                              ggTGGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGAAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTgg  >  1:601739/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    AAGTCAGCTTTCGGCATTTTCGCTCTCCTTTATGCTGCTGGCACTGCTGCTGGTGTTGCGTGGCGATCTGCTCGACCGCTGGCAGCAGCAGCTACCTCCAGAAAGCCCGAACTACTTTTTAATTAACATCGCCACAGAACAGGTAGCACCGCTAAAAGCGTTCCTCGCGGAACATCAGATAGTCCCGGAATCGTTTTATCCGGTGGTGCGGGCGCGGCTGACGGCGATTAACGATAAGCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGG  >  NC_000913/521777‑522068

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