Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,316,160	G→A	*891* (TAG→TAA)	rcsD →	phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,316,160	0	G	A	100.0%	136.1 / NA	42	*891* (TAG→TAA)	rcsD	phosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/20);  total (22/20)

TATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTA  >  NC_000913/2316027‑2316291                                                                                                                                      |                                                                                                                                    tATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                >  2:817929/1‑139 (MQ=255)        cATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTa                                                                                                                         >  1:361957/1‑139 (MQ=255)         aTGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTa                                                                                                                              <  2:363475/133‑1 (MQ=255)         aTGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTa                                                                                                                              >  1:363475/1‑133 (MQ=255)           gCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTa                                                                                                                              >  1:657837/1‑131 (MQ=255)           gCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTa                                                                                                                              <  2:657837/131‑1 (MQ=255)            cTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTAc                                                                                                                     <  2:951195/139‑1 (MQ=255)            cTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTAc                                                                                                                     <  1:817929/139‑1 (MQ=255)               aaTCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                  <  1:785581/139‑1 (MQ=255)                 tCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                <  2:361957/139‑1 (MQ=255)                 tCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                <  2:806379/139‑1 (MQ=255)                      tACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAata                                                                                                           >  2:418852/1‑139 (MQ=255)                           ggCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                      <  1:418852/139‑1 (MQ=255)                            gCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACg                                                                                                     >  1:458342/1‑139 (MQ=255)                               aGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTACAGGAATAGAAAAATGCATCAGCGACATTGACAGTTATGACAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTaa                                                                                                  >  1:708756/1‑139 (MQ=255)                                         tgAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGa                                                                                        >  1:473302/1‑139 (MQ=255)                                          gAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                       <  2:473302/139‑1 (MQ=255)                                          gAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                       <  2:940416/139‑1 (MQ=255)                                          gAAACGCTGGAACATCTGATTCGTGACAAGGATGTTCTAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                       <  2:708756/139‑1 (MQ=255)                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                               >  1:59485/1‑95 (MQ=255)                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                               <  2:59485/95‑1 (MQ=255)                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAcc                                                                                   <  1:599956/139‑1 (MQ=255)                                                 tGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATc                                                                                >  2:396740/1‑139 (MQ=255)                                                         cTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTc                                                                        <  2:458342/139‑1 (MQ=255)                                                                          tGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                       >  1:26573/1‑139 (MQ=255)                                                                            ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAAt                                                     >  1:200561/1‑139 (MQ=255)                                                                              ccAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCa                                                   >  2:427745/1‑139 (MQ=255)                                                                                    aTAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAg                                             >  1:300837/1‑139 (MQ=255)                                                                                             tACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    <  2:26573/139‑1 (MQ=255)                                                                                                       aCATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                >  2:484487/1‑37 (MQ=255)                                                                                                       aCATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                <  1:484487/37‑1 (MQ=255)                                                                                                                  tATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGaa               >  2:362333/1‑139 (MQ=255)                                                                                                                  tATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGaa               >  1:205518/1‑139 (MQ=255)                                                                                                                  tATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGaa               >  2:958942/1‑139 (MQ=255)                                                                                                                   aTGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAt              >  2:1021593/1‑139 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAc                                                                                    >  1:922872/1‑68 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAc                                                                                    <  2:922872/68‑1 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    <  2:75033/116‑1 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    >  1:75033/1‑116 (MQ=255)                                                                                                                         agagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTg                                  <  1:121210/113‑1 (MQ=255)                                                                                                                         agagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTg                                  >  2:121210/1‑113 (MQ=255)                                                                                                                               tGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTa  <  2:71946/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTA  >  NC_000913/2316027‑2316291

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