Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,012,700	C→T	*105* (TAG→TAA)	fliE ←	flagellar basal‑body component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,012,700	0	C	T	100.0%	155.9 / NA	48	*105* (TAG→TAA)	fliE	flagellar basal‑body component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/23);  total (25/23)

GCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTC  >  NC_000913/2012569‑2012828                                                                                                                                    |                                                                                                                                 gcgTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTacacat                                                                                                                           >  1:28465/1‑137 (MQ=255)           aTCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                 <  2:562788/139‑1 (MQ=255)             cGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTCCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCACGCtcatca                                                                                                               >  1:674161/1‑139 (MQ=255)               tGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACt                                                                                                             >  1:1035395/1‑139 (MQ=255)                     gCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      <  2:446890/124‑1 (MQ=255)                     gCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      >  1:446890/1‑124 (MQ=255)                                 cAGGATCAGGTATCCATTTTAAACGAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacc                                                                                           <  1:634168/139‑1 (MQ=255)                                          gTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                                                  <  1:948585/139‑1 (MQ=255)                                              ccATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCa                                                                              >  2:1067098/1‑139 (MQ=255)                                              ccATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCa                                                                              >  2:1134531/1‑139 (MQ=255)                                                    tAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGaa                                                                        >  2:517580/1‑139 (MQ=255)                                                           gAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTATAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAt                                                                 <  1:324659/139‑1 (MQ=255)                                                                           gCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAg                                                 <  1:1134531/139‑1 (MQ=255)                                                                                atccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACg                                                                                <  1:910327/103‑1 (MQ=255)                                                                                atccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACg                                                                                >  2:910327/1‑103 (MQ=255)                                                                                          cACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCg                                  >  1:619663/1‑139 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                         <  2:69446/50‑1 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                         >  1:69446/1‑50 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATccc                                                                    >  1:558070/1‑103 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATccc                                                                    <  2:558070/103‑1 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                      <  2:79425/101‑1 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                      >  1:79425/1‑101 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAttt                                                               <  1:910750/108‑1 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAttt                                                               >  2:910750/1‑108 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGc                                               <  1:600815/124‑1 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGc                                               >  2:600815/1‑124 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                >  2:910772/1‑139 (MQ=255)                                                                                            cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                >  1:327937/1‑139 (MQ=255)                                                                                                 tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAtt                                                                <  2:360875/102‑1 (MQ=255)                                                                                                 tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAtt                                                                >  1:360875/1‑102 (MQ=255)                                                                                                   agaTCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATcaca                         <  2:327937/139‑1 (MQ=255)                                                                                                    gaTCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACAt                        <  1:517580/139‑1 (MQ=255)                                                                                                         gATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcatca                                                                                                               <  1:19825/47‑1 (MQ=255)                                                                                                         gATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcatca                                                                                                               >  2:19825/1‑47 (MQ=255)                                                                                                           tGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            >  2:1103531/1‑64 (MQ=255)                                                                                                           tGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            <  1:1103531/64‑1 (MQ=255)                                                                                                               aTTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacca                                                                                          >  2:108695/1‑62 (MQ=255)                                                                                                               aTTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacca                                                                                          <  1:108695/62‑1 (MQ=255)                                                                                                               aTTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAacgcc             <  2:603757/139‑1 (MQ=255)                                                                                                                    caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            >  2:454785/1‑55 (MQ=255)                                                                                                                    caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            <  1:454785/55‑1 (MQ=255)                                                                                                                     aacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCg       <  2:619663/139‑1 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                                                  >  2:187760/1‑59 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                                                  <  1:187760/59‑1 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtcat                             >  1:82193/1‑112 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtcat                             <  2:82193/112‑1 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTc  >  1:883651/1‑139 (MQ=255)                                                                                                                          gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTc  >  2:1052582/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                 GCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTC  >  NC_000913/2012569‑2012828

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