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	98.0%	157.0 / ‑6.2	50	*891* (TAG→TAA)	rcsD	phosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (24/25);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGA  >  NC_000913/2316027‑2316286                                                                                                                                      |                                                                                                                               tATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                            >  1:749431/1‑138 (MQ=255) tATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                            <  2:749431/138‑1 (MQ=255)  aTTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGt                                                                                                                          >  2:1584277/1‑139 (MQ=255)     tGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGCGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGc                                                                                                                       <  1:988697/139‑1 (MQ=255)      gCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGcc                                                                                                                      <  1:239267/139‑1 (MQ=255)       ccATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     <  2:1006780/139‑1 (MQ=255)                     gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     <  1:1278072/125‑1 (MQ=255)                     gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     >  2:1278072/1‑125 (MQ=255)                      tACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                           >  2:1541156/1‑134 (MQ=255)                      tACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                           <  1:1541156/134‑1 (MQ=255)                      tACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAata                                                                                                      <  1:1584277/139‑1 (MQ=255)                           ggCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGCTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                 >  1:595109/1‑139 (MQ=255)                                gCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTa                                                                                                                    <  2:1289597/115‑1 (MQ=255)                                gCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTa                                                                                                                    >  1:1289597/1‑115 (MQ=255)                                  aGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAatta                                                                                          <  1:818621/139‑1 (MQ=255)                                            aaCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGa                                                                                >  2:987522/1‑139 (MQ=255)                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                            >  2:556783/1‑93 (MQ=255)                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                            <  1:556783/93‑1 (MQ=255)                                                    aaCATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGa                                                                        <  1:987522/139‑1 (MQ=255)                                                             ttCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCa                                                                             >  1:1199099/1‑125 (MQ=255)                                                             ttCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCa                                                                             <  2:1199099/125‑1 (MQ=255)                                                               cGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTc                                                             >  1:135318/1‑139 (MQ=255)                                                                 tGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                 <  2:1559616/101‑1 (MQ=255)                                                                 tGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                 >  1:1559616/1‑101 (MQ=255)                                                                   agaAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     <  1:1209860/79‑1 (MQ=255)                                                                   agaAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGACt                                                                                                                     >  2:1209860/1‑79 (MQ=255)                                                                          tGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAAC‑AGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAACAACCGATAGTCTTGTTCGGTATTCGCaaa                                                 >  1:790520/1‑139 (MQ=255)                                                                              ccAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  <  1:1408756/103‑1 (MQ=255)                                                                              ccAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  >  2:1408756/1‑103 (MQ=255)                                                                                  gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                             >  1:18812/1‑72 (MQ=255)                                                                                  gAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                             <  2:18812/72‑1 (MQ=255)                                                                                      aGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  >  2:278922/1‑95 (MQ=255)                                                                                      aGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  <  1:278922/95‑1 (MQ=255)                                                                                       gAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaa                                     >  1:365867/1‑139 (MQ=255)                                                                                        aaaaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaaa                                    <  2:595109/139‑1 (MQ=255)                                                                                        aaaaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaaa                                    >  1:757759/1‑139 (MQ=255)                                                                                             tACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTAt                                                                                                                   <  2:1039923/55‑1 (MQ=255)                                                                                             gaCATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTAt                                                                                                                   >  1:1039923/2‑55 (MQ=255)                                                                                                     cGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     >  2:572249/1‑45 (MQ=255)                                                                                                     cGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                     <  1:572249/45‑1 (MQ=255)                                                                                                     cGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGaa                       >  1:69799/1‑139 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                               <  1:1429645/116‑1 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                               >  2:1429645/1‑116 (MQ=255)                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAAtt        >  2:1071082/1‑139 (MQ=255)                                                                                                                     gTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  >  2:1074786/1‑64 (MQ=255)                                                                                                                     gTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                  <  1:1074786/64‑1 (MQ=255)                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa  <  2:69799/139‑1 (MQ=255)                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa  <  2:757759/139‑1 (MQ=255)                                                                                                                           agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGa                                >  2:375325/1‑108 (MQ=255)                                                                                                                           agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGa                                <  1:375325/108‑1 (MQ=255)                                                                                                                                      |                                                                                                                               TATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGA  >  NC_000913/2316027‑2316286

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