Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,925,446 A→T 19.4% intergenic (‑151/+414) PP_2563 ← / ← PP_2564 antibiotic biosynthesis protein/quercetin 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,925,4460AT19.4% 66.0 / 11.3 31intergenic (‑151/+414)PP_2563/PP_2564antibiotic biosynthesis protein/quercetin 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base T (3/3);  total (14/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.70e-01

TCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGAA  >  NC_002947/2925301‑2925595
                                                                                                                                                 |                                                                                                                                                     
tCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCggg                                                                                                                                                   >  2:475616/1‑150 (MQ=255)
    gTGCTCCGCAGGAACTTTGTGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCACAAACTTAATGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAg                                                                                                                                               <  1:195346/150‑1 (MQ=255)
      gCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGca                                                                                                                                             <  1:475616/150‑2 (MQ=255)
        tCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgt                                                                                                                                           >  1:283121/1‑150 (MQ=255)
          cGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAAACCTGCTAACCCGACGGGTTAGCggggt                                                                                                                                         <  2:329125/150‑4 (MQ=255)
          cGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAAACCTGCTAACCCGACGGGTTAGCggggt                                                                                                                                         >  1:329125/1‑147 (MQ=255)
               gAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCa                                                                                                                                    >  2:158703/1‑150 (MQ=255)
               gAACGTTGGGCTGGCTGGTAGCACAACGCTAGGGCATGCTGGCTGATGGCTCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCa                                                                                                                                    >  1:13908/1‑150 (MQ=255)
                    ttGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGga                                                                                                                               <  2:13908/150‑1 (MQ=255)
                                  aGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCt                                                                                                                 >  1:495856/1‑150 (MQ=255)
                                      caACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCg                                                                                                             <  2:247175/150‑1 (MQ=255)
                                               tGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgttg                                                                                                                                         <  1:242693/113‑1 (MQ=255)
                                               tGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgttg                                                                                                                                         >  2:242693/1‑113 (MQ=255)
                                                gCAGGCTGGCTAATCCAGCGAATATTGGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTg                                                                                                   >  2:55529/1‑150 (MQ=255)
                                                 cAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGt                                                                                                  <  2:561828/150‑1 (MQ=255)
                                                   ggctggctGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGcag                                                                                                                                            >  2:495125/1‑104 (MQ=255)
                                                   ggctggctGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGcag                                                                                                                                            <  1:495125/106‑3 (MQ=255)
                                                        gctgATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGcaca                                                                                            <  1:368832/149‑1 (MQ=255)
                                                        gctgATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGcaca                                                                                            >  2:368832/1‑149 (MQ=255)
                                                            aTGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCACATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCaaa                                                                                       <  2:25793/150‑1 (MQ=255)
                                                                 gacAATATTTGAACGCTGCCTTGCGATAGCCGGGCAGGTCCGGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGa                                                                                  >  2:24498/3‑150 (MQ=255)
                                                                     atatTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATgcg                                                                              >  1:157238/1‑150 (MQ=255)
                                                                             aGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGt                                                                      <  2:157238/150‑1 (MQ=255)
                                                                                               cGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATc                                                    <  1:362153/150‑1 (MQ=255)
                                                                                                           gggCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCCCCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGaaa                                        <  1:568107/150‑1 (MQ=255)
                                                                                                                 cAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATgg                                                                                                                                 <  1:476653/55‑1 (MQ=255)
                                                                                                                 cAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATgg                                                                                                                                 >  2:476653/1‑55 (MQ=255)
                                                                                                                  aaGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCCAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTc                                 <  2:398328/150‑1 (MQ=255)
                                                                                                                                      ccgCTAACCCGTCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTgc                                                                                                       <  2:454498/58‑1 (MQ=37)
                                                                                                                                      ccgCTAACCCGTCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTgc                                                                                                       >  1:454498/3‑60 (MQ=37)
                                                                                                                                            aCCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAgccg       <  2:495856/150‑1 (MQ=255)
                                                                                                                                                 aCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGaa  <  1:158703/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGAA  >  NC_002947/2925301‑2925595

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