Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,886,181 A→G 12.1% intergenic (+60/+63) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,886,1810AG12.1% 82.7 / 3.6 33intergenic (+60/+63)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base G (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

AGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTA  >  NC_002947/4886033‑4886325
                                                                                                                                                    |                                                                                                                                                
aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                                  <  2:200191/149‑1 (MQ=255)
    cAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                              <  4:59538/149‑1 (MQ=255)
       gCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                           <  1:153481/149‑1 (MQ=255)
         tGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACACGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         <  1:16067/149‑1 (MQ=255)
              gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTccc                                                                                                                                    >  7:260735/1‑149 (MQ=255)
                  gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                <  1:205957/149‑1 (MQ=255)
                  gTGATTAGGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                <  8:128364/149‑1 (MQ=255)
                        aTGACGCGGCGGCCGAGACCATCATTCGTATTGTCTTGAAAGGTTGTGAGCCTGAGTCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                           <  8:99742/148‑1 (MQ=255)
                         tGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         <  6:56511/133‑1 (MQ=255)
                         tGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         >  5:56511/1‑133 (MQ=255)
                           aCGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCggg                                                                                                                        <  6:257419/148‑1 (MQ=255)
                             gcggcgGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAAACATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                            <  3:184267/126‑1 (MQ=255)
                             gcggcgGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAAACATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                            >  4:184267/1‑126 (MQ=255)
                              cggcggCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTg                                                                                                                    <  3:28623/149‑1 (MQ=17)
                                 cggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGccc                                                                                                                 <  8:260735/149‑1 (MQ=255)
                                           catcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         >  1:173319/1‑115 (MQ=255)
                                           catcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         <  2:173319/115‑1 (MQ=255)
                                                                     gtgAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTACCACAGGATCACCACATGCCTGAAGGGTGTGGAGTCCCTGTGGGAGCGCGCTTGCCCGCGAATGGGCCAGAAGAGGCAATAGAGATTTAGGcc                                                                             >  5:307782/1‑149 (MQ=255)
                                                                                                       ccATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAcctcct                                           >  5:76713/1‑149 (MQ=255)
                                                                                                         aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGCATGGGGCAGAACAGGCAATAGAGCTTTAGGCCGCGACCCCGGCATCTGCCCGCCGCCCCCcctccgc                                          >  2:310924/1‑146 (MQ=38)
                                                                                                                   aaGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGa                               >  7:210794/1‑149 (MQ=255)
                                                                                                                       ccGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCg                            >  5:181830/1‑148 (MQ=18)
                                                                                                                       ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCg                            >  6:54998/1‑148 (MQ=255)
                                                                                                                         gCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCg                                     >  7:30873/1‑136 (MQ=17)
                                                                                                                         gCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCg                                     <  8:30873/136‑1 (MQ=17)
                                                                                                                                   gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtc               <  1:310924/149‑1 (MQ=255)
                                                                                                                                    aTCACAACATGCCTGAAGGTTGTGGAGTCACTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAATTAGACCGCTACCCCg                                                                    >  7:226981/1‑95 (MQ=38)
                                                                                                                                        ccacATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTAccc                                                                      >  5:283370/1‑89 (MQ=255)
                                                                                                                                        ccacATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTAccc                                                                      <  6:283370/89‑1 (MQ=255)
                                                                                                                                          acaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        <  6:76713/149‑1 (MQ=255)
                                                                                                                                            aTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa      >  2:179853/1‑149 (MQ=255)
                                                                                                                                                cTGAAGGTTTTGGAGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTa  <  6:307782/149‑1 (MQ=255)
                                                                                                                                                 tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                                <  5:221428/102‑1 (MQ=255)
                                                                                                                                                 tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                                >  6:221428/1‑102 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTA  >  NC_002947/4886033‑4886325

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: