Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,287,319 T→C 16.7% intergenic (+17/+24) PP_2889 → / ← xthA transmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,287,3190TC16.7% 64.6 / 6.5 30intergenic (+17/+24)PP_2889/xthAtransmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (3/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

GGACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGG  >  NC_002947/3287170‑3287460
                                                                                                                                                     |                                                                                                                                              
ggAGGGGCGGCGGGGGGGTGGCGGGCGGGGGGGGGGGTACTGGTGGTATGGGGTTTATATTTTGGCGGGGGTGAGGGCGGCGGGGGGGGGGCGTGGGGGGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGgtgt                                                                                                                                                <  1:1501639/150‑1 (MQ=255)
ggACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGgtgt                                                                                                                                                >  1:1847263/1‑150 (MQ=255)
      gCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGTGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCt                                                                                                                                          >  1:244777/1‑150 (MQ=255)
           gggTCGGTGCCCAGCGGGGGGCGCGGGACTGGTCGCGTGGGTCTTATCACTTGGCCTTGGTCAGCGCGGGGGACGGGGTGCGGGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGc                                                                                                                                     <  1:1704548/150‑1 (MQ=255)
                                ggggggaCTGGTGGCATGGGGCTTATCGCTTTGCCTTGGTCAGCGGGGCGGACGATGGGGGTGGGGGGGGGGTGGGGATGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCaa                                                                                                               <  1:1379769/144‑1 (MQ=255)
                                       cTGGTGGCATGGGCTTTTTATGTTTTCCTTGGTAGGCGGGGGGGACGGGGGGGGGGGGGCGGGGGTGGAGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTccagc                                                                                                        <  2:773995/150‑1 (MQ=255)
                                        tGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCTCGGCGGACGATGTGCGAGGGGCGGGGAAGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGCGCCCGCCACCCTAATCAGCTCAACTccagcc                                                                                                       >  2:628093/1‑150 (MQ=255)
                                          gTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTccagccag                                                                                                     >  2:1650397/1‑150 (MQ=255)
                                          gTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGAAGGCGTTGAGGTTGGGGTTGTACGTGTGGGCCCAAGGTGTGCACCTTGGTCCCGTACCCCTAATCAGCTCAATTccagccag                                                                                                     >  2:7721/1‑150 (MQ=255)
                                               aTGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCgg                                                                                                >  1:1749434/1‑150 (MQ=255)
                                                tGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGTGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGgc                                                                                               <  2:244777/150‑1 (MQ=255)
                                                        ataACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTc                                                                                       >  1:1582351/1‑150 (MQ=255)
                                                              gggCCTTGGGCAGCGCGGCGGACGTTTGTTGGGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACgg                                                                                 <  2:1365225/148‑1 (MQ=255)
                                                                            gcggcgGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCcgc                                                                   >  1:166143/1‑150 (MQ=255)
                                                                                cggGCGAGGTGCGGGGGGCGGGGGTGGGGGAGGTGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAg                                                               <  1:1251020/150‑1 (MQ=255)
                                                                                  gACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGACTCCAGTCACTATTTCAGCTCAACTCTAGCCAGATCGCCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCt                                                             >  1:1006054/1‑150 (MQ=255)
                                                                                         gcggtggtgGGGGGGTGGTGGAGGCGTTGAGGGGGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCATCTCCAGGTAGAGCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTc                                                      <  2:758300/147‑1 (MQ=255)
                                                                                               gggCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTACAAATTGATGACCCAAGGTGTGCACCTTGGGCCCTGCTCCCTAATACGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCACCTCATATACCACGcc                                                >  1:1107142/1‑150 (MQ=255)
                                                                                                  cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGc                                             >  2:1841984/1‑150 (MQ=255)
                                                                                                  cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGccca                                                                                                                                  <  2:289699/66‑2 (MQ=25)
                                                                                                  cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGccca                                                                                                                                  >  1:289699/1‑65 (MQ=25)
                                                                                                      ggggggcgAGGGGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCCTGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGcc                                         <  2:166143/147‑1 (MQ=255)
                                                                                                                   ttGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCa                            >  2:1816284/1‑150 (MQ=255)
                                                                                                                              gTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCtgt                 <  2:1582351/150‑1 (MQ=255)
                                                                                                                                    aGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGcc           <  2:1107142/150‑1 (MQ=255)
                                                                                                                                       gtgGTCCCAAGGTGTGCACCTTGGGCCCGTCATCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATg        <  2:144381/150‑1 (MQ=255)
                                                                                                                                         gggCCCAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAt      >  2:1157158/1‑150 (MQ=255)
                                                                                                                                         gggCCCAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAt      >  2:1010698/1‑150 (MQ=255)
                                                                                                                                            cccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa               <  2:156420/138‑1 (MQ=255)
                                                                                                                                            cccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa               >  1:156420/1‑138 (MQ=255)
                                                                                                                                            cccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAg   <  2:1444277/150‑1 (MQ=255)
                                                                                                                                             ccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAgg  >  1:954906/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGG  >  NC_002947/3287170‑3287460

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