Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 2,626,199 Δ2 bp coding (75‑76/273 nt) hupB_2 → DNA‑binding protein HU‑beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,626,1990G.97.9% 217.4 / ‑1.5 48coding (75/273 nt)hupB_2DNA‑binding protein HU‑beta
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (22/25);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01
*ACB1432,626,2000G.95.8% 207.8 / 1.8 48coding (76/273 nt)hupB_2DNA‑binding protein HU‑beta
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (21/25);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACAAACTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCTGGACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAACCCGCAAACCGGCAAGGC  >  ACB143/2626090‑2626327
                                                                                                             ||                                                                                                                               
gAAACAAACTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGTACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGgc                                                                                                       <  7:237343/135‑1 (MQ=255)
gAAACAAACTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGgcg                                                                                                      >  8:165440/1‑136 (MQ=255)
  aacaaacTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGcc                                                                                                    >  8:197485/1‑136 (MQ=255)
       acTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGaa                                                                                               <  8:151655/136‑1 (MQ=255)
          cAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGc                                                                                             >  1:62898/1‑135 (MQ=255)
          cAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGc                                                                                             >  4:25581/1‑135 (MQ=255)
               gagaAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTg                                                                                       <  7:197485/136‑1 (MQ=255)
               gagaAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTg                                                                                       <  6:88120/136‑1 (MQ=255)
               gagaAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTg                                                                                       <  6:224378/136‑1 (MQ=255)
                agaAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGa                                                                                      <  3:272057/136‑1 (MQ=255)
                   aaTACGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACg                                                                                    <  8:57308/135‑1 (MQ=255)
                   aaTAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGa                                                                                   >  6:42688/1‑136 (MQ=255)
                     tAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCCGCGCCCTGAAGCAAGGTGACGAtg                                                                                 >  1:163851/1‑136 (MQ=255)
                        ggggACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACCTCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACTCAGGCATCGCATAAGTCACCGGCGCCCGGAAGCAAGGTGACGATGTgg                                                                              >  1:135682/1‑136 (MQ=255)
                             cTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTg                                                                         <  2:170179/136‑1 (MQ=255)
                               tAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTggg                                                                       >  3:42307/1‑135 (MQ=255)
                                agagTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGgc                                                                                                       >  8:288527/1‑103 (MQ=255)
                                agagTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGgc                                                                                                       <  7:288527/103‑1 (MQ=255)
                                  agTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTgg                                                                    <  4:347568/136‑1 (MQ=255)
                                      aaCAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTg                                                                                       >  8:238108/1‑113 (MQ=255)
                                      aaCAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTg                                                                                       <  7:238108/113‑1 (MQ=255)
                                      aaCAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCttt                                                                <  4:148173/136‑1 (MQ=255)
                                           gTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTAc                                                           <  8:295918/136‑1 (MQ=255)
                                           gTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTAc                                                           <  2:206306/136‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  2:339548/135‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  2:181849/135‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  7:163980/135‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  2:6578/135‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  3:323483/135‑1 (MQ=255)
                                              ggAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCt                                                         <  4:180658/135‑1 (MQ=255)
                                               gAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAATCCAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTtc                                                       >  5:240863/1‑136 (MQ=255)
                                               gAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGC‑‑GACGCAGTCATCGAATACCTCACCGGCGCCCAGAAGCAAGGTGCAGATGTGGTCCTGGTTGGCGTTGGGACCTtc                                                       >  6:221062/1‑136 (MQ=255)
                                                   tGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCCTCGAATCCCTCACCGGCCCCCTGAAGCAAGCTGACGATGTGGTACTGGTTGGCGTTGGTACCTTCACgc                                                   >  3:158691/1‑135 (MQ=255)
                                                   tGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCgg                                                   >  8:276263/1‑136 (MQ=255)
                                                        gACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCa                                                                                                            <  4:316795/74‑1 (MQ=255)
                                                        gACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCa                                                                                                            >  3:316795/1‑74 (MQ=255)
                                                        gACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAg                                              >  5:94455/1‑136 (MQ=255)
                                                        gACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAg                                              >  7:50738/1‑136 (MQ=255)
                                                        gACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAg                                              >  8:133385/1‑136 (MQ=255)
                                                         aCGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAg                                              <  5:10223/135‑1 (MQ=255)
                                                               tCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAgcgcg                                       <  7:170588/136‑1 (MQ=255)
                                                                  ccgcATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCg                                    >  1:77317/1‑136 (MQ=255)
                                                                  ccgcATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCg                                    >  3:140145/1‑136 (MQ=255)
                                                                             gACATCCCGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGc                         <  1:200601/136‑1 (MQ=255)
                                                                                   ccGAAAGCTGTAGCCGGCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAAc                   >  6:112466/1‑136 (MQ=255)
                                                                                             tAGCGGTCCGTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGTCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAACCCGCAAACCg         <  6:191901/136‑1 (MQ=255)
                                                                                                ccggccgTGCGCT‑‑ACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTAGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAACCCGCAAACCggc       >  7:86749/1‑135 (MQ=255)
                                                                                                       tGCGCTGGACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAACCCGCAAACCggcaaggc  >  3:96668/1‑135 (MQ=255)
                                                                                                             ||                                                                                                                               
GAAACAAACTCAATAGAGAAATAAGGGGACTTAGAGTGAACAAGTCGGAACTGATTGACGCTATCGCCGCATCTGCAGACATCCCGAAAGCTGTAGCCGGCCGTGCGCTGGACGCAGTCATCGAATCCGTCACCGGCGCCCTGAAGCAAGGTGACGATGTGGTACTGGTTGGCTTTGGTACCTTCTCGGTCAAGGAGCGCGCCGAGCGTACTGGCCGCAACCCGCAAACCGGCAAGGC  >  ACB143/2626090‑2626327

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