Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,586,056:1 +C 100% intergenic (+139/+76) PP_4061 → / ← PP_4063 hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,586,0561.C91.7% 66.8 / 0.3 24intergenic (+139/+76)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (10/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC  >  NC_002947/4585910‑4586195
                                                                                                                                                      |                                                                                                                                              
tggCGTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑c                                                                                                                                               <  8:66013/148‑1 (MQ=39)
    ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg                                                                                                                                           >  4:117978/1‑145 (MQ=39)
    ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg                                                                                                                                           >  2:67947/1‑145 (MQ=39)
        gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCgg                                                                                                                                       >  1:245817/1‑149 (MQ=39)
              aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTg                                                                                                                                 >  5:152136/1‑148 (MQ=39)
              aCCAAGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTg                                                                                                                                 >  7:237844/1‑148 (MQ=38)
                  cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑c                                                                                                                           >  1:121069/1‑149 (MQ=39)
                    gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg                                                                                                                         <  2:121069/149‑1 (MQ=39)
                    gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg                                                                                                                         <  8:237844/149‑1 (MQ=39)
                                            cgcTGCTCCATTGAAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCAGGGCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                    <  2:247765/146‑1 (MQ=255)
                                            cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGG‑CAGCGCC‑GCA‑TTCAAGGCCGGCACGGTCCCTGGTAGGAGc                                                                                                      <  6:140589/144‑1 (MQ=14)
                                            cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGc                                                                                                      <  2:245817/144‑1 (MQ=38)
                                            cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                    <  1:67947/146‑1 (MQ=38)
                                            cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCg                                                                                                     <  3:117978/145‑1 (MQ=38)
                                             gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGttt                                                                                                <  6:152136/149‑1 (MQ=38)
                                                    cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAg                                                                                                          >  8:193665/1‑132 (MQ=38)
                                                    cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAg                                                                                                          <  7:193665/132‑1 (MQ=38)
                                                         aacggTTGAGGATGAGCGGGTGGTTGGCATGGTGAAGATCTTGCATTGCTCGGGCGGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGGTGCCGCTGGTAGGAGCGGGTTTACCCGCGAaga                                                                                     <  6:280481/145‑1 (MQ=25)
                                                         aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGt                                                                                                            >  1:239291/1‑125 (MQ=38)
                                                         aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGt                                                                                                            <  2:239291/125‑1 (MQ=38)
                                                                      gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCCCC‑GGTAAAACCCCACCAAACAACGAA‑C‑CAGCCGCCCATTAAATA‑C‑CCGCAGCCCCGGGTAGGAGCGGGTTTTACCGCCAAGAGTCCGGCCCCCCCa                                                                       >  6:150769/1‑149 (MQ=255)
                                                                            gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGcccc                                                                                                                >  7:280160/1‑102 (MQ=25)
                                                                            gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGcccc                                                                                                                <  8:280160/102‑1 (MQ=25)
                                                                                        gTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATg‑cgg                                                                                                                         <  6:104310/81‑1 (MQ=21)
                                                                                                                             gggTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCaa                                                              >  3:353800/2‑104 (MQ=255)
                                                                                                                                     ccGCTCC‑TACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTCATCTCACGCATGCGGAACtt         >  6:280758/1‑149 (MQ=255)
                                                                                                                                             tACCAGGGACC‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc  >  2:306997/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCC‑TTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC  >  NC_002947/4585910‑4586195

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