Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 433,771 2 bp→AC coding (266‑267/273 nt) hupB → HU, DNA‑binding transcriptional regulator, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606433,7710TA96.5% 216.5 / ‑0.8 58V89D (GTA→GAC) hupBHU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (30/26);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01
*REL606433,7720AC94.8% 229.5 / ‑0.5 58V89D (GTA→GAChupBHU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (30/25);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCTTTTTTA  >  REL606/433701‑433840
                                                                      ||                                                                    
ggTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACt                                                                  <  2:1067254/76‑1 (MQ=255)
ggTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACt                                                                  <  2:1039187/76‑1 (MQ=255)
 gTAAAGAGATCCCCATCGCTGCTGCTAAAGTATCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTa                                                                 <  1:1377801/76‑1 (MQ=255)
  tAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTaa                                                                <  1:912209/76‑1 (MQ=255)
   aaaGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAg                                                               <  1:799561/76‑1 (MQ=255)
   aaaGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAg                                                               >  2:1589274/1‑76 (MQ=255)
    aaGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGc                                                              >  2:45612/1‑76 (MQ=255)
      gagaTCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGt                                                            <  2:558489/76‑1 (MQ=255)
       agaTCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGtt                                                           <  1:29865/76‑1 (MQ=255)
        gaTCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTg                                                          <  2:197047/76‑1 (MQ=255)
         aTCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGt                                                         >  2:1615434/1‑76 (MQ=255)
         aTCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGt                                                         <  1:7615/76‑1 (MQ=255)
           cACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTcc                                                       >  1:147233/1‑76 (MQ=255)
           cACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTcc                                                       >  1:405627/1‑76 (MQ=255)
              cATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCa                                                    <  2:1003293/76‑1 (MQ=255)
               aTCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAg                                                   >  2:319087/1‑76 (MQ=255)
               aTCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAg                                                   >  1:525578/1‑76 (MQ=255)
                 cGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTg                                                 <  2:962578/76‑1 (MQ=255)
                 cGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTg                                                 >  1:414226/1‑76 (MQ=255)
                 cGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTg                                                 >  1:235703/1‑76 (MQ=255)
                        gctAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGAtgt                                          >  1:1348813/1‑76 (MQ=255)
                        gctAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGAtgt                                          >  1:408530/1‑76 (MQ=255)
                         ctAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGAtgtg                                         >  2:1585379/1‑76 (MQ=255)
                          tAAAGTACCGAGCGTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGa                                        <  2:886581/76‑1 (MQ=255)
                           aaaGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGAc                                       >  1:440317/1‑76 (MQ=255)
                           aaaGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGAc                                       >  1:187731/1‑76 (MQ=255)
                            aaGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCGCAGTGGGGATGTGACg                                      >  1:1535443/1‑76 (MQ=255)
                             aGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGa                                     >  2:319709/1‑76 (MQ=255)
                             aGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGa                                     >  2:561984/1‑76 (MQ=255)
                             aGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGa                                     >  1:952421/1‑76 (MQ=255)
                              gTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGaa                                    <  1:1346311/76‑1 (MQ=255)
                                aCCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGt                                  >  1:227686/1‑76 (MQ=255)
                                aCCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGt                                  >  2:618228/1‑76 (MQ=255)
                                   gAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCa                               <  1:1009529/76‑1 (MQ=255)
                                      cTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAgg                            <  2:959124/76‑1 (MQ=255)
                                         ccGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGgcg                         >  1:210061/1‑76 (MQ=255)
                                         ccGTGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGgcg                         <  1:286152/76‑1 (MQ=255)
                                            tGCAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCAt                      >  2:574369/1‑76 (MQ=255)
                                              cAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCt                    <  1:505359/76‑1 (MQ=255)
                                              cAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCt                    <  2:1328253/76‑1 (MQ=255)
                                              cAGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCt                    <  1:40045/76‑1 (MQ=255)
                                               aGGTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTa                   >  1:1041626/1‑76 (MQ=255)
                                                ggTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTan                  <  1:1071401/76‑2 (MQ=255)
                                                 gTAAAGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACt                 <  1:976310/76‑1 (MQ=255)
                                                   aaaGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGa               <  2:409907/76‑1 (MQ=255)
                                                   aaaGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGa               <  2:1383142/76‑1 (MQ=255)
                                                   aaaGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGa               >  1:924036/1‑76 (MQ=255)
                                                    aaGCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGAt              <  1:1406836/76‑1 (MQ=255)
                                                      gCACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGAtgt            <  2:660396/76‑1 (MQ=255)
                                                       cACTGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGAtgtg           <  1:38266/76‑1 (MQ=255)
                                                          tGAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCt        >  1:1162410/1‑76 (MQ=255)
                                                          tGAAAGACGCGGAAAGCAAAGCGTTGTCCCCAGTGGGGATGAGACGAAGTTCAAGGGCGCATCTACTGATGTGCCt        <  1:830451/76‑1 (MQ=255)
                                                           gAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCtt       >  1:1326892/1‑76 (MQ=255)
                                                           gAAAGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCtt       <  2:140893/76‑1 (MQ=255)
                                                            aaaGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCACGGGCGCATCTACTGATGTGCCttt      >  2:172986/1‑76 (MQ=255)
                                                             aaGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCtttt     >  2:148375/1‑76 (MQ=255)
                                                              aGACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCttttt    >  1:1194173/1‑76 (MQ=255)
                                                               gACGCGGACAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCtttttt   >  1:322759/1‑76 (MQ=255)
                                                                aCGCGGACAACTAAGCGTTGTCGCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCTTTTTTa  <  1:142851/76‑1 (MQ=255)
                                                                      ||                                                                    
GGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAGCGTTGTCCCCAGTGGGGATGTGACGAAGTTCAAGGGCGCATCTACTGATGTGCCTTTTTTA  >  REL606/433701‑433840

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