Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,541,580	A→G	intergenic (+21/‑457)	fimB → / → fimE	regulator for fimA/regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,541,580	0	A	G	100.0%	80.6 / NA	24	intergenic (+21/‑457)	fimB/fimE	regulator for fimA/regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/12);  total (12/12)

TCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACA  >  NC_000913/4541448‑4541710                                                                                                                                     |                                                                                                                                   tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGtt                                                                                                                             >  1:100609/1‑140 (MQ=255) tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGtt                                                                                                                             >  2:149535/1‑140 (MQ=255) tCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGtt                                                                                                                             >  1:242097/1‑140 (MQ=255)           cGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACa                                                                                                                   <  1:149535/140‑1 (MQ=255)             cAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACat                                                                                                                  >  2:216589/1‑139 (MQ=255)             cAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACat                                                                                                                  <  1:216589/139‑1 (MQ=255)                    cGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTa                                                                                                          <  1:263364/140‑1 (MQ=255)                    cGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTa                                                                                                          <  2:164228/140‑1 (MQ=255)                               cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCaa                                                                                               <  2:100609/140‑1 (MQ=255)                                  gTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCaaaaa                                                                                            <  2:242097/140‑1 (MQ=255)                                        cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGATCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTa                                                                                      >  1:263744/1‑140 (MQ=255)                                        cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTa                                                                                      >  1:248864/1‑140 (MQ=255)                                        cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTa                                                                                      >  2:32384/1‑140 (MQ=255)                                        cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTa                                                                                      >  2:198623/1‑140 (MQ=255)                                            aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCa                                                                                  >  1:276990/1‑140 (MQ=255)                                                     gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                                         <  2:276990/140‑1 (MQ=255)                                                         gTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                                                     <  1:40653/140‑1 (MQ=255)                                                                               agCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTa                                               >  1:125958/1‑140 (MQ=255)                                                                                   agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                                           <  2:263744/140‑1 (MQ=255)                                                                                               cGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACa                               >  1:174855/1‑140 (MQ=255)                                                                                                               agagATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                          >  1:342047/1‑129 (MQ=255)                                                                                                               agagATAGTTGAGATACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                          <  2:342047/129‑1 (MQ=255)                                                                                                                            aTACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACa  <  2:125958/140‑1 (MQ=255)                                                                                                                            aTACCAGGGGTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACa  <  2:345826/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   TCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACA  >  NC_000913/4541448‑4541710

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