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

TTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACC  >  NC_000913/2012564‑2012831                                                                                                                                         |                                                                                                                                    ttATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTa                                                                                                                                   >  2:947726/1‑139 (MQ=255)   aTCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTaca                                                                                                                                 <  1:147416/139‑1 (MQ=255)      gcgTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCt                                                                                                                              <  2:43160/139‑1 (MQ=255)        gTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGc                                                                                                                            <  2:557078/139‑1 (MQ=255)                aTCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                    <  2:1308339/139‑1 (MQ=255)                aTCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATACCCCACGCGGTGAGATCTGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                    >  1:1035293/1‑139 (MQ=255)                  cGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcatca                                                                                                                  >  2:919908/1‑139 (MQ=255)                   gTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTacac                                                                                                                                <  2:596208/124‑1 (MQ=255)                   gTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTacac                                                                                                                                >  1:596208/1‑124 (MQ=255)                          gCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACGACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGa                                                                                                          >  1:1458067/1‑139 (MQ=255)                                     tCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                               <  1:838121/139‑1 (MQ=255)                                     tCAGGATCAGATATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAAGAGGTTACAAAACTTACACCTGCATGCCCATCACTTCCTGATACGCCGccac                                                                                               <  1:919908/139‑1 (MQ=255)                                       aGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacca                                                                                             <  1:1193961/139‑1 (MQ=255)                                                 aTCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                    <  2:39859/106‑1 (MQ=255)                                                 aTCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                                                    >  1:39859/1‑106 (MQ=255)                                                   ccATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCa                                                                                 >  1:764080/1‑139 (MQ=255)                                                    cATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCAc                                                                                >  1:1410334/1‑139 (MQ=255)                                                       tttAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTg                                                                             >  1:759180/1‑139 (MQ=255)                                                         tAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGaa                                                                           >  2:487879/1‑139 (MQ=255)                                                         tAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGaa                                                                           >  1:106575/1‑139 (MQ=255)                                                           aaCCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                                                         <  1:1304485/139‑1 (MQ=255)                                                                   aTTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTg                                                                 <  2:12348/139‑1 (MQ=255)                                                                                 ccccATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTaca                                                                                                                                 >  2:130940/1‑61 (MQ=255)                                                                                 ccccATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTaca                                                                                                                                 <  1:130940/61‑1 (MQ=255)                                                                                     atccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCttt                                               <  2:1035293/139‑1 (MQ=255)                                                                                      tccatGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCtttt                                              >  2:1484956/1‑139 (MQ=255)                                                                                          tGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCa                                          >  2:1184635/1‑139 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                               <  1:638368/85‑1 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccac                                                                                               >  2:638368/1‑85 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                       <  2:720189/125‑1 (MQ=255)                                                                                           gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                       >  1:720189/1‑125 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTACTTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                >  1:1381325/1‑78 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                <  2:1381325/78‑1 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                >  2:148977/1‑78 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                <  1:148977/78‑1 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                   >  2:1481068/1‑139 (MQ=255)                                                                                                 cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt                                   >  1:1236459/1‑139 (MQ=255)                                                                                                  gcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtc                                  <  1:1184635/139‑1 (MQ=255)                                                                                                        agaTCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATcaca                            <  1:1484956/139‑1 (MQ=255)                                                                                                             ggATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGtt                       <  2:294767/139‑1 (MQ=255)                                                                                                              gATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGttt                      >  2:416798/1‑139 (MQ=255)                                                                                                                tGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTaa                    >  2:782200/1‑139 (MQ=255)                                                                                                                         caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                                                     <  1:1460761/65‑1 (MQ=255)                                                                                                                         caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                                                     >  2:1460761/1‑65 (MQ=255)                                                                                                                         caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAacgccacgcc           >  1:1115858/1‑139 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAg                                                    >  2:1273723/1‑92 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAg                                                    <  1:1273723/92‑1 (MQ=255)                                                                                                                               gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTc     >  2:553538/1‑139 (MQ=255)                                                                                                                                 tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                       >  2:53225/1‑87 (MQ=255)                                                                                                                                 tACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                                                       <  1:53225/87‑1 (MQ=255)                                                                                                                                  aCAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                >  1:1245608/1‑45 (MQ=255)                                                                                                                                  aCAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                                                <  2:1245608/45‑1 (MQ=255)                                                                                                                                  aCAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCAcc  <  1:487879/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    TTATCGCGTATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACC  >  NC_000913/2012564‑2012831

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