Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,288	Δ1 bp	coding (252/597 nt)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,287	0	A	.	100.0%	133.7 / NA	32	coding (251/597 nt)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (18/14);  total (18/14)

GGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGT  >  NC_000913/4542208‑4542363                                                                                |                                                                             ggTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGac                                                                             <  2:871026/81‑2 (MQ=255)  gTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGacg                                                                            >  2:871062/1‑79 (MQ=255)      aataaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTgg                                                                       <  2:871027/81‑1 (MQ=255)          aataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAccc                                                                   <  2:871029/81‑1 (MQ=255)          aataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAcc                                                                    >  1:871072/1‑80 (MQ=255)          aataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAcc                                                                    <  1:871029/80‑1 (MQ=255)               tCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGaa                                                              >  1:871081/1‑81 (MQ=255)                  ccGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGt                                                           >  2:871060/1‑81 (MQ=255)                   cGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTg                                                          >  1:871107/1‑81 (MQ=255)                    gACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGc                                                         <  1:871031/81‑1 (MQ=255)                     aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCt                                                        <  1:871032/81‑1 (MQ=255)                     aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCt                                                        <  2:871032/81‑1 (MQ=255)                       tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTaa                                                      >  1:871095/1‑81 (MQ=255)                        gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAAc                                                     >  1:871069/1‑81 (MQ=255)                        gaagaaCGGATTTTCTACCGCTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAAc                                                     >  1:871074/1‑81 (MQ=255)                         aagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACt                                                    <  1:871034/81‑1 (MQ=255)                                gATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAgg                                             >  1:871078/1‑81 (MQ=255)                                 aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGgc                                            >  1:871073/1‑81 (MQ=255)                                 aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGgc                                            <  2:871034/81‑1 (MQ=255)                                        aCCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc                                     >  1:871084/1‑81 (MQ=255)                                         ccGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg                                    >  2:871081/1‑81 (MQ=255)                                           gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                  >  1:871085/1‑81 (MQ=255)                                             tCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACt                                >  1:871086/1‑81 (MQ=255)                                              cACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACt                                <  1:871036/80‑1 (MQ=255)                                                 ccGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACg                            <  2:871036/81‑1 (MQ=255)                                                        gTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatt                     >  1:871068/1‑81 (MQ=255)                                                            tgatgaGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatttatt                 >  2:871088/1‑81 (MQ=255)                                                              atgaGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtc               <  2:871037/81‑1 (MQ=255)                                                                 aGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc             <  2:871038/80‑1 (MQ=255)                                                                   cGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgcc          >  1:871076/1‑81 (MQ=255)                                                                    gTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgccg         <  2:871039/81‑1 (MQ=255)                                                                       aaGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCgg      <  1:871041/81‑1 (MQ=255)                                                                        aGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCgg      >  2:871079/1‑80 (MQ=255)                                                                           cGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGt  >  1:871099/1‑81 (MQ=255)                                                                                |                                                                             GGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGT  >  NC_000913/4542208‑4542363

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