Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,332,148 (T)6→5 intergenic (+131/‑79) topA → / → cysB DNA topoisomerase I/DNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,332,1430T.96.2% 239.7 / 4.1 52intergenic (+126/‑84)topA/cysBDNA topoisomerase I/DNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (27/23);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGT  >  REL606/1332068‑1332214
                                                                            |                                                                       
tttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCt                                                                         <  1:1425333/76‑1 (MQ=255)
  tttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttt                                                                       >  2:1408538/1‑76 (MQ=255)
  tttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttt                                                                       >  1:1654222/1‑76 (MQ=255)
   ttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttt                                                                      <  1:1890920/76‑1 (MQ=255)
    tATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttttt                                                                     <  1:670039/76‑1 (MQ=255)
    tATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttttt                                                                     <  2:1782629/76‑1 (MQ=255)
          atatTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTttattatt                                                              >  1:1855963/1‑76 (MQ=255)
           tattatTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTttattatta                                                             >  2:890444/1‑76 (MQ=255)
           tattatTTTTTCTTCACACCTATACACTAAGGCGATAAATGATATAGTGGTTATAGTTAGCACCTTTttattatta                                                             <  2:124762/76‑1 (MQ=255)
            attattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTaa                                                            <  2:233875/76‑1 (MQ=255)
             ttattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTaa                                                            <  1:904842/76‑1 (MQ=255)
              tattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAAt                                                          <  1:512266/76‑1 (MQ=255)
                ttttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCg                                                        <  1:1536371/76‑1 (MQ=255)
                 tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       <  1:1743694/76‑1 (MQ=255)
                 tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       <  1:628073/76‑1 (MQ=255)
                 tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       <  1:584611/76‑1 (MQ=255)
                 tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       <  1:508400/76‑1 (MQ=255)
                 tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTAGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       >  1:1363178/1‑76 (MQ=255)
                   tttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAt                                                     <  1:46058/76‑1 (MQ=255)
                     tcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTa                                                   >  1:1189439/1‑76 (MQ=255)
                     tcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTa                                                   >  1:304051/1‑76 (MQ=255)
                       ttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCTCCTGNTTATTATTAAATCGTATTAGt                                                 >  2:316909/1‑76 (MQ=255)
                       ttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGt                                                 <  2:1177614/76‑1 (MQ=255)
                         cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGt                                                 <  1:380523/76‑1 (MQ=255)
                         cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCa                                               <  2:1566779/76‑1 (MQ=255)
                         cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCa                                               <  2:134773/76‑1 (MQ=255)
                         cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCa                                               <  1:1371066/76‑1 (MQ=255)
                           cacCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCAcc                                             >  2:1468340/1‑76 (MQ=255)
                             ccTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCg                                           <  1:1826071/76‑1 (MQ=255)
                               tataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGcc                                         <  2:615062/76‑1 (MQ=255)
                                ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCa                                        <  2:1894208/76‑1 (MQ=255)
                                 taCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAg                                       <  2:651774/76‑1 (MQ=255)
                                  acacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAgg                                      <  2:423535/76‑1 (MQ=255)
                                   cacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgt                                     >  2:1204042/1‑76 (MQ=255)
                                    acTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgtg                                    >  1:1085697/1‑76 (MQ=255)
                                      tAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                  >  1:955624/1‑76 (MQ=255)
                                      tAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                  >  2:608827/1‑76 (MQ=255)
                                       aaGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGa                                 >  1:1105963/1‑76 (MQ=255)
                                       aaGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGa                                 <  2:225871/76‑1 (MQ=255)
                                            tataAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaa                            >  2:752645/1‑76 (MQ=255)
                                              taAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaa                          <  2:241853/76‑1 (MQ=255)
                                               aaaTGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                         >  1:1456576/1‑76 (MQ=255)
                                               aaaTGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                         >  1:883475/1‑76 (MQ=255)
                                               aaaTGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                         >  1:499773/1‑76 (MQ=255)
                                               aaaTGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                         >  1:1742979/1‑76 (MQ=255)
                                                     tataGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGt                   >  1:553692/1‑76 (MQ=255)
                                                        aGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCt                >  1:642105/1‑76 (MQ=255)
                                                            gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            >  1:294130/1‑76 (MQ=255)
                                                            gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            <  2:1912440/76‑1 (MQ=255)
                                                              tataGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGc          >  1:164500/1‑76 (MQ=255)
                                                              tataGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGc          >  1:1256550/1‑76 (MQ=255)
                                                               ataGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCg         >  1:1121195/1‑76 (MQ=255)
                                                                   ttAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTa     >  1:890886/1‑76 (MQ=255)
                                                                    tAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    >  1:1014856/1‑76 (MQ=255)
                                                                    tAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    >  2:291689/1‑76 (MQ=255)
                                                                    tAGCACCTTTTTATTATTAAATCGTATTAGTCACCAGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    >  2:229327/1‑76 (MQ=255)
                                                                     aGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    <  1:4555/76‑1 (MQ=255)
                                                                     aGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAg   <  2:1493694/76‑1 (MQ=255)
                                                                      gCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGt  <  2:1189439/76‑1 (MQ=255)
                                                                            |                                                                       
TTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGT  >  REL606/1332068‑1332214

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