Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,087,485 G→T 11.1% R264M (AGG→ATG)  PP_2695 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,087,4850GT11.1% 100.9 / 3.7 36R264M (AGG→ATG) PP_2695LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base T (2/2);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GGCTGTTCATTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGATTCGTTGGCGGTAGGGGCATGGCGCGCCAT  >  NC_002947/3087356‑3087621
                                                                                                                                 |                                                                                                                                        
ggCTGTTCATTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCATGGTAGAACGGCATATCGGCa                                                                                                                      >  2:252326/1‑150 (MQ=255)
        aTTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAggg                                                                                                                                        <  2:527173/124‑1 (MQ=255)
        aTTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAggg                                                                                                                                        >  1:527173/1‑124 (MQ=255)
        aTTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCATGGTAGAACGGc                                                                                                                               <  2:837372/133‑1 (MQ=255)
        aTTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCATGGTAGAACGGc                                                                                                                               >  1:837372/1‑133 (MQ=255)
         tttCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACAGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCAccc                                                                                                             >  2:138652/1‑150 (MQ=255)
                   gATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCgg                                                                                                                        <  2:662063/129‑1 (MQ=255)
                   gATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCgg                                                                                                                        >  1:662063/1‑129 (MQ=255)
                     tttCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCa                                                                                                 <  2:428311/150‑1 (MQ=255)
                             aTCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCa                                                                                         <  2:305690/150‑1 (MQ=255)
                              tCGCCATGCTGGAGGAGGGGGTAGGGAGCAGTACCTTGTCGTCCCGCGGTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCGCTAGAGGAGCCCAGGGTAGATCGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCAc                                                                                        <  1:138652/150‑1 (MQ=255)
                                gCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCt                                                                                      >  2:807354/1‑150 (MQ=255)
                                                        aTCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCATGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGttt                                                              <  1:252326/150‑1 (MQ=255)
                                                         tCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCCCCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTc                                                             <  2:386596/150‑1 (MQ=255)
                                                         tCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGt                                                                                           >  1:170957/1‑120 (MQ=255)
                                                         tCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGt                                                                                           <  2:170957/120‑1 (MQ=255)
                                                              aCCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGAAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCaa                                                        <  2:302873/150‑1 (MQ=255)
                                                                cTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCtt                                                                                             <  1:162479/111‑1 (MQ=255)
                                                                cTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCtt                                                                                             >  2:162479/1‑111 (MQ=255)
                                                                cTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTAGCGTTCGTACCCCTTTAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTTTGCAAGt                                                      >  2:375877/1‑150 (MQ=255)
                                                                    ccGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTc                                                             >  2:635911/1‑139 (MQ=255)
                                                                    ccGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTc                                                             <  1:635911/139‑1 (MQ=255)
                                                                         ccTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGtt                                                               >  2:189753/1‑132 (MQ=255)
                                                                         ccTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGtt                                                               <  1:189753/132‑1 (MQ=255)
                                                                           tCGCTTTTCCGCAGGAGAGCGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGa                                           >  1:333163/1‑150 (MQ=255)
                                                                              cTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAgg                                        >  1:150869/1‑150 (MQ=255)
                                                                                      cAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGa                                <  2:150869/150‑1 (MQ=255)
                                                                                       aGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGAt                               <  1:391035/150‑1 (MQ=255)
                                                                                              cGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTNATGGCGNTTCTGCAAGTGCNGTTGGCGAAGGCTGGCTGATTCGTTgg                        <  1:807354/150‑1 (MQ=255)
                                                                                                      tGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCAcc                                                                                                              >  2:617810/1‑56 (MQ=255)
                                                                                                      tGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCAcc                                                                                                              <  1:617810/56‑1 (MQ=255)
                                                                                                                   ccTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGATTCGTTGGCGGTAGGGGCAt            <  2:273404/141‑1 (MQ=255)
                                                                                                                   ccTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGATTCGTTGGCGGTAGGGGCAt            >  1:273404/1‑141 (MQ=255)
                                                                                                                    cTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTCCATACCGGGTGCTCGATGGCTGGCTGATTCGTTGGCGGTAGGGGCATGGCGCGCCAt  >  2:450745/1‑150 (MQ=255)
                                                                                                                            gCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGaa                                          >  2:100691/1‑102 (MQ=255)
                                                                                                                            gCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGaa                                          <  1:100691/102‑1 (MQ=255)
                                                                                                                             cccAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATg                                                                    >  2:354835/1‑75 (MQ=255)
                                                                                                                             cccAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATg                                                                    <  1:354835/75‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                        
GGCTGTTCATTTCCAACATGATTTCGTTGATCGCCATGCTGGAGGAGGGGGTAGGGATCAGTACCTTGCCGTACCTCGCTTTTCCGCAGGAGAACGACAAACTGGCGTTCGTACCCCTTGAGGAGCCCAGGGTAGAACGGCATATCGGCATGCTCACCCGCAAAGGGCGCAGCTTGTCACCTGTGGCGGGGGCACTGATGGCGTTTCTGCAAGTGCAGTTGGCGAAGGCTGGCTGATTCGTTGGCGGTAGGGGCATGGCGCGCCAT  >  NC_002947/3087356‑3087621

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