Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,846,798 T→C 9.5% intergenic (‑102/‑113) PP_3395 ← / → PP_3396 LysR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,846,7980TC 9.5% 109.9 / 3.1 42intergenic (‑102/‑113)PP_3395/PP_3396LysR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base C (2/2);  total (20/22)
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.80e-01

GGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTGGCAGTGG  >  NC_002947/3846653‑3846944
                                                                                                                                                 |                                                                                                                                                  
ggCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTa                                                                                                                                                   <  8:198732/147‑1 (MQ=255)
   ttGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGcccc                                                                                                                                              >  5:326991/1‑149 (MQ=255)
        gAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCAc                                                                                                                                         >  2:249751/1‑149 (MQ=255)
         aaGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCAc                                                                                                                                         >  7:77577/1‑148 (MQ=255)
                  ccgTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCt                                                                                                                               <  8:267154/149‑1 (MQ=255)
                       gCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCaa                                                                                                                           <  1:249751/148‑1 (MQ=255)
                                cAAATACCGGATTGGCTTCTAAATCCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGTGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTc                                                                                                                 <  3:338304/149‑1 (MQ=255)
                                cAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGACCAACTCCAGAGGCTGGCATAAAGACCTACTCAGAATCTGTAGCCCCTGCACCAAAAGGCCTGCAAAGCAttttt                                                                                                                  >  3:151939/1‑148 (MQ=255)
                                 aaaaaCCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCa                                                                                                                <  1:318650/149‑1 (MQ=255)
                                 aaaaaCCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCa                                                                                                                <  7:352099/149‑1 (MQ=255)
                                 aaaaaCCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCa                                                                                                                <  4:137948/149‑1 (MQ=255)
                                         cATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTgg                                                                                                        <  5:8501/149‑1 (MQ=255)
                                                  cTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGACCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGctct                                                                                               >  7:341454/1‑149 (MQ=255)
                                                     aaaCCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGc                                                                                            <  5:348978/149‑1 (MQ=255)
                                                        ccGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGa                                                                                         >  4:128396/1‑149 (MQ=255)
                                                         cGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGAt                                                                                        <  4:178666/149‑1 (MQ=255)
                                                                 aaTTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCtcctcc                                                                                 >  6:119843/1‑148 (MQ=255)
                                                                  aTTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCtcctcc                                                                                 >  8:315363/1‑147 (MQ=255)
                                                                  aTTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCtcctcc                                                                                 >  6:139799/1‑147 (MQ=255)
                                                                     aCCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGa                                                                             <  3:128396/148‑1 (MQ=255)
                                                                       cGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGc                                                                          <  4:61443/149‑1 (MQ=255)
                                                                       cGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGc                                                                          <  6:326991/149‑1 (MQ=255)
                                                                                 tCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCa                                                                <  7:315363/149‑1 (MQ=255)
                                                                                     cAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCAAGACCCACCAAGTCCATggg                                                            >  2:320983/1‑149 (MQ=255)
                                                                                              gggAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTg                                                    >  6:368687/1‑148 (MQ=255)
                                                                                                 aaCGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCag                                                 >  3:173271/1‑148 (MQ=255)
                                                                                                   cGCTGCCCACTTCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCaggag                                              <  1:320983/149‑1 (MQ=255)
                                                                                                     cTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCaggagg                                             >  8:263149/1‑148 (MQ=255)
                                                                                                          cACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGCAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGa                                       >  7:186242/1‑149 (MQ=255)
                                                                                                          cACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGCAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGa                                       >  8:242335/1‑149 (MQ=255)
                                                                                                               ccAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCAt                                  >  3:370406/1‑149 (MQ=255)
                                                                                                               ccAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCAt                                  <  8:297289/149‑1 (MQ=255)
                                                                                                                          cATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCaccgaccg                        >  6:128645/1‑148 (MQ=255)
                                                                                                                              aatgccTACTGAGAATCTGTAGCACCTGCACCATAAGGCCTGCAAAGGATTTTTCTGGTGTTGTCAGATCTCTTGCCGATCTCCTCCCAGATGCCCCTAGTTCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACTGACCGTTAcc                   <  4:151939/145‑1 (MQ=255)
                                                                                                                                  ccTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCAc               <  8:77577/149‑1 (MQ=255)
                                                                                                                                    tACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCAccc             >  4:116868/1‑149 (MQ=255)
                                                                                                                                       tGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCTGACGCCCCAAGTCCATgg                                                             <  3:41028/98‑1 (MQ=255)
                                                                                                                                       tGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCTGACGCCCCAAGTCCATgg                                                             >  4:41028/1‑98 (MQ=255)
                                                                                                                                         agaATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTgg        <  5:128645/149‑1 (MQ=255)
                                                                                                                                           aaTCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGGTTGCAGGAGGATGTGCCTCCTGCACCGCCCGATACCCCCCCCGCggg        >  8:80445/1‑147 (MQ=255)
                                                                                                                                            aTCTGCAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTGGCa      <  7:242335/148‑1 (MQ=255)
                                                                                                                                            aTCTGCAGCCCATGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTGGCa      <  8:186242/148‑1 (MQ=255)
                                                                                                                                                 tAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTGGCAGTgg  <  3:116868/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCTTGGGGAAGGGACCGCCGTAGCACTGTAGCAAAAACCGCATTGGCTTCTAAAACCGGGTTTGAATTACCGGCCCCCAATCCCCAACCTTGGGGGAACGCTGCCCACCGCCAGCGGCTGGCATAAAGCCCTACTGAGAATCTGTAGCCCCTGCACCATAAGGCCTGCAAAGCATTTTTCAGGTGTTGGCACAGCTCTTGCCGATCTCCTCCCAGACGCCCCAAGTCCATGGGCCCGTTTGCAGGAGGAATTGACTCATGCACCGACCGATACCCCACCCGCTGGCAGTGG  >  NC_002947/3846653‑3846944

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