Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,873,099	T→C	D80G (GAC→GGC)	ygbE ←	DUF3561 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,873,099	0	T	C	100.0%	198.7 / NA	55	D80G (GAC→GGC)	ygbE	DUF3561 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/25);  total (30/25)

AGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACC  >  NC_000913/2872951‑2873234                                                                                                                                                     |                                                                                                                                        aGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGccc                                                                                                                                       <  1:356467/151‑1 (MQ=255)   aaCATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCAt                                                                                                                                     <  2:141259/151‑1 (MQ=255)    aCATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATc                                                                                                                                    >  2:465307/1‑151 (MQ=255)            tcATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCaatgaatg                                                                                                                            >  2:570089/1‑151 (MQ=255)            tcATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCaatgaatg                                                                                                                            >  1:220325/1‑151 (MQ=255)              attattAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACCTAAACAGAGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATCCTGTAGCGTAGCTTTCCGCCCATCAATGAATGc                                                                                                                           >  1:118036/1‑150 (MQ=255)              attattAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGc                                                                                                                           >  2:169431/1‑150 (MQ=255)                tattaGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAg                                                                                                                         <  1:31238/150‑1 (MQ=255)                    aGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGCACAAATTACCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCa                                                                                                                     >  1:388100/1‑150 (MQ=255)                        aCGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATccc                                                                                                                >  2:491905/1‑151 (MQ=255)                           gTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCccac                                                                                                             >  2:333336/1‑151 (MQ=255)                               gAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCccacca                                                                                                           >  1:657540/1‑149 (MQ=255)                               gAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACt                                                                                                         >  2:499099/1‑151 (MQ=255)                                  ttATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTAcc                                                                                                      >  2:776818/1‑151 (MQ=255)                                      gAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACcgcg                                                                                                   <  1:341022/150‑1 (MQ=255)                                       aaCAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGAc                                                                                                 >  2:210274/1‑151 (MQ=255)                                         cAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACgg                                                                                               >  2:206904/1‑151 (MQ=255)                                                 aCGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAgc                                                                                       >  2:482333/1‑151 (MQ=255)                                                  cGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAgcg                                                                                      <  2:63809/151‑1 (MQ=255)                                                      aaaGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCag                                                                                  <  1:219133/151‑1 (MQ=255)                                                       aaGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCaga                                                                                 <  1:774259/151‑1 (MQ=255)                                                         gTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaa                                                                                >  2:459677/1‑150 (MQ=255)                                                            cAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaaaga                                                                             >  1:644512/1‑150 (MQ=255)                                                              aaTTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaaagaa                                                                            >  2:27655/1‑149 (MQ=255)                                                               aTTAGCCCAGTAGCCACATAAACAGTGCGCTAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAgg                                                                         >  1:118369/1‑151 (MQ=255)                                                                  aGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGcc                                                                       <  1:700177/150‑1 (MQ=255)                                                                     ccAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGgt                                                                   <  2:625975/151‑1 (MQ=255)                                                                             cacaTAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATa                                                            >  1:220387/1‑150 (MQ=255)                                                                                 tAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTAcc                                                                                                      >  1:579108/1‑104 (MQ=255)                                                                                 tAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTAcc                                                                                                      <  2:579165/104‑1 (MQ=255)                                                                                  aaaCAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaaga                                                        >  2:587047/1‑149 (MQ=255)                                                                                   aaCAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaaga                                                        >  1:199117/1‑148 (MQ=255)                                                                                     cAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAAc                                                   >  2:146000/1‑151 (MQ=255)                                                                                     cAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAAc                                                   <  2:384318/151‑1 (MQ=255)                                                                                      aGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACa                                                  <  1:27624/151‑1 (MQ=255)                                                                                         gcgcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGc                                                <  1:775696/150‑1 (MQ=255)                                                                                          cgcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCgt                                              <  2:70277/151‑1 (MQ=255)                                                                                           gcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCgtg                                             <  1:586179/151‑1 (MQ=255)                                                                                             cAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                           >  1:329319/1‑151 (MQ=255)                                                                                              aaaCATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTa                                          <  1:717973/151‑1 (MQ=255)                                                                                                  aTAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGc                                      <  1:11349/151‑1 (MQ=255)                                                                                                   tAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAAGGCAGCGAAATTCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                     >  1:523884/1‑151 (MQ=255)                                                                                                    aaTGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                     >  2:327903/1‑150 (MQ=255)                                                                                                    aaTGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                     >  2:585629/1‑150 (MQ=255)                                                                                                      tGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCAt                                   >  1:611838/1‑150 (MQ=255)                                                                                                       gCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATa                                  <  2:619584/150‑1 (MQ=255)                                                                                                       gCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATa                                  <  1:597329/150‑1 (MQ=255)                                                                                                        ccTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATaaa                                <  1:72863/151‑1 (MQ=255)                                                                                                                 aCCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCg                       <  2:649627/151‑1 (MQ=255)                                                                                                                    aGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCa                    <  1:333303/151‑1 (MQ=255)                                                                                                                       gTGAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCg                 <  1:642913/151‑1 (MQ=255)                                                                                                                         gAAAACAATACTGTAGCGTAGCTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGcc               >  1:49880/1‑151 (MQ=255)                                                                                                                                     gtagcgtagcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATaa    <  2:484407/150‑1 (MQ=255)                                                                                                                                     gtagcgtagcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATa     <  2:568058/149‑1 (MQ=255)                                                                                                                                     gtagcgtagcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAAc   <  2:738845/151‑1 (MQ=255)                                                                                                                                      tagcgtagcTTTCCGCCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGGGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAAcc  >  2:586729/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                        AGAACATAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACC  >  NC_000913/2872951‑2873234

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