Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,385,836 Δ1 bp 11.8% coding (85/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,8360T.11.8% 97.0 / 8.1 34coding (85/576 nt)PP_2093two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base . (2/2);  total (16/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 = 9.32e-01

ACACGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGC  >  NC_002947/2385689‑2385939
                                                                                                                                                   |                                                                                                        
acacGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACACAGGCGGGTAGGTCGATGGTCAGGTCTGGGGCTTCgat                                                                                                          <  2:519968/148‑1 (MQ=255)
acacGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACACAGGCGGGTAGGTCGATGGTCAGGTCTGGGGCTTCgat                                                                                                          >  1:519968/1‑148 (MQ=255)
 cacGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCgatgat                                                                                                       <  2:269815/150‑1 (MQ=255)
   cGATGGCCCGCGGCTGGGCACGGCTGACCAGGACTACCTGTTCCATCACAGCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCCCCGTGTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtct                                                                                                    <  2:243984/150‑1 (MQ=255)
                    tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctcactc                                                                                   >  1:299735/1‑150 (MQ=255)
                    tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCctcactc                                                                                   >  2:346287/1‑150 (MQ=255)
                    tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGccc                                                                                            >  1:306069/1‑140 (MQ=255)
                    tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGccc                                                                                            <  2:306069/141‑2 (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                 <  1:90190/97‑1 (MQ=255)
                                                                                                                                           ggCTTCGAGATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTc                 >  2:90190/1‑97 (MQ=255)
                                                                                                                                                 gatgatCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                              >  2:105195/1‑78 (MQ=255)
                                                                                                                                                 gatgatCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                              <  1:105195/78‑1 (MQ=255)
                                                                                                                                                   |                                                                                                        
ACACGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGC  >  NC_002947/2385689‑2385939

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: