Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,385,797:1	+A	9.3%	coding (81/576 nt)	PP_2093 ←	two‑component system response regulator NasT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,385,797	1	.	A	9.3%	157.3 / 10.3	43	coding (81/576 nt)	PP_2093	two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base . (19/20);  new base A (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAA  >  NC_002947/2385655‑2385940                                                                                                                                                |                                                                                                                                                gCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGA‑GATCATCGAAg                                                                                                                                           >  2:25451/1‑149 (MQ=255)  ccGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGA‑GATCATCGAAGc                                                                                                                                          <  7:70154/149‑1 (MQ=255)   cgcgGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGA‑GATCATCGAAGcc                                                                                                                                         <  1:25451/149‑1 (MQ=255)      ggCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATTGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGc                                                                                                                                      <  3:30235/149‑1 (MQ=255)             cACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCACACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtc                                                                                                                                                  >  7:10941/1‑131 (MQ=255)             cACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCACACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtc                                                                                                                                                  <  8:10941/131‑1 (MQ=255)                          cACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCt                                                                                                                    >  4:360354/1‑147 (MQ=255)                              aCCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAg                                                                                                                 >  8:84850/1‑146 (MQ=255)                              aCCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAg                                                                                                                 <  7:84850/146‑1 (MQ=255)                                  gTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGAcc                                                                                                           >  6:193544/1‑148 (MQ=255)                                       aTCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCGGCCCTCAGGCGACCGAGCtt                                                                                                     <  4:300050/149‑1 (MQ=255)                                             tCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTg                                                                                               >  8:62530/1‑149 (MQ=255)                                                   tCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtc                                                                                         <  3:174332/149‑1 (MQ=255)                                                              ggTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACACAGCCCTCAGGCGACCGAGCTTGctct                                                                                                <  8:355906/131‑1 (MQ=255)                                                              ggTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACACAGCCCTCAGGCGACCGAGCTTGctct                                                                                                >  7:355906/1‑131 (MQ=255)                                                               gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                              <  2:366078/148‑1 (MQ=255)                                                               gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcag                                                                                >  1:366078/1‑146 (MQ=255)                                                                  tCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                                                                          >  3:319076/1‑149 (MQ=255)                                                                         gCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAAc                                                                                                                            >  2:311102/1‑92 (MQ=255)                                                                         gCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAAc                                                                                                                            <  1:311102/92‑1 (MQ=255)                                                                          caccacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACa                                                                  >  2:372614/1‑149 (MQ=255)                                                                             cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcg                                                                                        <  2:353717/124‑1 (MQ=255)                                                                             cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcg                                                                                        >  1:353717/1‑124 (MQ=255)                                                                             cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAAc                                                                   >  6:4526/1‑145 (MQ=255)                                                                             cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAAc                                                                   <  5:4526/145‑1 (MQ=255)                                                                                   cGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTc                                                          <  3:360354/148‑1 (MQ=255)                                                                                    gggCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTccc                                                        <  1:227068/149‑1 (MQ=255)                                                                                        gcACCGTTTCGACGCAGGCGGGCAGGGCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGctc                                                                                                 <  8:166523/104‑1 (MQ=255)                                                                                                  gACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGac                                          <  1:339731/149‑1 (MQ=255)                                                                                                  gACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGa                                           <  5:193544/148‑1 (MQ=255)                                                                                                  gACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGACGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGa                                           >  4:275110/1‑148 (MQ=255)                                                                                                     gCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                              >  5:69485/1‑94 (MQ=255)                                                                                                     gCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                              <  6:69485/94‑1 (MQ=255)                                                                                                                         tCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGccc                    >  5:130098/1‑148 (MQ=255)                                                                                                                         tCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGccc                    >  2:93219/1‑148 (MQ=255)                                                                                                                           aGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACgg                                      >  5:290745/1‑128 (MQ=255)                                                                                                                           aGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACgg                                      <  6:290745/128‑1 (MQ=255)                                                                                                                           aGGTCTGGGGCTTCGACGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCATCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAg                 <  3:275110/149‑1 (MQ=255)                                                                                                                            ggTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAgg                <  6:130098/149‑1 (MQ=255)                                                                                                                             gTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAgg                >  4:276031/1‑148 (MQ=255)                                                                                                                                 ccggCTTCGA‑GATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATg            >  2:113415/3‑148 (MQ=255)                                                                                                                                          atgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGaaa  <  3:276031/149‑1 (MQ=255)                                                                                                                                          atgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGaaa  >  5:203084/1‑149 (MQ=255)                                                                                                                                            gatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGaca                                         >  2:104479/1‑108 (MQ=255)                                                                                                                                            gatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGaca                                         <  1:104479/108‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                GCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAA  >  NC_002947/2385655‑2385940

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