Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 961,069 G→T intergenic (‑83/‑325) katA ← / → ssuB vegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome961,0690GT96.9% 95.3 / ‑3.3 32intergenic (‑83/‑325)katA/ssuBvegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (15/16);  total (15/17)
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.17e-01

TCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGTAAG  >  BS168_reference_genome/960942‑961197
                                                                                                                                |                                                                                                                                
tcCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTa                                                                                                                         >  2:328649/1‑136 (MQ=255)
  gaaCAGGAGCTCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAAGATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       <  8:112776/135‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATctct                                                                                                                          <  3:256902/133‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       <  1:207790/136‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       <  7:130253/136‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       <  1:61592/136‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAt                                                                                                                        <  8:16927/135‑1 (MQ=255)
  cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAt                                                                                                                        >  4:75114/1‑135 (MQ=255)
    aCCGGAGCGCCCCAGCTAGTTGTCAGTTTAATTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAtt                                                                                                                       <  6:111195/136‑1 (MQ=255)
           cgcCCCAGCTAGTTGTCAGTTTATTTTAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       >  2:378943/1‑127 (MQ=255)
           cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTAtt                                                                                                                       <  1:378943/127‑1 (MQ=255)
                      gTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGa                                                                                                    <  2:143786/135‑1 (MQ=255)
                       ttGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGaaa                                                                                                  <  1:328649/136‑1 (MQ=255)
                          tCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAAttg                                                                                               >  2:193588/1‑136 (MQ=255)
                          tCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCAAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAAttg                                                                                               >  5:389596/1‑136 (MQ=255)
                            aGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAAttgtt                                                                                             <  5:250976/136‑1 (MQ=255)
                                   tttaaCTCATGTTATCACCTCTTGGAATTTATAATAACATAATATAATATCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATATCTATTAATAGATATATATTATAGAAATTGTTTaaaaaa                                                                                      >  6:20890/4‑135 (MQ=255)
                                                          ggAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGa                                                               >  1:31062/1‑136 (MQ=255)
                                                                tataAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACtt                                                          >  4:101689/1‑135 (MQ=255)
                                                                tataAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACtt                                                          >  6:264733/1‑135 (MQ=255)
                                                                                   ctctAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAATTTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa                                      >  8:160305/1‑136 (MQ=255)
                                                                                   ctctAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATAGATAGTATCGAAATTGTTGAAAAACATTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa                                      >  4:43926/1‑136 (MQ=255)
                                                                                     ctAAGT‑AAAAAGTCAATATTATTTTATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAtt                                    >  5:31177/1‑136 (MQ=255)
                                                                                                           ttttATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGaa                <  5:264733/135‑1 (MQ=255)
                                                                                                           ttttATAATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGaa                <  3:101689/135‑1 (MQ=255)
                                                                                                               attattattA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACaa            <  6:389596/132‑1 (MQ=255)
                                                                                                               ataATTATTA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAg           >  6:126138/1‑136 (MQ=255)
                                                                                                                  attattA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCag        <  1:193588/136‑1 (MQ=255)
                                                                                                                  attattA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCag        <  6:31177/136‑1 (MQ=255)
                                                                                                                  attattA‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCag        <  5:20890/136‑1 (MQ=255)
                                                                                                                   ttatta‑AAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCa         >  2:239930/1‑134 (MQ=255)
                                                                                                                         aaaaaaaTTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGTAAg  >  7:246458/2‑136 (MQ=255)
                                                                                                                                |                                                                                                                                
TCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGT‑AAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGTAAG  >  BS168_reference_genome/960942‑961197

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