Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 924,066:1 +A 100% intergenic (‑46/‑233) U712_04490 ← / → U712_04495 Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881924,0661.A82.4% 48.0 / 3.3 17intergenic (‑46/‑233)U712_04490/U712_04495Vegetative catalase/Aliphatic sulfonates import ATP‑binding protein SsuB
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (1/13);  total (3/14)
Fisher's exact test for biased strand distribution p-value = 6.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTT  >  CP006881/923939‑924175
                                                                                                                                 |                                                                                                              
gTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCGCTTGGAATTTATAAGAACATAATTTCATCTCTAAGTaaaaaagtcac                                                                                                    >  5:35452/1‑133 (MQ=255)
  ccGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAata                                                                                                   <  7:83553/139‑1 (MQ=255)
              gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                       <  4:135819/139‑1 (MQ=255)
              gCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAat                                                                                       <  8:133346/139‑1 (MQ=255)
                     ttGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATtata                                                                                <  1:39571/139‑1 (MQ=255)
                               aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                      <  3:115809/139‑1 (MQ=255)
                               aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                                                                      <  5:127626/139‑1 (MQ=255)
                                         aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                            <  2:283818/139‑1 (MQ=255)
                                         aCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt                                                            <  6:253875/139‑1 (MQ=255)
                                           cGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCg                                                          <  3:143698/139‑1 (MQ=255)
                                               gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAGAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                             >  5:116099/1‑132 (MQ=255)
                                               gcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAACAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTaa                                                             <  6:116099/132‑1 (MQ=255)
                                                cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata                                                     <  5:218372/139‑1 (MQ=255)
                                                   cccAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAg                                                  <  4:233246/139‑1 (MQ=255)
                                                           gTTGTAAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGa                                            <  1:127289/137‑1 (MQ=255)
                                                            ttGTCAGTTTATTTGAACTCATGTTATCACATCTTGGAATTTATAAGAACATAATATAAACTATAATAAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAATAATCTATATATAGTATCTAAAt                                         >  8:252635/1‑139 (MQ=255)
                                                                   tttATTTGAACTCATGTTATCTCCTCTTGTAATTTATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATAGAAATTGTTGaa                                  <  7:252635/139‑1 (MQ=255)
                                                                                                    ttATAAGAACATAATATCATCTCTAAGTAAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACtt  >  8:77410/1‑138 (MQ=255)
                                                                                                                                 |                                                                                                              
GTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATATAGTATCGAAATTGTTGAAAAACCTTGTTGTAACAGCTTTTTATGACACTT  >  CP006881/923939‑924175

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: