Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,278,913 T→C 12.1% D23D (GAT→GAC PP_3750 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,278,9130TC12.1% 80.9 / 3.3 33D23D (GAT→GACPP_3750GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (13/16);  new base C (2/2);  total (15/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 = 8.11e-01

TCTTGGTAGCCTGTGACAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGT  >  NC_002947/4278769‑4279050
                                                                                                                                                 |                                                                                                                                         
tctTGGTAGCCTGTGACAAGCACCTGAGGCGCAGCCAATCCGGTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAAGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTAt                                                                                                                                      >  2:458110/1‑150 (MQ=255)
 ctTGGTAGCCTGTGACAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATc                                                                                                                                     <  2:386740/150‑1 (MQ=255)
          cTGTGACAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGc                                                                                                                            >  1:200056/1‑150 (MQ=255)
           tgtgACAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGcc                                                                                                                           >  1:425890/1‑150 (MQ=255)
               aCAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCt                                                                                                                       >  1:93703/1‑150 (MQ=255)
                        tGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCgg                                                                                                              <  2:241224/150‑1 (MQ=255)
                         gAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGACGGTATCGACAGCGGCCACCTGGCCGCCGGc                                                                                                             >  2:151527/1‑150 (MQ=255)
                                   cAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTgc                                                                                                   >  1:186413/1‑150 (MQ=255)
                                     aTCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTgg                                                                                                                     >  1:355256/1‑130 (MQ=255)
                                     aTCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTgg                                                                                                                     <  2:355256/130‑1 (MQ=255)
                                     aTCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTCCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTgccg                                                                                                 >  2:558667/1‑150 (MQ=255)
                                             tCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCg                                                                                         >  2:490050/1‑150 (MQ=255)
                                                tGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCTGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTct                                                                                      >  2:91938/1‑150 (MQ=255)
                                                 gAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGACGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTctg                                                                                     <  1:151527/150‑1 (MQ=255)
                                                          tccgctccgGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGaa                                                                            >  2:63803/1‑150 (MQ=255)
                                                                   gACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGc                                                                   <  1:63803/150‑1 (MQ=255)
                                                                        gCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCAACCTGGCCGCCGGCTCCCAACTGCCCCCTATCTGTCTTCTGGCTGAACAACTGTGCGTATc                                                              <  2:93646/150‑1 (MQ=255)
                                                                          ccATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATcgc                                                            <  1:560639/150‑1 (MQ=255)
                                                                           cATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATcgcg                                                           >  1:354749/1‑150 (MQ=255)
                                                                             tGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCa                                                         >  2:54824/1‑150 (MQ=255)
                                                                                aaCTGCACATGCGCCTGGAGGGCCGCAAAGGCCTGGCCGCTACGCTATGCTGAATGAATGATCAGTTGAGCCACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGa                                                      <  1:490050/150‑1 (MQ=255)
                                                                                      acaTTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTgg                                                <  1:558667/150‑1 (MQ=255)
                                                                                          tCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGa                                            <  2:200056/150‑1 (MQ=255)
                                                                                           cGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAg                                           <  2:425890/150‑1 (MQ=255)
                                                                                                       cGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATACTGCGCTGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCgg                               <  1:458110/150‑1 (MQ=255)
                                                                                                                      gATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTgg                                                <  2:129867/118‑1 (MQ=255)
                                                                                                                      gATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTgg                                                >  1:129867/1‑118 (MQ=255)
                                                                                                                      gATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATAGCGCACGACCCTCGCCTAGGTCGATGTGCGGCTTGCCTTCGATatc                <  2:93703/150‑3 (MQ=255)
                                                                                                                       aTCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAAcc               <  1:525501/150‑1 (MQ=255)
                                                                                                                         cAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTg             <  1:54824/150‑1 (MQ=255)
                                                                                                                                    aGCAATTGCGCGACGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCt           >  2:166735/1‑142 (MQ=255)
                                                                                                                                    aGCAATTGCGCGACGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCt           <  1:166735/142‑1 (MQ=255)
                                                                                                                                     gCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGt  <  2:538488/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
TCTTGGTAGCCTGTGACAAGCACCTGAGTCGCAGCCAATCCGCTATCCTGAGCAATCCTCCGCTCCGGACCCGCCCATGGAACTGCACATTCGCCTGGAGGGCCGCAAAGGCCTGGCCGATCAGCTTTACCAGCAATTGCGCGATGGTATCGACAGCGGCCACCTGGCCGCCGGCTCCCAACTGCCGCCTACCCGTCTGCTGGCTGAACAACTGGGCGTATCGCGCAAGACCGTGGCCGAGGCCTATTCGCGGCTGACCTACGATAACCTGCTCAGCGGCGT  >  NC_002947/4278769‑4279050

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: