Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,339,322 T→C V54A (GTA→GCA)  fis → DNA‑binding protein Fis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,339,3220TC97.8% 178.3 / ‑4.1 46V54A (GTA→GCA) fisDNA‑binding protein Fis
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/21);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TCGGTTAAACAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGTACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCT  >  REL606/3339249‑3339392
                                                                         |                                                                      
tCGGTTAAACAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCAc                                                                      >  2:1015258/1‑76 (MQ=255)
 cGGTTAAACAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACt                                                                     <  2:333894/76‑1 (MQ=255)
  ggTTAAACAGGCACCGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTg                                                                    <  2:675334/76‑1 (MQ=255)
    ttAAACAGGCACTGAAGAACTATTTTTCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGc                                                                  >  1:455632/1‑76 (MQ=255)
     tAAACAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCt                                                                 >  1:782563/1‑76 (MQ=255)
       aaCAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCActgactga                                                               >  1:220022/1‑76 (MQ=255)
          aGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGaaat                                                            >  1:342907/1‑74 (MQ=255)
            gCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAg                                                          <  1:1025532/76‑1 (MQ=255)
             naCNGNAGTACGATNTTTCTGACGTGAATGGTNAGGATGTGGATGACCTCTATGAGCTGACGCTGGCTGAAGTAGa                                                         <  1:915223/75‑1 (MQ=255)
             ctcTGAAGAACTATTTTGCTCAACTGAATCGTCAGGATGTCAATGACCGCGATGATCTGGCACTGGCTGAAGTAGa                                                         <  1:117866/74‑1 (MQ=255)
              aCTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGTACTGGTTGAAGTAGaa                                                        <  1:521122/76‑1 (MQ=255)
               cTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAAc                                                       <  2:606926/76‑1 (MQ=255)
                tGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACa                                                      <  1:355377/76‑1 (MQ=255)
                tGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACa                                                      <  1:342926/76‑1 (MQ=255)
                tGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACa                                                      <  2:820253/76‑1 (MQ=255)
                  aagaaCTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGc                                                    >  2:612290/1‑76 (MQ=255)
                      aCTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCt                                                <  2:650106/76‑1 (MQ=255)
                       cTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTg                                               <  2:790222/76‑1 (MQ=255)
                        tATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGt                                              >  1:255582/1‑76 (MQ=255)
                           tttGCTCAACTGAATGGTCAGGATGTGAATGCCCTCTATGAGCT‑GCACTGGCTGAAGTAGAACAGCCCCTGTTGGa                                          <  2:396897/76‑1 (MQ=255)
                           tttGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGAAGAACAGCCCCTGTTgg                                           >  2:148495/1‑76 (MQ=255)
                           tttGCTCAACTGAATGGTCAGGATGTGAATGACATCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCACTGTTgg                                           <  1:713516/76‑1 (MQ=255)
                           tttGCTCAACTGAATGGTCAGGATGTCAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTgg                                           >  1:956147/1‑76 (MQ=255)
                               cTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACAt                                       <  2:823244/76‑1 (MQ=255)
                               cTCAACTGAAGTGTCAGGCTGTGAATGACCTCTAGGGGCGGGCACTGGCTGAACGGGAACAGGCCCAGTTGGACAt                                       <  1:106703/76‑1 (MQ=255)
                                tCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATg                                      >  2:634537/1‑76 (MQ=255)
                                    cTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGa                                  >  2:667318/1‑76 (MQ=255)
                                    cTGAATGGTCAGGATGTGAATGACATCTATGAGCTGGCACTGGCTGACGTAGAACAGCCCCTGTTGTAGATGGGGa                                  >  2:667314/1‑76 (MQ=255)
                                    cTGAATGGTCAGGATGTGAATCACGCATATGCGCTGGCACCGACTGAAGTAGTACAGGCCCTGCTCGACATGGTgc                                  >  1:992797/1‑75 (MQ=255)
                                     tGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGAt                                 >  2:625250/1‑76 (MQ=255)
                                     tGAATGGCCAGGACGTGGATGACCTCTCTGAGCTGGCACTGGCCGAAGTAGAACAGCCCCCGTTGGACATGGTGAt                                 >  1:965869/1‑76 (MQ=255)
                                          ggTCAGGATGTGAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAAt                            >  2:510313/1‑76 (MQ=255)
                                                   gtgAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTg                   <  2:45492/76‑1 (MQ=255)
                                                    tgAATGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTgg                  >  2:546771/1‑76 (MQ=255)
                                                      aaTGACCTCTATGAGCTGGCACTGGCTAAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTa                <  2:450576/76‑1 (MQ=255)
                                                       aTGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTaa               <  2:380459/76‑1 (MQ=255)
                                                       aTGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTaa               <  1:718948/76‑1 (MQ=255)
                                                        tGACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAAc              <  2:756435/76‑1 (MQ=255)
                                                         gACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAAcc             <  2:313567/76‑1 (MQ=255)
                                                         gACCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAAcc             >  1:308212/1‑76 (MQ=255)
                                                          aCCTCTATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGATAACCa            >  1:307640/1‑76 (MQ=255)
                                                             tctATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGAAACCAGAc         >  1:260347/1‑76 (MQ=255)
                                                               tATGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGAccc       >  1:15655/1‑76 (MQ=255)
                                                                 tGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCgc     >  2:826589/1‑75 (MQ=255)
                                                                 tGAGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGt     >  2:826600/1‑76 (MQ=255)
                                                                   aGCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTAGTAACCAGACCCGtgc   >  2:442030/1‑76 (MQ=255)
                                                                    gCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGtgct  <  2:554188/76‑1 (MQ=255)
                                                                    gCTGGCACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGtgct  >  1:878811/1‑76 (MQ=255)
                                                                         |                                                                      
TCGGTTAAACAGGCACTGAAGAACTATTTTGCTCAACTGAATGGTCAGGATGTGAATGACCTCTATGAGCTGGTACTGGCTGAAGTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCT  >  REL606/3339249‑3339392

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