Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,941,335 C→A 11.7% R189R (CGC→CGA PP_2574 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,941,3350CA11.7% 95.1 / 3.5 34R189R (CGC→CGAPP_2574LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (16/14);  new base A (2/2);  total (18/16)
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.53e-01

ACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATGGTCGGTCAGG  >  NC_002947/2941198‑2941480
                                                                                                                                         |                                                                                                                                                 
aCATGCGACCAGGCATTTCTCACTGACCACGACCTGGATGGCACTCTCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGccc                                                                                                                                        <  4:218815/149‑1 (MQ=255)
              aTTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCAGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTCTGCCCGGAGCCGATGATc                                                                                                                           >  5:350402/1‑148 (MQ=255)
                  ctctATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCctgc                                                                                                                       <  6:354172/148‑1 (MQ=255)
                     tATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGa                                                                                                                   >  5:248386/1‑149 (MQ=255)
                         aCCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTg                                                                                                               >  8:148221/1‑149 (MQ=255)
                           cACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGc                                                                                                             >  1:322452/1‑149 (MQ=255)
                                                gAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGt                                                                                        <  1:303107/149‑1 (MQ=255)
                                                        gCCGCTGATGCCGCCTCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTg                                                                                 <  6:248386/148‑1 (MQ=255)
                                                          cGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCACACGTACTTTGTCAGCTTGTTc                                                                              <  8:216887/149‑1 (MQ=255)
                                                               aTGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAAc                                                                           <  6:350402/147‑1 (MQ=255)
                                                                   cgccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGtt                                                                               <  6:177523/139‑1 (MQ=255)
                                                                   cgccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGtt                                                                               >  5:177523/1‑139 (MQ=255)
                                                                    gccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATgg                                                                    >  6:324704/1‑149 (MQ=255)
                                                                    gccgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATgg                                                                    >  3:155577/1‑149 (MQ=255)
                                                                      cgcAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGc                                                                  <  5:324704/149‑1 (MQ=255)
                                                                                            aCCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTg                                                               >  2:149424/1‑130 (MQ=255)
                                                                                            aCCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTg                                                               <  1:149424/130‑1 (MQ=255)
                                                                                                 agaACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGccc                                                                                                           <  7:140822/81‑1 (MQ=255)
                                                                                                 agaACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGccc                                                                                                           >  8:140822/1‑81 (MQ=255)
                                                                                                     ccaccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCtt                                              >  1:208710/1‑138 (MQ=255)
                                                                                                     ccaccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCtt                                              <  2:208710/138‑1 (MQ=255)
                                                                                                        ccGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcga                                >  1:345678/1‑149 (MQ=255)
                                                                                                          gCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcgatc                              >  2:150072/1‑149 (MQ=255)
                                                                                                            tttGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGtcgatc                              >  1:81418/1‑147 (MQ=255)
                                                                                                              tGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCgag                           <  4:155577/148‑1 (MQ=255)
                                                                                                                   tccgggCACAGGTCGCGCAGCGAGACCTGGGCCTGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTg                     >  3:326162/5‑149 (MQ=38)
                                                                                                                       ggCACAGGTCGCGCAGCGAGACCTGGGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGcgcg                 >  3:198392/1‑149 (MQ=255)
                                                                                                                          ccAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCa              >  5:232796/1‑149 (MQ=255)
                                                                                                                          acAGGTCGCGCAGCGAGACCTGGGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCa              <  4:198392/148‑1 (MQ=255)
                                                                                                                            aGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATg            >  4:16371/1‑149 (MQ=255)
                                                                                                                             ggTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATgg           >  2:282795/1‑149 (MQ=255)
                                                                                                                             ggTCGCGCAGCGAGACCTGGGCCTGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATgg           <  4:326162/149‑1 (MQ=255)
                                                                                                                              gTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATggt          <  2:345678/149‑1 (MQ=255)
                                                                                                                                   gCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATggtcg        <  3:16371/146‑1 (MQ=255)
                                                                                                                                       cgcgACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATGGTCGGTCAgg  <  2:81418/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
ACATTCGACCTGGCATTTCTCTATGACCACGACCTGGATGGCACTATCGAGGCCGAGCCGCTGATGCCGCCGCAGAAGCCCTATGTGCTGCTACCCGAGAACCACCGCTTTGCCAGCCAGGCCCAGGTCTCGCTGCGCGACCTGTGCCCGGAGCCGATGATCCTGCTGGACGTGGCCCCTAGCCGCACGTACTTTGTCAGCTTGTTCAACGAGATGGGCCTGGCACCGAACATCGTCTTCAGTTCACCGTCGATCGAGATGGTGCGCGGCATGGTCGGTCAGG  >  NC_002947/2941198‑2941480

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