Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	522,143	A→G	E576E (GAA→GAG)	ybbP →	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	522,143	0	A	G	100.0%	170.5 / NA	48	E576E (GAA→GAG)	ybbP	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

GCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAAGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGCCTAATTT  >  NC_000913/522014‑522271                                                                                                                                  |                                                                                                                                 gCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTg                                                                                                                         >  1:702586/1‑139 (MQ=255)    gACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCg                                                                                                                      <  2:1272588/139‑1 (MQ=255)      cAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGaa                                                                                                                    <  2:669993/139‑1 (MQ=255)                   aGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGtt                                                                                                       >  2:1084540/1‑139 (MQ=255)                    gATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATggagga                                                                                                                                >  2:592134/1‑113 (MQ=255)                    gATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATggagga                                                                                                                                >  2:1127882/1‑113 (MQ=255)                    gATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGtcg                                                                                                      <  2:702586/139‑3 (MQ=255)                      tGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATCATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  1:902440/139‑1 (MQ=255)                      tGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  1:1427419/139‑1 (MQ=255)                       gAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGccc                                                                                                   >  2:1454897/1‑139 (MQ=255)                       gAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGccc                                                                                                   >  2:1341987/1‑139 (MQ=255)                                     cgAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCg                                                                                                                      >  1:1140012/1‑106 (MQ=255)                                     cgAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCg                                                                                                                      <  2:1140012/106‑1 (MQ=255)                                     cgAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTg                                                                                     <  1:1454897/139‑1 (MQ=255)                                       aaCTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    >  1:130239/1‑122 (MQ=255)                                       aaCTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  2:130239/122‑1 (MQ=255)                                            aaTCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTa                                                                              >  1:915226/1‑139 (MQ=255)                                                     tGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATa                                                                     <  1:1341987/139‑1 (MQ=255)                                                      ggCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    >  2:163677/1‑107 (MQ=255)                                                      ggCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  1:163677/107‑1 (MQ=255)                                                        cAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATAccc                                                                  >  2:1522331/1‑139 (MQ=255)                                                            aTACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGGGAAACGCTTAAACGt                                                                                                        <  2:1129674/97‑1 (MQ=255)                                                            aTACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGt                                                                                                        >  1:1129674/1‑97 (MQ=255)                                                                              aTCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGaa                                                                                                                    >  2:1183455/1‑67 (MQ=255)                                                                              aTCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGaa                                                                                                                    <  1:1183455/67‑1 (MQ=255)                                                                                  gATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGc                                        >  2:1304960/1‑139 (MQ=255)                                                                                   aTTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACg                                                                                                         >  1:859596/1‑73 (MQ=255)                                                                                   aTTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACg                                                                                                         <  2:859596/73‑1 (MQ=255)                                                                                   aTTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCt                                                                                                  <  2:1269325/80‑1 (MQ=255)                                                                                   aTTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCt                                                                                                  >  1:1269325/1‑80 (MQ=255)                                                                                   aTTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGcc                                       >  2:608673/1‑139 (MQ=255)                                                                                    ttGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTa                                                                              >  2:209285/1‑99 (MQ=255)                                                                                    ttGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTa                                                                              <  1:209285/99‑1 (MQ=255)                                                                                    ttGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTaa                                                  >  2:711624/1‑127 (MQ=255)                                                                                    ttGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTaa                                                  <  1:711624/127‑1 (MQ=255)                                                                                       tCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTgcg                                   >  2:890664/1‑139 (MQ=255)                                                                                        cGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTgcgc                                  >  2:70132/1‑139 (MQ=255)                                                                                          ccGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  1:1428342/71‑1 (MQ=255)                                                                                          ccGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    >  2:1428342/1‑71 (MQ=255)                                                                                           cGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    <  2:1310054/70‑1 (MQ=255)                                                                                           cGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGcc                                                                                                    >  1:1310054/1‑70 (MQ=255)                                                                                           cGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCaaa                               <  2:915226/139‑1 (MQ=255)                                                                                                aCTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGa                          <  1:154321/139‑1 (MQ=255)                                                                                                     ccgccAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTggg                     >  2:1444324/1‑139 (MQ=255)                                                                                                     ccgccAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTggg                     <  1:1522331/139‑1 (MQ=255)                                                                                                        ccAAAAGCCGATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGaaa                  >  2:1402386/1‑139 (MQ=255)                                                                                                                 gATGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGc         <  1:608673/139‑1 (MQ=255)                                                                                                                   tGAAGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGCCt       >  2:1388693/1‑139 (MQ=255)                                                                                                                      aGTGTCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGCCTAAt    <  1:116821/139‑1 (MQ=255)                                                                                                                        tgtCGATGGAGGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGCCTAAttt  <  1:890664/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                 GCCGACAGAAGGTAATGAAGATGAGGCGCTTAACCGCGAACTCAATCTTACCTGGCAAAATACGCGGCCCGATCATAATCCGATTGTCGCCGGTAACTGGCCGCCAAAAGCCGATGAAGTGTCGATGGAAGAGGGGCTGGCGAAACGCTTAAACGTTGCCCTCGGTGATACCGTGACTTTTATGGGCGATACCCAGGAGTTCCGCGCTAAAGTGACCAGCCTGCGCAAAGTGGACTGGGAAAGTCTGCGGCCTAATTT  >  NC_000913/522014‑522271

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