Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 961,032:1 +A 100% intergenic (‑46/‑362) 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,0321.A100.0% 55.6 / NA 18intergenic (‑46/‑362)katA/ssuBvegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (7/11);  total (7/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAT  >  BS168_reference_genome/960901‑961162
                                                                                                                                    |                                                                                                                                  
gTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTaaaaaag                                                                                                                             <  2:91027/139‑3 (MQ=255)
    gTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCaa                                                                                                                          >  4:109125/1‑139 (MQ=255)
              aCCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTta                                                                                                                >  7:166308/1‑139 (MQ=255)
                  gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                                            <  3:109125/139‑1 (MQ=255)
                    aGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAatta                                                                                                          <  7:126099/139‑1 (MQ=255)
                         ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATtata                                                                                                     <  3:234768/139‑1 (MQ=255)
                                            aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  8:166308/139‑1 (MQ=255)
                                            aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  4:39334/139‑1 (MQ=255)
                                            aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  <  5:256166/139‑1 (MQ=255)
                                             aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                                                  >  7:81338/1‑138 (MQ=255)
                                             aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  3:149053/139‑1 (MQ=255)
                                             aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                                                 <  5:139472/139‑1 (MQ=255)
                                                   gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGata                                                                            >  8:111714/1‑138 (MQ=255)
                                                                              gAACTAATGTTATAAACTCTTGGAATTTATAAGAAAATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGaaatagat                                                 <  8:81338/138‑6 (MQ=255)
                                                                                     tGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAATTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTa                                          >  6:5937/1‑138 (MQ=255)
                                                                                          tCACCTCTTTGAATTTATAAGAACATAATATGATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCt                                    <  1:239208/139‑1 (MQ=255)
                                                                                                  tGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGa                            >  6:247105/1‑139 (MQ=255)
                                                                                                                            ctctAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat  >  2:3127/1‑139 (MQ=255)
                                                                                                                               tAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat  <  5:74730/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAT  >  BS168_reference_genome/960901‑961162

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: