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%	133.8 / NA	42	*105* (TAG→TAA)	fliE	flagellar basal‑body component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/21);  total (21/21)

TTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTC  >  NC_000913/2012564‑2012828                                                                                                                                         |                                                                                                                                 ttATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTa                                                                                                                                <  2:1185671/139‑1 (MQ=255)   aTCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTaca                                                                                                                              >  1:1048774/1‑139 (MQ=255)         tataCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                         >  2:42164/1‑138 (MQ=255)         tataCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                         <  1:42164/138‑1 (MQ=255)                aTCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                 <  1:680676/139‑1 (MQ=255)                   gTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGc                                                                                                                     >  2:1238157/1‑132 (MQ=255)                   gTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGc                                                                                                                     <  1:1238157/132‑1 (MQ=255)                   gTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCAc                                                                                                              >  1:929247/1‑139 (MQ=255)                    tGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACt                                                                                                             <  2:1048774/139‑1 (MQ=255)                                  cGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcgccgc                                                                                               <  2:929247/139‑1 (MQ=255)                                    ttCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAAATTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                             >  2:1040195/1‑139 (MQ=255)                                          aTCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGTTACGCCGCCACCAGCt                                                                                       <  2:1205153/139‑1 (MQ=255)                                                       tttAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      >  2:677273/1‑95 (MQ=255)                                                       tttAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      <  1:677273/95‑1 (MQ=255)                                                         tAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGaa                                                                        >  2:1214139/1‑139 (MQ=255)                                                             ccAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATccc                                                                    <  1:1040195/139‑1 (MQ=255)                                                                 aaaTTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAtt                                                                >  1:573553/1‑139 (MQ=255)                                                                        gtgagtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAg                                                         >  2:4277/1‑139 (MQ=255)                                                                        gtgagtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAg                                                         <  2:406098/139‑1 (MQ=255)                                                                             tttGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                    <  1:426145/139‑1 (MQ=255)                                                                                gCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAg                                                 <  1:1407182/139‑1 (MQ=255)                                                                                      tccatGCCCGACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                 >  1:1029949/1‑69 (MQ=39)                                                                                      tccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACGCCTGCATGCtcat                                                                                                                 <  2:1029949/69‑1 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCt                                                                                                                           <  1:805583/54‑1 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCt                                                                                                                           >  2:805583/1‑54 (MQ=255)                                                                                              ccACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATc                                   >  1:682239/1‑139 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      >  1:911406/1‑53 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                      <  2:911406/53‑1 (MQ=255)                                                                                                                  tAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCGACCAGCTTATTACGCACCTGAATCCGCATTTGCATAGaaa                                                      <  1:935565/100‑1 (MQ=255)                                                                                                                  tAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGaaa                                                      >  2:935565/1‑100 (MQ=255)                                                                                                                         caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACg                                                                                <  2:1419412/67‑1 (MQ=255)                                                                                                                         caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACg                                                                                >  1:1419412/1‑67 (MQ=255)                                                                                                                            caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCatat                                    <  1:62360/108‑1 (MQ=255)                                                                                                                            caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCatat                                    >  2:62360/1‑108 (MQ=255)                                                                                                                            caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCggg     >  1:4672/1‑139 (MQ=255)                                                                                                                              ggTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacca                                                                                          <  2:553506/52‑1 (MQ=255)                                                                                                                              ggTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacca                                                                                          >  1:553506/1‑52 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCtttt                                           >  2:1327776/1‑98 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCtttt                                           <  1:1327776/98‑1 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTc  >  1:913318/1‑139 (MQ=255)                                                                                                                                 tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            >  1:819186/1‑47 (MQ=255)                                                                                                                                 tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                            <  2:819186/47‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 TTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTC  >  NC_000913/2012564‑2012828

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