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%	155.7 / NA	37	coding (251/597 nt)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (22/15);  total (22/15)

GGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAG  >  NC_000913/4542208‑4542362                                                                                |                                                                            ggTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGac                                                                            >  2:941731/1‑80 (MQ=255) ggTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGac                                                                            <  2:941663/81‑2 (MQ=255)  gTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGacg                                                                           >  2:941714/1‑79 (MQ=255)     gAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGacgct                                                                         <  2:941664/80‑5 (MQ=255)     gAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGacgct                                                                         <  2:941665/80‑5 (MQ=255)     gAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTg                                                                       >  2:941721/1‑81 (MQ=255)        taaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAc                                                                    <  2:941666/81‑1 (MQ=255)         aaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAcc                                                                   >  1:941732/1‑81 (MQ=255)          aataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGAccc                                                                  <  1:941667/81‑1 (MQ=255)                 gCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACg                                                           >  2:941716/1‑81 (MQ=255)                 gCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACg                                                           >  1:941724/1‑81 (MQ=255)                    gACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGTAACGTGc                                                        >  1:941733/1‑81 (MQ=255)                     aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGc                                                        >  1:941711/1‑80 (MQ=255)                     aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGc                                                        >  1:941727/1‑80 (MQ=255)                       tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTaa                                                     >  2:941720/1‑81 (MQ=255)                         aagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACt                                                   >  2:941717/1‑81 (MQ=255)                          agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTg                                                  >  1:941755/1‑81 (MQ=255)                            aaCGGATTTTCTACCGTTCACCCGTTACGTTTTGGTGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGa                                                <  1:941668/81‑1 (MQ=255)                            aaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTgg                                                 >  2:941734/1‑80 (MQ=255)                             aCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGa                                                >  1:941739/1‑80 (MQ=255)                               ggATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAg                                             >  1:941715/1‑81 (MQ=255)                                gATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAgg                                            >  1:941722/1‑81 (MQ=255)                                 aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGgc                                           <  2:941669/81‑1 (MQ=255)                                  ttttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGgcg                                          >  1:941693/1‑81 (MQ=255)                                   tttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGgcgc                                         >  2:941726/1‑81 (MQ=255)                                          cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCgg                                  <  1:941670/81‑1 (MQ=255)                                               aCCCGTTACGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGa                             >  1:941740/1‑81 (MQ=255)                                                      aCGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtat                       <  2:941671/80‑1 (MQ=255)                                                       cGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtata                      <  1:941672/80‑1 (MQ=255)                                                       cGTTTTGATGAGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTatat                     >  1:941725/1‑81 (MQ=255)                                                             gatgaGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATatttattt               >  2:941723/1‑81 (MQ=255)                                                               tgaGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtct             <  2:941673/81‑1 (MQ=255)                                                                gaGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc            >  1:941708/1‑81 (MQ=255)                                                                 aGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc            <  2:941677/80‑1 (MQ=255)                                                                 aGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc            <  1:941676/80‑1 (MQ=255)                                                                 aGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc            <  2:941674/80‑1 (MQ=255)                                                                 aGCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcg           <  2:941675/81‑1 (MQ=255)                                                                  gCGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgc          <  1:941678/81‑1 (MQ=255)                                                                   cGTGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgcc         >  2:941736/1‑81 (MQ=255)                                                                     tGAAGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgccgc       <  1:941679/81‑1 (MQ=255)                                                                        aGCCGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCgg     >  2:941753/1‑80 (MQ=255)                                                                          ccGTGG‑ACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAg  <  1:941680/81‑1 (MQ=255)                                                                                |                                                                            GGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAG  >  NC_000913/4542208‑4542362

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