Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,845,703:1 +G 100% intergenic (‑70/+28) PP_3394 ← / ← PP_3395 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,845,7031.G95.8% 74.7 / ‑2.8 24intergenic (‑70/+28)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (14/9);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845558‑3845842
                                                                                                                                                   |                                                                                                                                           
ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCggg                                                                                                                                          >  6:243762/1‑148 (MQ=255)
  gctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgggct                                                                                                                                        <  2:229305/149‑4 (MQ=255)
      tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAt                                                                                                                                      >  3:176368/1‑148 (MQ=255)
          cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                                 <  5:35174/149‑1 (MQ=255)
             gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                               >  8:188222/1‑148 (MQ=255)
              aTCACTTGCAGGCTCTGCTAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc                                                                                                                              <  4:32250/148‑1 (MQ=255)
               tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCa                                                                                                                             <  7:53594/148‑1 (MQ=255)
                 aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAggc                                                                                                                          >  2:144575/1‑147 (MQ=255)
                        ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                                   <  8:250865/149‑1 (MQ=255)
                                              gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc                                                                                             >  1:80457/1‑149 (MQ=255)
                                                        ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAg                                                                                                              <  5:309013/122‑1 (MQ=255)
                                                        ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAg                                                                                                              >  6:309013/1‑122 (MQ=255)
                                                             tattattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTTATACATGGTGAGCCACCATGGATGAGCCCCCCCCCTGCCCCGCCGGGGGCCgcagca                                                                               >  2:149177/3‑148 (MQ=255)
                                                                       gtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCCGCACCAggtggtt                                                                     >  5:17824/1‑147 (MQ=255)
                                                                                    ctctGTCTTAAGCGTTGCGCCAGAATAATAGTTGTATGAGGGGAAGAAAAGACGTGTTTCGGGGGCTTATCCATCGTCAGCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                        <  2:80457/148‑1 (MQ=255)
                                                                                               gCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                               >  3:278246/1‑146 (MQ=255)
                                                                                                      gcCATACTAAGAGCTGTATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAACCATCGTCACCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                     <  7:34859/149‑1 (MQ=255)
                                                                                                        cAGACTAAGAGCTGGAGGACGGGCAGAAACGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCACCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGGCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                   <  7:210954/149‑1 (MQ=255)
                                                                                                               agagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                       >  7:284611/1‑138 (MQ=255)
                                                                                                                gagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg                           >  8:219816/1‑149 (MQ=255)
                                                                                                                     ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGc                      <  5:243762/149‑1 (MQ=255)
                                                                                                                      gATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                       >  2:302363/1‑147 (MQ=255)
                                                                                                                      gATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                       >  2:297573/1‑147 (MQ=255)
                                                                                                                         gACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                  >  3:30895/1‑149 (MQ=255)
                                                                                                                               cAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGCCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTATGCCGGTCGCGCATGCAGATGACGaa            <  7:188222/149‑1 (MQ=255)
                                                                                                                                         gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  >  7:42866/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845558‑3845842

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