Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,138,494 T→G 43.0% intergenic (‑108/+33) pigA_1 ← / ← hemR Heme oxygenase PigA/Hemin receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,138,4940TG43.0% 1.3 / 4.1 25intergenic (‑108/+33)pigA_1/hemRHeme oxygenase PigA/Hemin receptor
Reads supporting (aligned to +/- strand):  ref base T (3/11);  new base G (8/3);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCG  >  ACB122/1138367‑1138623
                                                                                                                                |                                                                                                                                 
aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCCCGTAGCGGGCGGGccccc                                                                                                                             >  2:153713/1‑132 (MQ=255)
  ggCTGGGCGGACGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCACATCGCGGGCGGGGCCGCaaa                                                                                                                          >  8:93857/1‑135 (MQ=255)
          gggcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTGGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCgg                                                                                                                  <  7:56961/132‑1 (MQ=255)
                   tGGGGGCGGCCTTGGAAAGAGGCCGGCTAGAAACAAGACGAACCCGCGGCCGCGAGACAGTAAAAAACCCCGGTATTTGCGTTGCCAGACCCGGCCTCTTCGCGGGCGGGCCCGAAACGAGGCCGGTACAGACg                                                                                                          <  4:158564/134‑1 (MQ=255)
                      gggAGTCCTTGGAAAGAGGGCGTCTACATACACGACGAAACAGCGGGAGCGAGACAGTAAAAAACGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGGGCCCGCAAAGAGGCCGGTACAGACGTAAc                                                                                                      <  7:465614/135‑1 (MQ=255)
                         aGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGGGGGCCCCCAAAGAGGCCGGGTCAGGCCTTACGcc                                                                                                   >  7:461182/1‑135 (MQ=255)
                           cccTTGCAACGAGGGCGGCTAGATACAAGACGAAACCGCGGGAGCGAGACAGTAAAAAACGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                 <  5:464921/134‑1 (MQ=255)
                             cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTAGCGCGCGTGCCCCCAAAGAGGCCCGTTCAAACCGAACGcc                                                                                                   >  5:239068/1‑131 (MQ=255)
                                     gcgggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCa                                                                                       <  8:151436/133‑1 (MQ=255)
                                                         cGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGGGCCCGCAAAGAGGCCGGTACAGACCTGACCCCGCCGATATCCCCAAACACATTTATAACGCAGt                                                                   >  3:356488/1‑135 (MQ=255)
                                                               aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCGCGTCGCGGGCGGGCCCGCAAAGAGGCCGGGTCCGCCGGTACGCCGCCGATATCCCACACCCCATTTATACCGCAGGTGCGcc                                                             >  1:30414/1‑135 (MQ=255)
                                                                    gAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGGGCCCGCCAAGAGGCCCGGTCCGACGGAACGCCGCAGATATCCCAAACCCCATTGATACCGAAGGTGCGCCCCGGc                                                        >  3:236621/1‑135 (MQ=255)
                                                                    gAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCGCTTCGCGGGGGGGGCCGCCAAGAGGGCGGGACACACGTAACGCCGCCGCTATCCCCAACCACAATGCTACCGAAGTTGCGcccccg                                                         >  8:239418/1‑132 (MQ=255)
                                                                    gAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCGAGTAGCGGGGGGGGCCGCCAAGAGGCCGGGTCAGCCGTTACGCCGCCGCTATACCCAACCCCCTTTATTCCCAACTTGCGcccccg                                                         >  4:231879/1‑132 (MQ=255)
                                                                          gaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTATCGGGCGTGGCCGCAAAGAGGGCGGTACATACGTAAAGCCTCAGATATCCCAAACCACATTGATACCGCAGTTGCGCCCCCGCAAgg                                                    >  8:236899/1‑133 (MQ=255)
                                                                              gaaaaaaacGCCGGTATTGGCGTGGCCAGTCCCGGCGTCTGCGCGGGCGGGCCCGAAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGc                                              <  3:463586/133‑1 (MQ=255)
                                                                               gtaaaaaaCCCGGGATTTTGCGTTGCCAGACCCGGCATCTTCTCGGGCGTGCCCGAAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGc                                              <  8:122667/133‑1 (MQ=255)
                                                                                            tATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                <  5:182150/135‑1 (MQ=255)
                                                                                            tATTGGCGTTCCCAGTCCCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                <  2:442018/135‑1 (MQ=255)
                                                                                                        ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTca                      <  2:308099/134‑1 (MQ=255)
                                                                                                            tACCGGCCTGTTCGCGGGCGTGCCGGAAAAGATGGCGGTACAGACGTAAGGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGGATGGTCTGGCGGTCGATGTTGGCTGGCTGGGTCACAgcc                 <  7:63776/135‑1 (MQ=255)
                                                                                                                    tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAcc         <  6:239068/135‑1 (MQ=255)
                                                                                                                      tacgcgGGCGTGCCCGCAAAGAGGCGGGTACAGACGTAACGCCTCAGATATGCCAAACCACATTGATACCGATCTTGCGCCCCCGCATGGTCAGGCGGTCGAGGTGGGCTGGCTGGGTCACAGCCGGCTCACCCt       <  4:111921/133‑1 (MQ=255)
                                                                                                                         gcgGGCGTGGCCGCCAAGAGGCCGGTCCAGACGGAACGCCTGAGATATCCCACCCCCCCTTGGTACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTGGATGTTGGCTGGGTGGGTCCCCGCCGCCTACCCCGGGc    >  8:331925/1‑135 (MQ=255)
                                                                                                                           gggCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg  <  8:151102/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                 
AGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCG  >  ACB122/1138367‑1138623

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