Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,332,148 (T)6→7 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,1421.T97.3% 134.2 / ‑3.0 37intergenic (+125/‑85)topA/cysBDNA topoisomerase I/DNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (19/17);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAA  >  REL606/1332067‑1332212
                                                                             |                                                                      
ttttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCAcc                                                                         >  1:1506734/1‑76 (MQ=255)
  ttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtt                                                                      >  2:1182805/1‑76 (MQ=255)
  ttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtt                                                                      >  1:1124259/1‑76 (MQ=255)
    ttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttt                                                                    >  2:318751/1‑76 (MQ=255)
      aTTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttttt                                                                  <  1:1460355/76‑1 (MQ=255)
      aTTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttttt                                                                  >  1:169521/1‑76 (MQ=255)
               tattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTaa                                                          <  1:27521/76‑1 (MQ=255)
                 ttttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAAt                                                        <  2:761877/76‑1 (MQ=255)
                     ttcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTa                                                    >  2:1403163/1‑76 (MQ=255)
                      tcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTAt                                                   >  1:901262/1‑76 (MQ=255)
                         tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAg                                                <  1:119319/76‑1 (MQ=255)
                         tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAg                                                >  2:871377/1‑76 (MQ=255)
                          cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGt                                               >  1:351421/1‑76 (MQ=255)
                            cacCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCa                                             <  2:345394/76‑1 (MQ=255)
                                tataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCg                                         >  1:844258/1‑76 (MQ=255)
                                 ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGc                                        >  1:1464399/1‑76 (MQ=255)
                                 ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGc                                        >  2:839410/1‑76 (MQ=255)
                                 ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGc                                        >  2:695687/1‑76 (MQ=255)
                                 ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGc                                        >  1:1464411/1‑76 (MQ=255)
                                         aGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                >  2:202065/1‑76 (MQ=255)
                                         aGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                >  2:189695/1‑76 (MQ=255)
                                         aGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                >  2:202075/1‑76 (MQ=255)
                                          ggCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGa                               >  1:896062/1‑76 (MQ=255)
                                          ggCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGa                               <  2:805075/76‑1 (MQ=255)
                                             tataAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGa                            >  1:282211/1‑76 (MQ=255)
                                              ataAATGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaa                           <  1:1503708/76‑1 (MQ=255)
                                                  aTGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                       <  2:341801/76‑1 (MQ=255)
                                                  aTGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaaa                       <  2:600864/76‑1 (MQ=255)
                                                   tGATATAGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAAc                      >  1:603858/1‑76 (MQ=255)
                                                         aGTGTTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGCCCCCCGCCAGGTGTGATGAAAAAACGATGtt                <  1:602024/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCGCGTTTTTTTATTATTAATCCGTATTAGCCGCCCGCCAGGTGTGACGAAAAAACGATGtt                <  2:586399/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCCCTTTTTTTTATTATTAAATCGTATTAGTCCCCCGCCAGGTGTGACGAAAAAACGATGtt                <  1:996905/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACTTTTTTTTATTATTAAATCGTATTAGTCCCCCGCCAGGTGTGACGAAAAAACGATGtt                <  1:383765/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCCCCCGCCAGGTATGACGAAAAAACGATGtt                <  2:675446/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGtt                <  2:974605/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGtt                <  2:1051904/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGtt                <  2:1159436/76‑1 (MQ=255)
                                                         aGTGGTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGtt                <  2:1048289/76‑1 (MQ=255)
                                                             gTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGa            <  2:1185065/76‑1 (MQ=255)
                                                             gTTATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGa            >  1:694054/1‑76 (MQ=255)
                                                              ttATAGTTAGCACCTTTTTTTATTATTAAATCGTATTAGTCCCCCGCCAGGTGTGACGAAAAAACGATGTTCTGAt           <  1:1188010/76‑1 (MQ=255)
                                                                     tAGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCt    >  2:299495/1‑76 (MQ=255)
                                                                      aGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTa   <  2:282211/76‑1 (MQ=255)
                                                                      aGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTa   <  2:1061718/76‑1 (MQ=255)
                                                                      aGCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTa   >  1:477024/1‑76 (MQ=255)
                                                                       gCACCTTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa  >  1:7409/1‑76 (MQ=255)
                                                                             |                                                                      
TTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAA  >  REL606/1332067‑1332212

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