Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,083,515 A→T 8.2% I72N (ATT→AAT)  hmoB ← heme‑degrading monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,083,5150AT 8.2% 135.6 / 3.7 49I72N (ATT→AAT) hmoBheme‑degrading monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (27/18);  new base T (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GGTCAGAATGACATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCATCAAAA  >  BS168_reference_genome/1083379‑1083653
                                                                                                                                        |                                                                                                                                          
ggTCAGAATGACATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAata                                                                                                                                          <  7:206062/139‑1 (MQ=255)
ggTCAGAATGACATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAata                                                                                                                                          <  5:324158/139‑1 (MQ=255)
ggTCAGAATGACATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAata                                                                                                                                          <  2:333159/139‑1 (MQ=255)
           cATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATAc                                                                                                                               >  7:364547/1‑139 (MQ=255)
           cATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTCCCTGAGTAACGGCAATATTGTTCAATa                                                                                                                                >  4:70435/1‑138 (MQ=255)
            aTAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTTCTTCAAAGCCAGGCTCGTTTTCGACTTTTACCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACt                                                                                                                              >  2:320265/1‑139 (MQ=255)
               ggTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTCCCTGAGTAACGGCAATATTGTTCAATACTGCa                                                                                                                           <  3:70435/139‑1 (MQ=255)
                   tCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAAc                                                                                                                       >  8:185142/1‑139 (MQ=255)
                   tCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACATTTCCCGCTCTGTATTTAAATCTGTTTTAAAACATCGGACGCCCTTCCTGAGTAACGGCAATAATGTTCAATACAGCAAAAc                                                                                                                       >  6:99900/1‑139 (MQ=255)
                     aCTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTATTACCGCTCTGTTTTTTAATATTTTTTCAAACAGCGGAAGCCCTTCCTGAGTAAAGGCAATATTGTTAAATACAGCAAAAACg                                                                                                                     >  5:56616/1‑139 (MQ=255)
                     aCTGTCAAGCGGACGCAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCg                                                                                                                     >  7:361571/1‑139 (MQ=255)
                      cTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCgg                                                                                                                    <  8:364547/139‑1 (MQ=255)
                         tCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCggg                                                                                                                   <  5:121614/137‑1 (MQ=255)
                           aaGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGt                                                                                                               >  8:227983/1‑139 (MQ=255)
                               ggacggaGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGa                                                                                                           <  2:172784/139‑1 (MQ=255)
                                gacggaGCTCTCGAATTGCTTCAAAGCCTGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTTAATGTGTTTTAAAACAGCGGACGCCCTTCCTGAGTAACGGCAATTTTGTTCAATACTGCAAAACCGGGATGTTTGAt                                                                                                          <  6:56616/139‑1 (MQ=255)
                                gacgcaGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGAt                                                                                                          <  8:361571/139‑1 (MQ=255)
                                 acggaGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGAtc                                                                                                         >  7:15144/1‑139 (MQ=255)
                                   ggaGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACGTTTCCCGCTCTGTTTTTAAATCTGTTTTCAACCAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGAtctc                                                                                                       >  8:158582/1‑139 (MQ=255)
                                      gCACTCGAATTGCTTCAAAGCCAGGCGCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTcc                                                                                                    <  6:98762/139‑1 (MQ=255)
                                            gAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACt                                                                                                <  2:61370/137‑1 (MQ=255)
                                            gAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACt                                                                                                <  7:227983/137‑1 (MQ=255)
                                                  cTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAAt                                                                                        <  7:186748/139‑1 (MQ=255)
                                                     cAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGa                                                                                      <  3:63588/138‑1 (MQ=255)
                                                            aGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCAt                                                                                >  2:243027/1‑137 (MQ=255)
                                                            aGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCAt                                                                                >  3:89164/1‑137 (MQ=255)
                                                             ggCTCGTTTTCGCCTTTTCCCGCTCTGTTTTTAAATCTGTTTTGAAACATCGGAAGCCCTTACTGAGTAACGGCCATATTGTTCAATAAGTCAAAACCGGGATGTTTGTTCTCTGCGACTTGATAAATGAAGTCATAAg                                                                             >  4:308848/1‑139 (MQ=255)
                                                                          cTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGAt                                                                                                          <  1:275606/97‑1 (MQ=255)
                                                                          cTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGAt                                                                                                          >  2:275606/1‑97 (MQ=255)
                                                                                   agctGTTTTTATATCTGTGTTCATACAGGGGAGGGCCTTCCTGAGTAAAGTCAATATTGTTCAATACTGCAAAAACGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTGATAAGCGTGTGGAGCCTGGAAAACGGt                                                       <  7:271687/137‑1 (MQ=255)
                                                                                                tCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGCCAATATTGTTCAATACTACAAAACCGGGAAGTTTGATCTCTCCGACTAGGTCAATGACTTCATAACCGTGAGGTGCCTGGAAAACGGTATCTCCgtttgct                                          >  8:271687/1‑137 (MQ=255)
                                                                                                  tGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCa                                                                                 <  6:66336/98‑1 (MQ=255)
                                                                                                  tGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCa                                                                                 >  5:66336/1‑98 (MQ=255)
                                                                                                  tGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTACGACTAGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTCTCTCCgtgtgttt                                         >  3:85605/1‑138 (MQ=255)
                                                                                                              gCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCa                                                                                                                           <  5:153426/44‑1 (MQ=255)
                                                                                                              gCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCa                                                                                                                           >  6:153426/1‑44 (MQ=255)
                                                                                                                         tCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTGATGAATTAAAATGGCGGTTTCCt                 >  1:127234/1‑139 (MQ=255)
                                                                                                                         tCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCt                 >  5:169878/1‑139 (MQ=255)
                                                                                                                              aGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAAt                                                                                        >  6:244513/1‑63 (MQ=255)
                                                                                                                              aGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAAt                                                                                        <  5:244513/63‑1 (MQ=255)
                                                                                                                              aGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGcccc            >  5:211429/1‑139 (MQ=255)
                                                                                                                              aGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGcccc            >  2:56497/1‑139 (MQ=255)
                                                                                                                                tAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAAt                                                                                        <  1:104113/61‑1 (MQ=255)
                                                                                                                                tAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAAt                                                                                        >  2:104113/1‑61 (MQ=255)
                                                                                                                                  aCGGCATTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTaaaa                             >  7:266767/1‑118 (MQ=255)
                                                                                                                                  aCGGCATTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTaaaa                             >  7:273232/1‑118 (MQ=255)
                                                                                                                                  aCGGCATTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTaaaa                             <  8:266767/118‑1 (MQ=255)
                                                                                                                                  aCGGCATTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTaaaa                             <  8:273232/118‑1 (MQ=255)
                                                                                                                                  aCGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCa        >  8:25692/1‑139 (MQ=255)
                                                                                                                                   cGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCAt       >  5:220037/1‑139 (MQ=255)
                                                                                                                                   cGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCAt       >  2:7886/1‑139 (MQ=255)
                                                                                                                                      cAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCATCaa    <  7:158582/139‑1 (MQ=255)
                                                                                                                                       aaTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTTCATCaaa   <  4:89164/139‑1 (MQ=255)
                                                                                                                                       aaTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCATCaaa   <  1:31166/139‑1 (MQ=255)
                                                                                                                                       aaTATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCATCaaa   <  1:56497/139‑1 (MQ=255)
                                                                                                                                        atatTGTTCAATACTGCAAAACCGGGATGTTAGATCTCTCCGACTTGGTCAATGACTTCATCAGCGTGTGGTGCCTGGAAAACGGTATCTCCGGTTGTTTCATGAAGTATCCTGGCGTTTTCCTGCCCCGGCAGCaagc  >  6:351682/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                          
GGTCAGAATGACATAGGTGTCACTGTCAAGCGGACGGAGCACTCGAATTGCTTCAAAGCCAGGCTCGTTTTCGACTTTTCCCGCTCTGTTTTTAAATCTGTTTTCAAACAGCGGACGCCCTTCCTGAGTAACGGCAATATTGTTCAATACTGCAAAACCGGGATGTTTGATCTCTCCGACTTGGTCAATGACTTCATAAGCGTGTGGTGCCTGGAAAACGGTATCTCCGTTTGTTTCATGAATTAAAATGGCGTTTTCCTGCCCCTGCATCAAAA  >  BS168_reference_genome/1083379‑1083653

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