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% 133.2 / ‑3.8 29intergenic (‑91/+144)yciR/rnbmodulator of Rnase II stability/exoribonuclease II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (6/22);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAA  >  REL606/1345128‑1345272
                                                                           |                                                                     
gTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGt                                                                       <  2:749027/76‑1 (MQ=255)
  tAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttt                                                                     <  2:789912/76‑1 (MQ=255)
   aaCAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGGTAAACGCATCGtttt                                                                    >  2:726520/1‑76 (MQ=255)
   aaCAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttt                                                                    <  2:1164029/76‑1 (MQ=255)
     cAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttt                                                                  <  2:558598/76‑1 (MQ=255)
        cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTGAACGCGGTGtgttttttt                                                               >  2:1565666/1‑75 (MQ=255)
        cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttttt                                                               >  1:419927/1‑75 (MQ=255)
          gAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttttat                                                             >  1:1122664/1‑73 (MQ=255)
             tCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATAtatt                                                         <  1:863833/76‑1 (MQ=255)
             tCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATAtatt                                                         <  1:1439971/76‑1 (MQ=255)
               cgcgGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATAtattta                                                       <  1:331986/76‑1 (MQ=255)
               cgcgGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATAtattta                                                       <  1:1224854/76‑1 (MQ=255)
                  ggTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTATatttattt                                                    >  2:623142/1‑76 (MQ=255)
                     taaataaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTa                                                 <  2:1092046/76‑1 (MQ=255)
                         taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATcc                                             <  2:366652/76‑1 (MQ=255)
                         taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATcc                                             <  2:236449/76‑1 (MQ=255)
                         taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATcc                                             <  2:1435028/76‑1 (MQ=255)
                         taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATcc                                             <  2:1283911/76‑1 (MQ=255)
                         taaaAGGAACAACATGACTTACATGAAATTAACGGCGGCCAAACGCCTTATTTTTTTATTTATTTATTTTTAATcc                                             >  2:817546/1‑76 (MQ=255)
                          aaaaGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCATTGCTTTTTTTTATATTTTTTTGTAATTCa                                            >  2:1093848/1‑76 (MQ=255)
                           aaaGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAg                                           <  2:107859/76‑1 (MQ=255)
                                 acaacaTGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTgg                                     <  2:625030/76‑1 (MQ=255)
                                       tGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAAc                               <  1:807204/76‑1 (MQ=255)
                                         aCTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGc                             <  1:1481337/76‑1 (MQ=255)
                                         aCTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGc                             <  1:188598/76‑1 (MQ=255)
                                         aCTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGc                             <  1:436674/76‑1 (MQ=255)
                                          cTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGcc                            <  1:1216618/76‑1 (MQ=255)
                                           ttACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCa                           <  1:1331128/76‑1 (MQ=255)
                                            tACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAg                          <  2:583898/76‑1 (MQ=255)
                                              cATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGtt                        <  1:163431/76‑1 (MQ=255)
                                                tGAAATTAACGGCGGCTAAACGCATTTTTTTTTTTTTTTTTTTTTTTTAATTCAGTTTTTGTAAAAAACCAGtttt                      >  1:844466/1‑76 (MQ=255)
                                                 gAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTc                     <  1:679333/76‑1 (MQ=255)
                                                        aCGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGaa              <  1:1201514/76‑1 (MQ=255)
                                                           gcggcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGt           <  2:770194/76‑1 (MQ=255)
                                                              gcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAg        >  1:107206/1‑76 (MQ=255)
                                                              gcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAg        >  1:692001/1‑76 (MQ=255)
                                                               cTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGt       >  1:804747/1‑76 (MQ=255)
                                                                 aaaCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTa     >  2:860347/1‑76 (MQ=255)
                                                                 aaaCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTa     >  1:493974/1‑76 (MQ=255)
                                                                    cGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGCACAACGCCAGTTTTCAAACGAAAGTCAGTTaaaa  >  2:570717/1‑76 (MQ=255)
                                                                           |                                                                     
GTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAA  >  REL606/1345128‑1345272

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