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.94.6% 328.3 / 12.9 74intergenic (+126/‑84)topA/cysBDNA topoisomerase I/DNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base . (45/25);  total (49/25)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGT  >  REL606/1332069‑1332214
                                                                          |                                                                       
ttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtt                                                                        <  2:645208/76‑1 (MQ=255)
ttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtt                                                                        >  2:1315622/1‑76 (MQ=255)
ttttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtt                                                                        <  1:1446306/76‑1 (MQ=255)
 tttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttt                                                                       <  1:1521820/76‑1 (MQ=255)
 tttATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttt                                                                       >  2:1449557/1‑76 (MQ=255)
   tATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttttt                                                                     >  1:1666040/1‑76 (MQ=255)
   tATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCttttt                                                                     >  2:1469583/1‑76 (MQ=255)
     ttGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttttat                                                                   >  2:1006164/1‑74 (MQ=255)
     ttGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttttat                                                                   >  2:1752326/1‑74 (MQ=255)
      tGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCtttttatt                                                                  <  2:1390191/76‑4 (MQ=255)
         atatTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTttattatt                                                              <  1:1370942/76‑1 (MQ=255)
           attattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTaa                                                            >  1:1771238/1‑76 (MQ=255)
           attattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTaa                                                            >  2:1487567/1‑76 (MQ=255)
           attattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTaa                                                            >  1:1025053/1‑76 (MQ=255)
            ttattTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTaa                                                            >  2:442382/1‑76 (MQ=255)
              attTTTTCTTCACAACTATACACTAAGGCTCTAAATGATATAGTGGTTATAGTTAGCACATTTTTATTATTACATc                                                         >  2:1775172/1‑76 (MQ=255)
               ttttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCg                                                        <  2:1136296/76‑1 (MQ=255)
                tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       >  2:997812/1‑76 (MQ=255)
                tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       >  1:1632079/1‑76 (MQ=255)
                tttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGt                                                       >  1:1431690/1‑76 (MQ=255)
                 ttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGt                                                       >  1:486893/1‑76 (MQ=255)
                 ttttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTa                                                      >  1:883707/1‑76 (MQ=255)
                  tttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAt                                                     >  2:1888757/1‑76 (MQ=255)
                  tttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAt                                                     <  1:995730/76‑1 (MQ=255)
                  tttCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAt                                                     >  2:627511/1‑76 (MQ=255)
                   ttcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAtt                                                    >  1:1178313/1‑76 (MQ=255)
                   ttcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTAtt                                                    >  1:1010465/1‑76 (MQ=255)
                    tcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTa                                                   >  2:436944/1‑76 (MQ=255)
                    tcttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTa                                                   >  1:336861/1‑76 (MQ=255)
                     cttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAg                                                  >  1:461079/1‑76 (MQ=255)
                     cttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAg                                                  <  1:1194007/76‑1 (MQ=255)
                      ttcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGt                                                 >  2:728090/1‑76 (MQ=255)
                       tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTc                                                <  1:1528206/76‑1 (MQ=255)
                       tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTc                                                >  2:1819877/1‑76 (MQ=255)
                       tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTc                                                <  1:113169/76‑1 (MQ=255)
                       tcACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTc                                                <  1:1727603/76‑1 (MQ=255)
                        cacaCCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCa                                               >  2:1523985/1‑76 (MQ=255)
                           acCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCAcca                                            <  2:1696406/76‑2 (MQ=255)
                            ccTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCg                                           <  2:58848/76‑1 (MQ=255)
                             cTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGc                                          <  1:1766224/76‑1 (MQ=255)
                             cTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGc                                          <  2:1138094/76‑1 (MQ=255)
                             cTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGc                                          <  2:1427490/76‑1 (MQ=255)
                              tataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGcc                                         >  2:191119/1‑76 (MQ=255)
                               ataCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCa                                        <  2:702821/76‑1 (MQ=255)
                                tcnacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAg                                       >  2:1331120/4‑76 (MQ=255)
                                taCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAg                                       >  1:897437/1‑76 (MQ=255)
                                taCACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAg                                       <  1:553890/76‑1 (MQ=255)
                                  cacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgt                                     <  1:954191/76‑1 (MQ=255)
                                  cacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgt                                     >  2:1790433/1‑76 (MQ=255)
                                  cacTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgt                                     <  1:1580625/76‑1 (MQ=255)
                                    cTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGgtgt                                   >  1:1097923/1‑76 (MQ=255)
                                     tAAGGCTATAAATGATGTAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                  <  2:1052561/76‑1 (MQ=255)
                                     tAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGtgtg                                  >  2:1668432/1‑76 (MQ=255)
                                         gCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACg                               >  2:1715718/1‑76 (MQ=255)
                                          cTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaa                             <  1:965483/76‑1 (MQ=255)
                                          cTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaa                             >  1:264025/1‑76 (MQ=255)
                                          cTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaa                             <  2:1752406/76‑1 (MQ=255)
                                             taAATGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGaaaaa                          <  2:1761453/76‑1 (MQ=255)
                                                 tGATATAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGa                      >  2:764628/1‑76 (MQ=255)
                                                   atatAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATg                    >  2:1189365/1‑76 (MQ=255)
                                                   atatAGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATg                    >  2:1657271/1‑76 (MQ=255)
                                                    tataGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGt                   >  1:828698/1‑76 (MQ=255)
                                                     attgTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAAGTGTGACGAAAAAACGATGt                   >  2:1531885/4‑76 (MQ=255)
                                                      taGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTc                 >  1:548073/1‑76 (MQ=255)
                                                       aGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCt                >  2:1531123/1‑76 (MQ=255)
                                                       aGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCt                >  1:755689/1‑76 (MQ=255)
                                                       aGTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCt                >  2:1009090/1‑76 (MQ=255)
                                                        gTGGTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTg               <  1:1429915/76‑1 (MQ=255)
                                                          gtttATAGTTACCCCCTTTTTATTATTAAATCGTATTAGTCCCCCGCCAGGTGTGACGAAAAAACGATGTTCTGAt             <  2:755689/74‑1 (MQ=255)
                                                          ggTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGAt             >  2:1207267/1‑76 (MQ=255)
                                                           gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            <  2:1023673/76‑1 (MQ=255)
                                                           gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            >  2:519601/1‑76 (MQ=255)
                                                           gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            >  2:686504/1‑76 (MQ=255)
                                                           gTTATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATg            >  2:1252584/1‑76 (MQ=255)
                                                            ttATAGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATgg           >  2:1607115/1‑76 (MQ=255)
                                                             tataGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGc          >  1:1762116/1‑76 (MQ=255)
                                                               taGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGt        >  2:51595/1‑76 (MQ=255)
                                                               taGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGt        <  1:1861947/76‑1 (MQ=255)
                                                                aGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTc       <  1:38261/76‑1 (MQ=255)
                                                                aGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTc       >  2:36113/1‑76 (MQ=255)
                                                                aGTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTc       >  2:1531805/1‑76 (MQ=255)
                                                                 gTTAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCt      >  1:807571/1‑76 (MQ=255)
                                                                   tAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    >  1:1582262/1‑76 (MQ=255)
                                                                   tAGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTaa    >  2:1266494/1‑76 (MQ=255)
                                                                    aGCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAg   >  1:1731319/1‑76 (MQ=255)
                                                                     gCACCTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGt  <  2:1151199/76‑1 (MQ=255)
                                                                          |                                                                       
TTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAGTGGTTATAGTTAGCACCTTTTTTATTATTAAATCGTATTAGTCACCCGCCAGGTGTGACGAAAAAACGATGTTCTGATGGCGTCTAAGT  >  REL606/1332069‑1332214

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