Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 830,378 T→C Q89R (CAG→CGG)  cecR ← DNA‑binding transcriptional dual regulator CecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913830,3780TC96.3% 85.7 / NA 27Q89R (CAG→CGG) cecRDNA‑binding transcriptional dual regulator CecR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/8);  minor base A (1/0);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGG  >  NC_000913/830263‑830504
                                                                                                                    |                                                                                                                              
gagCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAg                                                                                                               >  2:605834/1‑133 (MQ=255)
gagCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAg                                                                                                               >  1:395114/1‑133 (MQ=255)
      ctCACGGGAGATAAACTTGCTGAGGTTGGCGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGGTCAGTTCACGGGTGGCAGCCGGATTCAGACGTGGTTGAGCGACAAAGCGttag                                                                                                         >  2:299126/1‑131 (MQ=255)
           gggAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCCGCCCGATCCGGCAGTGCTTGCGCGAACAAGGGTTCGGCTTc                                                                                                    >  2:199353/1‑133 (MQ=255)
                   aaCTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCAt                                                                                            <  2:570285/133‑1 (MQ=255)
                   aaCTTGCTGAGGTTGACGGTGTCACCCTGGGTCAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCAt                                                                                            <  2:186367/133‑1 (MQ=255)
                        gCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcgga                                                                                       >  2:577138/1‑133 (MQ=255)
                        gCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcgga                                                                                       >  1:578262/1‑133 (MQ=255)
                             ggTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGaa                                                                                  >  1:324695/1‑133 (MQ=255)
                                  aCGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTCCGGCTTCGTCGTCATGCGGACGGAACTGCt                                                                             >  1:249633/1‑133 (MQ=255)
                                        tCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCaa                                                                       >  2:85419/1‑133 (MQ=255)
                                        tCATCCAGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCaa                                                                       >  2:455808/1‑133 (MQ=255)
                                              tGGGTGAGCAGCTTATTCTTGTTCCTGCAGTCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAAt                                                                 <  1:199353/133‑1 (MQ=255)
                                                ggTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCg                                                               >  1:39743/1‑133 (MQ=255)
                                                ggTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCg                                                               >  1:376701/1‑133 (MQ=255)
                                                ggTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCCGTTGGAACTGTTCGCCAATAAAAtag                                                               >  2:287930/1‑131 (MQ=255)
                                                  tGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGc                                                             >  1:399983/1‑133 (MQ=255)
                                                         cTTAATCACGTTCCTGCATGCGTGAAGGATAAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCAc                                                      <  1:287930/133‑1 (MQ=255)
                                                                   ttCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACAGAACTGCTCGCCGATAAAATCGGCAGTCCACTGGGCGCAgg                                            >  2:260481/1‑133 (MQ=255)
                                                                    tCCTGCAGGCACGAAGGATCAGTTTACGGATGGCAGCCCGAACAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGc                                           <  2:246498/133‑1 (MQ=255)
                                                                               cgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAt                                >  1:342708/1‑133 (MQ=255)
                                                                               cgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAt                                >  1:20653/1‑133 (MQ=255)
                                                                               cgAAGGATCAATTCACGGATGGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAt                                >  1:495294/1‑133 (MQ=255)
                                                                                         gTTCACGGCTGGCAGCCCAATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGa                      <  1:299126/133‑1 (MQ=255)
                                                                                              cGGATGGCTGCCCGATCATGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                 <  1:257764/133‑1 (MQ=255)
                                                                                              cGGATAGCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                 >  2:78972/1‑133 (MQ=255)
                                                                                                    gCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTTGAACCGAag               >  2:394180/1‑129 (MQ=255)
                                                                                                    gCAGCCCGATCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAag               <  1:394180/129‑1 (MQ=255)
                                                                                                              tCAGGCCGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATgg  <  1:605834/133‑1 (MQ=255)
                                                                                                                    |                                                                                                                              
GAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGG  >  NC_000913/830263‑830504

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