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%	192.4 / NA	59	*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 (23/36);  total (23/36)

GTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAC  >  NC_000913/2316026‑2316292                                                                                                                                       |                                                                                                                                     gTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                                  <  2:740720/139‑1 (MQ=255)  tATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                 <  2:274651/139‑1 (MQ=255)  tATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAAAAAgg                                                                                                                                 >  2:793437/1‑139 (MQ=255)   aTTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                                <  2:859173/139‑1 (MQ=255)   aTTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                                <  1:501301/139‑1 (MQ=255)   aTTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                                <  1:965798/139‑1 (MQ=255)    tttGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTa                                                                                                                               <  2:221741/139‑1 (MQ=255)      tGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGc                                                                                                                             <  1:503624/139‑1 (MQ=255)        ccATGCTAAATCTGGTTCCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           <  2:835389/139‑1 (MQ=255)        ccATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           <  1:524950/139‑1 (MQ=255)              tAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACa                                                                                                                     <  2:140942/139‑1 (MQ=255)              tAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACa                                                                                                                     <  1:85984/139‑1 (MQ=255)                aaTCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                   <  1:591083/139‑1 (MQ=255)                aaTCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                   <  1:204916/139‑1 (MQ=255)                  tCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                 <  1:317003/139‑1 (MQ=255)                  tCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                 <  1:793437/139‑1 (MQ=255)                   cTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAAc                                                                                                                <  2:709191/139‑1 (MQ=255)                      gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                 >  2:897836/1‑119 (MQ=255)                      gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                 <  1:897836/119‑1 (MQ=255)                       tACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAata                                                                                                            >  1:233046/1‑139 (MQ=255)                            ggCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                 <  1:220148/129‑1 (MQ=255)                            ggCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                 >  2:220148/1‑129 (MQ=255)                                           gAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                        <  2:507721/139‑1 (MQ=255)                                            aaaCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATg                                                                                       >  1:989280/1‑139 (MQ=255)                                             aaCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGa                                                                                      >  2:517670/1‑139 (MQ=255)                                               cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAcc                                                                                    >  1:385496/1‑139 (MQ=255)                                                     aaCATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGa                                                                              <  1:319746/139‑1 (MQ=255)                                                              ttCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGt                                                                     >  1:110767/1‑139 (MQ=255)                                                                 gTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCg                                                                  <  2:233046/139‑1 (MQ=255)                                                                  tGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           <  1:603456/81‑1 (MQ=255)                                                                  tGAGAAGGATGTTCCAGCAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           >  2:603456/1‑81 (MQ=255)                                                                           tGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                        >  1:1010527/1‑139 (MQ=255)                                                                             ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAAt                                                      <  1:517670/139‑1 (MQ=255)                                                                                  ggAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACtt                                                 >  2:717899/1‑139 (MQ=255)                                                                                   gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           >  2:698894/1‑64 (MQ=255)                                                                                   gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                           <  1:698894/64‑1 (MQ=255)                                                                                   gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCgg                                                                 >  2:1043953/1‑122 (MQ=255)                                                                                   gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCgg                                                                 <  1:1043953/122‑1 (MQ=255)                                                                                    aaTAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGGCCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGa                                               >  2:1041240/1‑139 (MQ=255)                                                                                        gAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaa                                           >  1:48279/1‑139 (MQ=255)                                                                                         aaaaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAc                                                                                     <  1:701068/96‑1 (MQ=255)                                                                                         aaaaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAc                                                                                     >  2:701068/1‑96 (MQ=255)                                                                                             aTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGGCCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGa                                      <  1:1041240/139‑1 (MQ=255)                                                                                                                 gTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCg                  >  2:451169/1‑139 (MQ=255)                                                                                                                     tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                   <  1:720058/38‑1 (MQ=255)                                                                                                                     tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                   >  2:720058/1‑38 (MQ=255)                                                                                                                      gTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAttt             <  1:1038099/139‑1 (MQ=255)                                                                                                                         aaGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGaa          <  1:18151/139‑1 (MQ=255)                                                                                                                         aaGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGaa          <  2:989280/139‑1 (MQ=255)                                                                                                                           gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa        <  1:676345/139‑1 (MQ=255)                                                                                                                           gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa        <  1:451169/139‑1 (MQ=255)                                                                                                                            agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGaa          <  2:134377/136‑1 (MQ=255)                                                                                                                            agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGaa          >  1:134377/1‑136 (MQ=255)                                                                                                                            agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGAc       >  1:68531/1‑139 (MQ=255)                                                                                                                                tGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTa   >  1:923654/1‑139 (MQ=255)                                                                                                                                 gCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTg                                                <  2:172222/93‑1 (MQ=255)                                                                                                                                 gCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTg                                                >  1:172222/1‑93 (MQ=255)                                                                                                                                 gCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAc  <  2:110767/139‑1 (MQ=255)                                                                                                                                 gCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAc  >  1:172784/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     GTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTAC  >  NC_000913/2316026‑2316292

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