Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,345,210 (T)8→7 intergenic (‑98/+137) yciR ← / ← rnb modulator of Rnase II stability/exoribonuclease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,345,2030T.96.5% 128.2 / ‑2.7 29intergenic (‑91/+144)yciR/rnbmodulator of Rnase II stability/exoribonuclease II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (8/20);  total (9/20)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAA  >  REL606/1345130‑1345274
                                                                         |                                                                       
tAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttt                                                                       <  2:39392/76‑1 (MQ=255)
 aaCAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttt                                                                      <  1:908481/76‑1 (MQ=255)
    aGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttt                                                                   <  2:601813/76‑1 (MQ=255)
      cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttat                                                                 <  2:1131862/76‑3 (MQ=255)
       ggAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCAATGtttttttttt                                                                >  1:472974/1‑74 (MQ=255)
          gTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTatatat                                                            <  1:285867/76‑1 (MQ=255)
           tCCGCGGTTAAATANAAGGAACAACATGACTTACATGAAATTTACGGCGGCTAAACGCATCGTCTTTTTatatatn                                                           <  1:1044210/76‑2 (MQ=255)
             cgcgGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATAtattta                                                         <  1:694821/76‑1 (MQ=255)
                ggTTAAATAAANGGNACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGGCTTTTTATCTatttattt                                                      <  2:190847/76‑1 (MQ=255)
                  ttAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGt                                                    <  2:598913/76‑1 (MQ=255)
                    aaatagaAGGAACAACATGACTTACATGAAATTAACGGCGGGCACAAGATTTGTTTTTTTTTATTTTTTTTTGtta                                                  >  1:80092/1‑74 (MQ=255)
                       taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATGTATTTGTAATcc                                               <  2:46706/76‑1 (MQ=255)
                         aaaGGAACAACATGACTTACATGAAATTAACGGCGGCTAAATGCATTGTTTTTTNATATATTTTTNTGTTATCCAg                                             <  2:1196116/76‑1 (MQ=255)
                             gAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGtttt                                         <  1:708894/76‑1 (MQ=255)
                             gAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGtttt                                         <  1:312649/76‑1 (MQ=255)
                             gAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGtttt                                         <  1:333301/76‑1 (MQ=255)
                                 aacaTGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGaa                                     <  2:431445/76‑1 (MQ=255)
                                 aacaTGACTTACATGAAATTAACGGCGGCCAAAAGCATTGGTTTTTTTTTTTTTTATTTTTTAATCAATTTTTGGa                                      >  1:347697/1‑76 (MQ=16)
                                    aTGACTTACATGAAATTAACGGCGGCTAAAAGCATTGTTTTTTTTTTTTTTTTTTTTTTATTCCGTTTTTGGaaaa                                   >  1:924612/1‑76 (MQ=255)
                                         ttACATGAAATNAANGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCTAGTTTTGGAAATACGCCa                             <  2:185605/76‑1 (MQ=255)
                                             aTGAAATTAACGGCGGCTANANGCATTGTNTTTTTATATATTTGTTTGTAANCCNGTTTTGGAAAAACGCCAGttt                         <  2:175379/76‑1 (MQ=255)
                                              tGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGtttt                        <  2:376937/76‑1 (MQ=255)
                                               gAAATTAACGGNGGNTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTCGGAAAAACTCCAGTTTTc                       <  2:186277/76‑1 (MQ=255)
                                                aaaTTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCa                      <  2:832670/76‑1 (MQ=255)
                                                        ggcggcTAAACGGATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAg              >  2:48039/1‑76 (MQ=255)
                                                          cggcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTc            >  1:557071/1‑76 (MQ=255)
                                                           ggcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa           >  1:909298/1‑76 (MQ=255)
                                                            gcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAg          >  1:408785/1‑76 (MQ=255)
                                                            gcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAg          >  2:510408/1‑76 (MQ=255)
                                                              tAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGtt        >  2:86763/1‑76 (MQ=255)
                                                              tAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGtt        <  1:705189/76‑1 (MQ=255)
                                                               aaaCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTa       >  1:899820/1‑76 (MQ=255)
                                                               aaaCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTa       <  1:625659/76‑1 (MQ=255)
                                                               aaaCCCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGGAAGTCAGTTa       >  1:1114423/1‑76 (MQ=255)
                                                                  cGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTaaaa    <  2:1116924/76‑1 (MQ=255)
                                                                    cATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTaaaaaa  <  1:615269/76‑1 (MQ=255)
                                                                         |                                                                       
TAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAA  >  REL606/1345130‑1345274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: