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

ATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGA  >  NC_000913/2012566‑2012833                                                                                                                                       |                                                                                                                                      aTCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTaca                                                                                                                                   >  2:181234/1‑139 (MQ=255)         taCGGATCGTGGATGTCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGTAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                           <  1:526735/139‑1 (MQ=255)         taCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                           <  2:318527/139‑1 (MQ=255)         taCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                           <  1:651482/139‑1 (MQ=255)         taCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                                                           <  2:561054/139‑1 (MQ=255)          actgATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGc                                                                                                                          <  1:491337/136‑1 (MQ=255)          aCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGc                                                                                                                          <  2:528061/139‑1 (MQ=255)                            aaaaCGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcg                                                                                                        <  1:291710/139‑1 (MQ=255)                                  ttCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                  <  1:125145/139‑1 (MQ=255)                                  ttCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                  <  2:90383/139‑1 (MQ=255)                                  ttCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                  <  2:619928/139‑1 (MQ=255)                                                         aaCCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                           <  2:460942/139‑1 (MQ=255)                                                           ccAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATccc                                                                         <  1:547141/139‑1 (MQ=255)                                                                         agtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                              >  1:264989/1‑72 (MQ=255)                                                                         agtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                              <  2:264989/72‑1 (MQ=255)                                                                          gtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                              >  2:615028/1‑71 (MQ=255)                                                                          gtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                              <  1:615028/71‑1 (MQ=255)                                                                           tttGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                         >  2:145247/1‑139 (MQ=255)                                                                                   atccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTa                                                                                                                                     <  2:601034/55‑1 (MQ=255)                                                                                   atccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTa                                                                                                                                     >  1:601034/1‑55 (MQ=255)                                                                                               cgcgGTGAGATCGGATGTAATTAACAACAGGTTACACAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                                 >  1:749542/1‑79 (MQ=255)                                                                                               cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                                 <  2:749542/79‑1 (MQ=255)                                                                                               cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                     >  1:429963/1‑139 (MQ=255)                                                                                               cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                     >  2:346118/1‑139 (MQ=255)                                                                                                              tGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                     <  2:37169/124‑1 (MQ=255)                                                                                                              tGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                     >  1:37169/1‑124 (MQ=255)                                                                                                                   ttAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATcaca                              <  2:295262/126‑1 (MQ=255)                                                                                                                   ttAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATcaca                              >  1:295262/1‑126 (MQ=255)                                                                                                                     aacaacAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCata                                          >  1:175061/1‑112 (MQ=255)                                                                                                                     aacaacAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCata                                          <  2:175061/112‑1 (MQ=255)                                                                                                                       caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcg                                                                                                        <  2:603869/48‑1 (MQ=255)                                                                                                                       caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcg                                                                                                        >  1:603869/1‑48 (MQ=255)                                                                                                                          caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                           <  2:436337/74‑1 (MQ=255)                                                                                                                          caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                           >  1:436337/1‑74 (MQ=255)                                                                                                                            ggTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCAGGtt        >  1:646276/1‑139 (MQ=255)                                                                                                                             gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCt                                                                                            <  1:631437/54‑1 (MQ=255)                                                                                                                             gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCt                                                                                            >  2:631437/1‑54 (MQ=255)                                                                                                                             gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTc       >  1:574713/1‑139 (MQ=255)                                                                                                                               tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAacg                    <  1:697111/124‑1 (MQ=255)                                                                                                                               tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAacg                    >  2:697111/1‑124 (MQ=255)                                                                                                                                aCAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCAcc    >  1:340345/1‑139 (MQ=255)                                                                                                                                  aaaaCTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCga  >  1:94296/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                      ATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGA  >  NC_000913/2012566‑2012833

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: ‑