Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,385,840:1 +A 11.4% 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_0029472,385,8401.A11.4% 115.6 / 11.7 35coding (81/576 nt)PP_2093two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base . (15/16);  new base A (2/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

CGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAACTGCG  >  NC_002947/2385692‑2385988
                                                                                                                                                     |                                                                                                                                                    
cGATGGCCCGCGGCTGGGCACGGCTGACCAGGACTACCTGTTCCATCACAGCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCCCCGTGTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtct                                                                                                                                                     <  2:243984/150‑1 (MQ=255)
                 tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctcactc                                                                                                                                    >  1:299735/1‑150 (MQ=255)
                 tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctcactc                                                                                                                                    >  2:346287/1‑150 (MQ=255)
                 tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGccc                                                                                                                                             <  2:306069/141‑2 (MQ=255)
                 tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGccc                                                                                                                                             >  1:306069/1‑140 (MQ=255)
                   aCGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCaa                                                                                                                                  >  1:427086/1‑150 (MQ=255)
                      ccTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGc                                                                                                                               <  1:45010/149‑1 (MQ=255)
                                  tACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCg                                                                                                                   >  1:598683/1‑150 (MQ=255)
                                     cTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctca                                                                                                                                       <  1:795367/127‑1 (MQ=255)
                                     cTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctca                                                                                                                                       >  2:795367/1‑127 (MQ=255)
                                         tCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGc                                                                                                            <  2:50460/150‑1 (MQ=255)
                                            aTCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTg                                                                                                         >  1:385235/1‑150 (MQ=255)
                                               acaTCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGctc                                                                                                      >  1:738840/1‑150 (MQ=255)
                                                caTCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGctct                                                                                                     <  2:299735/150‑1 (MQ=255)
                                                     gggTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGgtg                                                                                                <  1:36772/149‑1 (MQ=255)
                                                        tCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcg                                                                                             >  2:463595/1‑150 (MQ=255)
                                                          cggCGAGTCGGTATCGCCAAGCACCACGTCCAGGCGCACCGTTTCGACGCAGGCGGGAAGGTCGATGGTAAGGTCTGGGGCCTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcgtc                                                                                           <  1:242993/149‑1 (MQ=255)
                                                                  cGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                                   <  2:385235/150‑1 (MQ=255)
                                                                    gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCaggagg                                                                                 >  1:349974/1‑150 (MQ=255)
                                                                    gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAgg                                                                                 >  1:146542/1‑150 (MQ=255)
                                                                     tATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGa                                                                                >  1:512655/1‑150 (MQ=255)
                                                                       tCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATg                                                                              <  1:580333/150‑1 (MQ=255)
                                                                         gATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATgcg                                                                            <  2:738840/150‑1 (MQ=255)
                                                                          ataaGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCGGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                                   >  2:491902/4‑126 (MQ=255)
                                                                          aTCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCGGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                                   <  1:491902/126‑1 (MQ=255)
                                                                                accacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCa                                                                          <  2:312531/145‑1 (MQ=255)
                                                                                accacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCa                                                                          >  1:312531/1‑145 (MQ=255)
                                                                                  cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                                   >  1:551349/1‑118 (MQ=255)
                                                                                  cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                                   <  2:551349/118‑1 (MQ=255)
                                                                                   acGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTc                                                                  <  2:427086/150‑1 (MQ=255)
                                                                                                ccGCTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCAcgc                                                     <  1:346287/150‑1 (MQ=255)
                                                                                                  gTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACgcgc                                                   <  2:598683/150‑1 (MQ=255)
                                                                                                                                        ggCTTCGAGATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTc                                                                  >  2:90190/1‑97 (MQ=255)
                                                                                                                                        ggCTTCGAGATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTc                                                                  <  1:90190/97‑1 (MQ=255)
                                                                                                                                              gatgatCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                                                                               <  1:105195/78‑1 (MQ=255)
                                                                                                                                              gatgatCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                                                                               >  2:105195/1‑78 (MQ=255)
                                                                                                                                                   tctcGTGGCCGGCCTCACTCAATGCTGCCCTCAGTCGACCGAGCTTGCTCTGGGTGTCGTCGTTCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAACTGCg  <  2:349974/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAACTGCG  >  NC_002947/2385692‑2385988

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