Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,975,085	C→A	E80D (GAG→GAT)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,975,085	0	C	A	100.0%	181.0 / NA	52	E80D (GAG→GAT)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/24);  total (28/24)

GTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTC  >  NC_000913/1974949‑1975213                                                                                                                                         |                                                                                                                                 gTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCa                                                                                                                               <  2:191100/140‑1 (MQ=255)  ttGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCa                                                                                                                               >  2:256871/1‑139 (MQ=255)  ttGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCa                                                                                                                               <  1:256871/139‑1 (MQ=255)          aGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTg                                                                                                                      >  2:326326/1‑140 (MQ=255)               tGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCT‑‑TCTGGCTTCAt                                                                                                               >  2:1041314/1‑140 (MQ=255)                     atatAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAg                                                                                                           <  2:1034701/140‑1 (MQ=255)                     atatAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAg                                                                                                           >  2:862386/1‑140 (MQ=255)                       ataATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCa                                                                                                         <  1:822396/140‑1 (MQ=255)                       ataATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCa                                                                                                         >  2:364270/1‑140 (MQ=255)                         aaTCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAgg                                                                                                       <  2:431660/140‑1 (MQ=255)                           tCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGtt                                                                                                     >  1:1092766/1‑140 (MQ=255)                               aGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTc                                                                                                                 >  2:349854/1‑124 (MQ=255)                               aGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTc                                                                                                                 <  1:349854/124‑1 (MQ=255)                               aGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTcc                                                                                                 >  2:350129/1‑140 (MQ=255)                               aGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTcc                                                                                                 >  1:108661/1‑140 (MQ=255)                                gCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCGAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCa                                                                                                         >  1:130639/1‑131 (MQ=255)                                gCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCa                                                                                                         <  2:130639/131‑1 (MQ=255)                                 cTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCAt                                                                                               <  1:350129/140‑1 (MQ=255)                                   ggcggcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCa                                                                                             >  2:502521/1‑140 (MQ=255)                                         aTCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCg                                                                                       <  1:862386/140‑1 (MQ=255)                                           ccGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTg                                                                                     <  2:108661/140‑1 (MQ=255)                                                tctTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTc                                                                                <  1:326326/140‑1 (MQ=255)                                                           tttATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGc                                                                     >  2:85977/1‑140 (MQ=255)                                                                aaGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATAACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAg                                                                <  1:297855/140‑1 (MQ=255)                                                                  gCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGcc                                                              >  2:345235/1‑140 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAg                                                                >  1:993094/1‑135 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAg                                                                <  2:993094/135‑1 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCaaa                                                           >  1:728456/1‑140 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCaaa                                                           >  1:31318/1‑140 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCaaa                                                           >  1:151153/1‑140 (MQ=255)                                                                     tCGAGCTGTTCTTGCATGATGTCCTTCGTTCCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCaaa                                                           >  2:717887/1‑140 (MQ=255)                                                                           tGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTcc                                                     <  1:85977/140‑1 (MQ=255)                                                                                ttGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGccc                                                <  1:364270/140‑1 (MQ=255)                                                                                   cATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctc                                             <  2:1092766/140‑1 (MQ=255)                                                                                   cATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCcctc                                             <  2:988019/140‑1 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGccc                                                <  1:840696/133‑1 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGccc                                                >  2:840696/1‑133 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCttt                                         >  2:971704/1‑140 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCttt                                         >  1:1015427/1‑140 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCttt                                         >  1:1020966/1‑140 (MQ=255)                                                                                       atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCttt                                         >  1:1000454/1‑140 (MQ=255)                                                                                                  ttCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTggg                              <  2:151153/140‑1 (MQ=255)                                                                                                  ttCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTggg                              <  2:728456/140‑1 (MQ=255)                                                                                                   tCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGc                             >  1:535231/1‑140 (MQ=255)                                                                                                           cAGATTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGaa                     <  2:1000454/140‑1 (MQ=255)                                                                                                              aTTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGa                  <  1:502521/140‑1 (MQ=255)                                                                                                              aTTGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGa                  >  1:121275/1‑140 (MQ=255)                                                                                                                tGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATg                <  2:1015427/140‑1 (MQ=255)                                                                                                                tGATAATGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATg                <  1:717887/140‑1 (MQ=255)                                                                                                                      tGTCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTc          >  1:639787/1‑140 (MQ=255)                                                                                                                        tCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCaa        <  2:535231/140‑1 (MQ=255)                                                                                                                        tCGGTGTTGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCaa        >  1:1166361/1‑140 (MQ=255)                                                                                                                              ttGAGTTGCATATCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTc  >  2:1133357/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                 GTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTC  >  NC_000913/1974949‑1975213

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