Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,339,489 Δ1 bp intergenic (+31/‑55) fis → / → yhdJ DNA‑binding protein Fis/predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,339,4880G.87.5% 168.1 / 24.5 49intergenic (+30/‑56)fis/yhdJDNA‑binding protein Fis/predicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (12/30);  minor base A (4/2);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 8.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGA  >  REL606/3339414‑3339560
                                                                          |                                                                        
cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAgc                                                                         >  2:297744/1‑75 (MQ=255)
  tGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAgcgc                                                                       >  1:780583/1‑73 (MQ=255)
  tGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAgcgc                                                                       >  1:782322/1‑73 (MQ=255)
    gTACGCTGCGTAAAAAATTGAAAAAATNCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCt                                                                     >  1:1529674/1‑76 (MQ=255)
    gTACGCTGCGTAAAAAATCGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTa                                                                    <  1:270215/76‑1 (MQ=255)
        gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCg                                                                <  1:942522/76‑1 (MQ=255)
          tGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGc                                                              <  1:606528/76‑1 (MQ=255)
             gTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATg                                                           <  1:733410/76‑1 (MQ=255)
              tAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATgg                                                          <  1:1166566/76‑1 (MQ=255)
                  aaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGaa                                                      <  2:1078003/76‑1 (MQ=255)
                  aaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGaa                                                      <  2:611966/76‑1 (MQ=255)
                  aaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGaa                                                      <  1:683553/76‑1 (MQ=255)
                  aaaGTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGNCATGGGGaa                                                      <  1:36347/76‑1 (MQ=255)
                    aTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAgc                                                    >  2:512072/1‑76 (MQ=255)
                        aaaaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCt                                                <  1:442913/76‑1 (MQ=255)
                         aaaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCtt                                               <  2:583581/76‑1 (MQ=255)
                          aaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCttt                                              >  2:853027/1‑76 (MQ=255)
                          aaaaTACGGCATGAACTAATTCAGGTTAGCGAAATGCTTGATTAACAAGCGCTACTCGGCATNGGGAAGCGCCttt                                              >  2:1392656/1‑76 (MQ=255)
                            aaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCttttt                                            >  2:993374/1‑76 (MQ=255)
                            aaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCttttt                                            >  1:209595/1‑76 (MQ=255)
                               aCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTtat                                         <  1:167464/76‑1 (MQ=255)
                                   cATGAAGTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGgt                                     <  1:1514842/76‑1 (MQ=255)
                                   cATGAACTAATTCAGGTTAGCGAAATGCTTGATTAAAAAGCGCTACTCGGCATTGGGAAGCGCCTTTTTTATAGgt                                     <  2:224897/76‑1 (MQ=255)
                                     tGAACTTATTTAGGTTAGCTATGTGCTTGATTAAAAAGCGCTACGCGGCATGGGGAAGCGCCTCCTTTATAGgtgt                                   >  2:914764/1‑76 (MQ=255)
                                     tGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGgtgt                                   <  1:1229439/76‑1 (MQ=255)
                                         cTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATANGTGTcaca                               >  1:1350410/1‑76 (MQ=255)
                                         cTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTcaca                               <  1:255138/76‑1 (MQ=255)
                                         cTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTcaca                               >  1:1481453/1‑76 (MQ=255)
                                          tAATTCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACaa                              <  2:1530816/76‑1 (MQ=255)
                                              tCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                          <  1:1235278/76‑1 (MQ=255)
                                              tCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                          <  2:1100684/76‑1 (MQ=255)
                                              tCAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAggg                          <  2:48703/76‑1 (MQ=255)
                                               cAGGTTAGCTAAATGCTTGATTAAACAGNGCTACTCGGCCTCGGGAATCGCCTTTTTTATAGGTGTCACAAAGGGa                         <  1:136544/76‑1 (MQ=255)
                                               cAGGTTAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGa                         <  2:331225/76‑1 (MQ=255)
                                                   ttAGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGa                     >  1:702485/1‑76 (MQ=255)
                                                    tAGCTAAATGCTTGATTAAAAAGCGCTACTCGGNATGGGGAAGCGCCTTTTTTATANGTGTCACATANGGAGTGAc                    >  2:1279578/1‑76 (MQ=255)
                                                     aGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGAcc                   <  2:1008814/76‑1 (MQ=255)
                                                     aGCTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGAcc                   <  2:185838/76‑1 (MQ=255)
                                                       cTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCAt                 <  1:1060753/76‑1 (MQ=255)
                                                       cTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCAt                 <  2:134815/76‑1 (MQ=255)
                                                       cTAAATGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCAt                 <  1:864307/76‑1 (MQ=255)
                                                           aTGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATgaga             <  2:1215913/76‑1 (MQ=255)
                                                            tGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGaa            >  1:469346/1‑76 (MQ=255)
                                                            tGCTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGAGAGTGACCATGAGaa            >  2:51476/1‑76 (MQ=255)
                                                              cTTGATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACa          <  2:1164537/76‑1 (MQ=255)
                                                                 gATTAAAAAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGa       <  1:382095/76‑1 (MQ=255)
                                                                    ttaaaaAGCGCTGCTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtg     <  2:1152953/74‑1 (MQ=255)
                                                                    ttaaaaAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtg     <  2:1152970/74‑1 (MQ=255)
                                                                     taaaaAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAGGGGAGTGACCATGAGAACAGGAtgt    >  1:458548/2‑76 (MQ=255)
                                                                      aaaaAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtgtg   >  2:544943/1‑76 (MQ=255)
                                                                      aaaaAGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtgtg   >  2:698887/1‑76 (MQ=255)
                                                                      aaaaAGCGCTACTCGGCATGAGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAtgtg   >  2:480422/1‑76 (MQ=255)
                                                                       aaaaGCGCTACTCGGCATGGGGAAGCGCCTTTNTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATNTGa  >  2:1434817/1‑76 (MQ=255)
                                                                          |                                                                        
CGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGA  >  REL606/3339414‑3339560

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