Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 433,653 T→G F50V (TTT→GTT)  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,6530TG95.6% 184.0 / ‑1.3 45F50V (TTT→GTT) hupBHU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (25/18);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCT  >  REL606/433582‑433724
                                                                       |                                                                       
aTGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGc                                                                     <  1:1434828/76‑1 (MQ=255)
  gCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCg                                                                   <  2:58608/76‑1 (MQ=255)
     attattGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTa                                                                <  2:878571/76‑1 (MQ=255)
      ttattGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTaa                                                               <  1:444298/76‑1 (MQ=255)
       tattGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGCTaaa                                                              >  2:1017581/1‑76 (MQ=255)
          tGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAagag                                                           >  2:418414/1‑76 (MQ=255)
           gCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGc                                                          >  2:126093/1‑76 (MQ=255)
           gCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGc                                                          >  1:1024828/1‑76 (MQ=255)
               ccGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGtgc                                                      >  1:335545/1‑76 (MQ=255)
               ccGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGtgc                                                      <  1:639745/76‑1 (MQ=255)
                 gTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGtgctg                                                    >  2:690502/1‑76 (MQ=255)
                 gTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGtgctg                                                    >  2:726764/1‑76 (MQ=255)
                  tAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGtgctgc                                                   >  2:701186/1‑76 (MQ=255)
                   aaCTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGcc                                                  >  2:1436876/1‑76 (MQ=255)
                    aCTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTATAGAGCGTGCTGcca                                                 <  1:846090/76‑2 (MQ=255)
                     cTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCg                                                <  2:276333/76‑1 (MQ=255)
                      tGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGt                                               <  2:397689/76‑1 (MQ=255)
                        aaTCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTAc                                             <  1:25712/76‑1 (MQ=255)
                          tctcTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTg                                           >  2:597322/1‑76 (MQ=255)
                               gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGc                                      >  1:144736/1‑76 (MQ=255)
                               gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGc                                      >  2:1425863/1‑76 (MQ=255)
                                aaaGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCa                                     >  2:983174/1‑76 (MQ=255)
                                aaaGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCa                                     <  2:1026904/76‑1 (MQ=255)
                                aaaGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCa                                     <  1:1403126/76‑1 (MQ=255)
                                 aagaagGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCaa                                    >  2:1267365/1‑76 (MQ=255)
                                 aagaagGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCaa                                    >  2:609809/1‑76 (MQ=255)
                                 aagaagGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCaa                                    <  2:402884/76‑1 (MQ=255)
                                  agaagGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAAc                                   <  1:39784/76‑1 (MQ=255)
                                     agGGGATGATGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCg                                <  2:1331825/76‑1 (MQ=255)
                                         gatgatGTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGa                            >  1:345962/1‑76 (MQ=255)
                                          atgatgTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGAc                           <  2:1178671/76‑1 (MQ=255)
                                            gatgTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCg                         <  2:890112/76‑1 (MQ=255)
                                             atgTAGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGCCCgg                        >  2:1233261/1‑76 (MQ=255)
                                                 aGCACTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaaa                    >  2:585594/1‑76 (MQ=255)
                                                     cTGGTAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAgaga                <  2:250583/76‑1 (MQ=255)
                                                         tAGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc            >  1:44434/1‑76 (MQ=255)
                                                          aGGTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAcc           <  1:678071/76‑1 (MQ=255)
                                                           ggTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCa          >  1:264838/1‑76 (MQ=255)
                                                           ggTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCa          <  1:960481/76‑1 (MQ=255)
                                                            gTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCAt         >  1:800624/1‑76 (MQ=255)
                                                            gTTTTGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCAt         <  1:85653/76‑1 (MQ=255)
                                                              tttGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCg       >  2:1361870/1‑76 (MQ=255)
                                                                tGGTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCgct     >  2:891419/1‑76 (MQ=255)
                                                                 ggTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCgctg    >  2:83818/1‑76 (MQ=255)
                                                                  gTACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCgctgc   >  1:904872/1‑76 (MQ=255)
                                                                   tACTGTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCgctgct  >  1:1331575/1‑76 (MQ=255)
                                                                       |                                                                       
ATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCT  >  REL606/433582‑433724

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