Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,860,407 A→G 15.1% R33R (CGT→CGC)  gntR_2 ← HTH‑type transcriptional regulator GntR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,860,4070AG15.1% 80.4 / 6.9 33R33R (CGT→CGC) gntR_2HTH‑type transcriptional regulator GntR
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (2/3);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGA  >  ACB143/3860273‑3860540
                                                                                                                                      |                                                                                                                                     
gacgGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGCCCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCAc                                                                                                                                      >  8:62537/1‑136 (MQ=255)
                   ccGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGt                                                                                                                   <  2:278027/136‑1 (MQ=255)
                      aCTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTaa                                                                                                                 >  1:114630/1‑135 (MQ=255)
                                 gCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCGCGCAAGGCGCGCGAAGCGGTAATCGGGGAAAccc                                                                                                     >  6:272008/1‑136 (MQ=255)
                                 gCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCGCGCAAGGCGCGCGAAGCGGTAATCGGGGAAAccc                                                                                                     >  2:92867/1‑136 (MQ=255)
                                               ccGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGc                                                                                        >  6:48029/1‑135 (MQ=255)
                                               ccGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCa                                                                                       <  8:77340/136‑1 (MQ=255)
                                                cGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCaa                                                                                      <  8:256863/136‑1 (MQ=255)
                                                    cggcCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTc                                                                                  <  6:132234/136‑1 (MQ=255)
                                                    cggcCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTc                                                                                  <  4:126944/136‑1 (MQ=255)
                                                    cggcCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTc                                                                                  <  1:259511/136‑1 (MQ=255)
                                                              ggCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGt                                                                         <  2:208126/135‑1 (MQ=255)
                                                                     tAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTATTCGGGTAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGgacgac                                                                 <  2:114630/136‑1 (MQ=255)
                                                                     tAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGgacgac                                                                 <  1:46056/136‑1 (MQ=255)
                                                                     tAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGgacgac                                                                 <  1:41434/136‑1 (MQ=255)
                                                                              ttCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCGACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTc                                                        >  1:286584/1‑136 (MQ=255)
                                                                               tCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCg                                                       >  5:183768/1‑136 (MQ=255)
                                                                                 gccgcTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCgg                                                      >  6:74436/1‑135 (MQ=255)
                                                                                   cgcTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAg                                                   >  5:45106/1‑136 (MQ=255)
                                                                                     cTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGcc                                                 >  7:54339/1‑136 (MQ=255)
                                                                                     cTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGcc                                                 >  4:159198/1‑136 (MQ=255)
                                                                                      tGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCGCGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCg                                                <  5:117660/136‑1 (MQ=255)
                                                                                      tGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGGGCTGACCCCGCGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCg                                                <  5:20052/136‑1 (MQ=255)
                                                                                        cggcgATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCg                                                <  5:55529/134‑1 (MQ=255)
                                                                                          gcgATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGc                                            >  3:221946/1‑136 (MQ=255)
                                                                                                        caccaGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAAtgtt                              >  3:138410/1‑136 (MQ=255)
                                                                                                        caccaGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAAtgt                               <  3:78205/135‑1 (MQ=255)
                                                                                                               tCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAAttt                       <  8:87720/136‑1 (MQ=255)
                                                                                                               tCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAAttt                       >  4:43179/1‑136 (MQ=255)
                                                                                                               tCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAAttt                       >  4:224790/1‑136 (MQ=255)
                                                                                                                     gcAACGGTGCTGACCCCGCGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTtactac                 <  1:92867/136‑1 (MQ=255)
                                                                                                                                    ccACGCAAGGGGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCGTGGCAACTTGGGCCTGTGTGTGCGGACCTGTAGTTCGGGAGCCGGTGCGGGACATTAATGTTGTAATTTTACTACTTGCAAGGAGAGGGa  <  4:43211/136‑1 (MQ=255)
                                                                                                                                    ccACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGa  <  5:154453/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGA  >  ACB143/3860273‑3860540

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: ‑