Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,138,494 T→G 39.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
*KT24401,138,4940TG39.0% 3.3 / 5.7 34intergenic (‑108/+33)pigA_1/hemRHeme oxygenase PigA/Hemin receptor
Reads supporting (aligned to +/- strand):  ref base T (7/13);  new base G (7/6);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCG  >  KT2440/1138360‑1138623
                                                                                                                                       |                                                                                                                                     
ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGGCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTg                                                                                                                                      <  7:298540/136‑1 (MQ=255)
 gAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTg                                                                                                                                      >  4:380209/1‑135 (MQ=255)
                 ggtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTAGTAG‑CCGGGGTGGCCGCCAA‑GA‑GGCCgg                                                                                                                    >  6:88706/1‑135 (MQ=255)
                 ggtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCGCTTAG‑CCGGGGGGCACGCCAA‑GA‑GGGCgc                                                                                                                    >  2:113140/1‑134 (MQ=255)
                 ggtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCACGTGG‑CGGGGGGGGCCCCCAA‑GA‑GGCCgg                                                                                                                    >  3:217433/1‑135 (MQ=255)
                      gtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGGGCCCGCAAA‑GA‑GGCCgggtcag                                                                                                               >  2:300221/1‑130 (MQ=255)
                      gtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCACATAG‑CGGGCGGGCCCCCCAA‑GA‑GGCCGGGACAg                                                                                                               >  2:96318/1‑135 (MQ=255)
                       tcctGGGGGAGTCCTGGGAAAGCGGGCGACTAGAAACACGCCGAAACCGCGGGAGCGAGACAGAAAAAACAGCCGGTATTTGCGTTCCCAGCCCCGGCCTCTTCG‑CGGGCGCGCCCGCAAA‑GA‑GGCCGGTACAGa                                                                                                              <  2:320412/132‑1 (MQ=255)
                       gcATGGGGGCGTCCTGGGAACGAGGGCGCCTAGATACACCACGAACCAGCGGCAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTGGCCAGACCCGGCCGCTTCG‑CGGGCGGGCCCGCAAA‑GA‑GGCCGGTACAGa                                                                                                              <  2:285451/134‑1 (MQ=255)
                            ggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGAAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGGGCCCGAAAA‑GA‑GGCCGGTACAGACGTaa                                                                                                         <  4:406934/135‑1 (MQ=255)
                                 gTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTc                                                                                                   <  8:152078/136‑1 (MQ=255)
                                 gTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTc                                                                                                   <  7:34329/136‑1 (MQ=255)
                                    cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCa                                                                                                  >  4:104672/1‑134 (MQ=255)
                                    cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGGGCCCGCCAA‑TA‑GGCCGGTTCAGACGGAACGccg                                                                                                    >  8:35357/1‑131 (MQ=255)
                                      tGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCGTCG‑CGGGCGTGCCCTCAAA‑GC‑GGCCGGTACCGGCGTTACGCCTCAGat                                                                                               >  6:208770/1‑135 (MQ=255)
                                        gAAAGAGGGCGTCTAGATACACGACGAAACAGCGGGCGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGat                                                                                               <  5:88706/133‑1 (MQ=255)
                                            gcgggCGTCTAGAAACAACACGAAACCGCGGGCGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACGGGCCTCTTCG‑CGGGCGGGCCCGAAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCa                                                                                        <  4:217433/133‑1 (MQ=255)
                                                  gTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTAGCCGGGGGGGCCCGCAAG‑AG‑GGCCGGT‑CCGCCGGTACGCCGCAGATATGCCCAAAccc                                                                                   >  8:182576/1‑134 (MQ=255)
                                                      aGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCGCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGGTCAGCCGTACCGCCCGAGAGAT‑CCCAAACCACATTg                                                                              >  2:76818/1‑135 (MQ=255)
                                                             aGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGaa                                                                     <  7:348632/136‑1 (MQ=255)
                                                                 gAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTAG‑CGGGCGTGGCCGCAAA‑GA‑GGCCGGTTCAGACGTAACGCCTCAGCTAT‑CCCCAACCCCATTGAT‑ACCGAAGTTg                                                                 >  8:109469/1‑136 (MQ=255)
                                                                   aaCCGCGGGAGCCAGAACGTAAAAAACGCGGGTTTTGGCGTTGCACTTCCGGGCCGCTCGG‑CGGGGGTGCCGGAAAT‑CA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCTCATTGAT‑ACCGAAGTTgc                                                                <  2:380852/135‑1 (MQ=255)
                                                                      aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGGGCCCGCAAA‑TA‑GGCCCGGGCAGCCGGTACGCCTCAGCTTT‑CCCAAACCACAATGAA‑TCCGAAGGTGCGcc                                                             >  4:319328/1‑135 (MQ=255)
                                                                                  aCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCa                                                 <  3:380209/135‑1 (MQ=255)
                                                                                     gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGc                                              <  2:300795/135‑1 (MQ=255)
                                                                                         aaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCGCATAG‑CGGGGGTGGCCGCCAA‑GACGGCCGG‑ACACACCGTACGCCGCAGAGAT‑CCCCAACCCCAATTATAACCAAGGTTGGCCCCCGGC‑TGGTGAGGCggggg                                         >  8:314864/1‑133 (MQ=255)
                                                                                              gCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                     <  3:104672/135‑1 (MQ=255)
                                                                                                   tATTTGCGTTGCCAGTCCCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                               <  8:24936/136‑1 (MQ=255)
                                                                                                   tATTTGCGTTGCCAGTCCCGGCATCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggct                                <  4:124652/135‑1 (MQ=255)
                                                                                                              cctgTACCGGCCTCTTCC‑GGGGGGTGCCCGAAAA‑GA‑GGCCTGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGTTGTTGGCTGGCTGGGTca                      <  4:358094/131‑1 (MQ=255)
                                                                                                                 gcaCCGGCCTCTGCG‑CGGGCGGCCCCGAAAA‑CA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCGGGCTGGGTCACAgc                  <  5:71480/133‑1 (MQ=255)
                                                                                                                  gaCCGGCCCTGTCC‑CGGGGGGGCCCGGAAACGC‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgcc                 <  7:417017/135‑1 (MQ=255)
                                                                                                                          tctTCG‑CGGGCGTGCCCGAAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGTTGGGTCACAGCCGCCTCAcc         <  2:58351/135‑1 (MQ=255)
                                                                                                                          tctACG‑CGGGCGGGCCCGAAAC‑CA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCAAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAcc         <  6:238000/135‑1 (MQ=255)
                                                                                                                               g‑cgggcggGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCt       <  4:349501/132‑1 (MQ=255)
                                                                                                                                  gggcgggcCGGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcg  <  6:207908/135‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCG‑CGGGCGTGCCCGCAAA‑GA‑GGCCGGTACAGACGTAACGCCTCAGATAT‑CCCAAACCACATTGAT‑ACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCG  >  KT2440/1138360‑1138623

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