Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 961,062 Δ1 bp 100% intergenic (‑76/‑332) 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,0620T.96.9% 110.1 / ‑1.7 32intergenic (‑76/‑332)katA/ssuBvegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (16/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGT  >  BS168_reference_genome/960935‑961194
                                                                                                                                |                                                                                                                                    
gATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTaaaaaaatttat                                                                                                                            <  3:206483/138‑7 (MQ=255)
gATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTaaaaaaatttat                                                                                                                            <  5:101912/138‑7 (MQ=255)
      ctcCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTAt                                                                                                                     >  3:150471/1‑138 (MQ=255)
         cAAGCGGAGCGCACCAGCTAGCTGTCAGTTTAGTAGAACTCATGTTATCACCTGTAGGAATTTATAAGAACAGAATATCATCTCTCAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAACTTTATCTCTATTaa                                                                                                                  <  4:263354/138‑1 (MQ=255)
         cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTAtt                                                                                                                    <  7:51002/136‑1 (MQ=255)
         cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCACCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAAt                                                                                                                 <  6:166778/139‑1 (MQ=255)
         cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGATATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAAt                                                                                                                 <  7:378456/139‑1 (MQ=255)
           aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATc                                                                                                                >  6:178655/1‑138 (MQ=255)
           aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCg                                                                                                               <  8:17286/139‑1 (MQ=255)
                aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTaa                                                                                                                  >  1:131761/1‑131 (MQ=255)
                aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTaa                                                                                                                  <  2:131761/131‑1 (MQ=255)
                aGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGatata                                                                                                          <  4:374120/139‑1 (MQ=255)
                                        aTTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAAc                                                                                   <  4:150471/138‑1 (MQ=255)
                                         tttGATCTCATGTTATCACCTCTTTGAATTTATAAGCACATAATATCATCTCTAAGTAAAAAGTGAATAATATTTTATAATTATTAAATAAATTTATCTCAATTAATCGATATATAGTATCGAAATTGTTGAAAAACCt                                                                                 <  3:367426/139‑1 (MQ=255)
                                               cTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCTATATATAGTATCGAATTTGATGAAAAACAtttttg                                                                            >  2:164984/1‑138 (MQ=255)
                                               cTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCttgttg                                                                            >  3:65125/1‑138 (MQ=255)
                                                     ttATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAg                                                                      <  8:134168/138‑1 (MQ=255)
                                                            ctctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTAt                                                              >  4:3418/1‑139 (MQ=255)
                                                              ctTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGa                                                            >  4:297153/1‑139 (MQ=255)
                                                                       tataAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTtaat                                                    <  3:3418/138‑1 (MQ=255)
                                                                                       catcTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa                                   >  7:95094/1‑139 (MQ=255)
                                                                                       catcTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa                                   >  1:384961/1‑139 (MQ=255)
                                                                                       catcTCTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataa                                   >  6:93878/1‑139 (MQ=255)
                                                                                         tctcTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAtt                                 >  1:249400/1‑139 (MQ=255)
                                                                                         tctcTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAtt                                 >  8:222295/1‑139 (MQ=255)
                                                                                         tctcTAAGTAAAAAGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAAtt                                 >  2:357035/1‑139 (MQ=255)
                                                                                          ctctAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAaataat                                  >  5:83253/1‑139 (MQ=255)
                                                                                               aGTAAAAAGTAAATATTAATTTATAATTATTAAAAAAATTTATCTCTATTAATCTATATATAGTATCGAAATTGTTGAAAAACATAGATGTAACTGCTTTTTATGAGACTAAATAATTTTTGGAGAAATAATATTTCCt                           >  8:355138/1‑139 (MQ=255)
                                                                                                     aaaGTCAATATTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTtctc                       >  6:148725/1‑138 (MQ=255)
                                                                                                            atatTATTTTATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACATTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGTAGAAATAATTTTTCCTTCTCTATTTGc                >  8:20840/1‑138 (MQ=255)
                                                                                                                    ttATAATTATTAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAg        <  2:249400/138‑1 (MQ=255)
                                                                                                                         attattAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGt  <  5:178655/139‑1 (MQ=255)
                                                                                                                         attattAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGt  <  5:148725/139‑1 (MQ=255)
                                                                                                                         attattAAAAAAATTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGt  <  8:95094/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                    
GATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTTAATAATTTTTGGAGAAATAATTTTTCCTTCTCTATTTGCGAAACAAGCAGAGT  >  BS168_reference_genome/960935‑961194

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: