Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,586,035 +G intergenic (+118/+97) 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,0331.G97.1% 108.1 / ‑3.0 34intergenic (+116/+99)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (13/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

GAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAG‑GGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAA  >  NC_002947/4585897‑4586110
                                                                                                                                           |                                                                               
gAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGc                                                                      <  4:103633/149‑1 (MQ=255)
  gCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTc                                                                    <  2:36550/149‑1 (MQ=255)
           gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATTACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAg‑g                                                           >  3:227536/1‑148 (MQ=39)
            aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGa                                                         >  2:96197/1‑149 (MQ=255)
             aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGA‑c                                                       <  6:130049/149‑1 (MQ=255)
                cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCTCG‑GGTAAACCCGCTCCTACCAG‑GGA‑ccgt                                                    <  4:227536/149‑4 (MQ=39)
                 ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGA‑ccgtg                                                   >  8:9597/1‑145 (MQ=39)
                   aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGcc                                                  >  5:397941/1‑149 (MQ=39)
                   aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGA‑ccgtgc                                                  <  6:392321/148‑6 (MQ=39)
                    aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCg                                                 <  3:119332/149‑1 (MQ=39)
                              aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATg                                       <  2:340909/149‑1 (MQ=39)
                                          tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGcccc                           <  7:9597/149‑1 (MQ=39)
                                           gAACAGCTGATCTGCGCTGCTCCATTGAAAGGCTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                          >  3:386500/1‑149 (MQ=38)
                                            aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTg                         >  8:72959/1‑149 (MQ=39)
                                             aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGATCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg                        <  5:394608/149‑1 (MQ=39)
                                             aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg                        <  1:96197/149‑1 (MQ=39)
                                             aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTg                         <  3:106355/148‑1 (MQ=39)
                                              cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                       <  6:397941/149‑1 (MQ=39)
                                                       tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                       >  1:48210/1‑141 (MQ=38)
                                                       tGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                       <  2:48210/140‑1 (MQ=38)
                                                         cgcTGCTCCATTGAAAGGCTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg               <  4:386500/146‑1 (MQ=38)
                                                          gctgctCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt           <  7:72959/149‑1 (MQ=38)
                                                              ctCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc         <  1:60507/147‑1 (MQ=38)
                                                                ccATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc         >  7:363823/1‑145 (MQ=38)
                                                                ccATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc         <  8:363823/145‑1 (MQ=38)
                                                                 cATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcg      <  7:41479/147‑1 (MQ=38)
                                                                 cATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCGTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcg      >  8:41479/1‑147 (MQ=38)
                                                                 cATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATAGATCGGGCCGGCCTC‑TTCGCGGGTTAAAACCCTCCTAGCAA‑GGACCGGGCAGGCCATGAAATCGGCGCGGCCCCAGGGAAGAGCGGGGTTAACagcg    >  1:303172/1‑149 (MQ=255)
                                                                   tttaaaGGTTGTGGATGACC‑GGGTGCTTGGCAAGGTGAACATCTTGCATTGATCGGCCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCCG‑GGACCGTGCCGGCCTTTAATGCGGCGCTGCCCCTGGTAGGTGCGGTTTTACCCGCGaa  <  4:123960/146‑1 (MQ=37)
                                                                   ttGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg               >  2:36351/1‑136 (MQ=38)
                                                                   ttGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg               <  1:36351/136‑1 (MQ=38)
                                                                   ttGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTA‑TTAGCGGGTAAAACCGCAACAACCAGTGCC‑CGGGCCGGCCTTGCATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTTACCGCGaa  >  8:300616/1‑149 (MQ=255)
                                                                            tGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg               >  2:230373/1‑127 (MQ=38)
                                                                            tGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAG‑GGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg               <  1:230373/127‑1 (MQ=38)
                                                                                                                                           |                                                                               
GAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAG‑GGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAA  >  NC_002947/4585897‑4586110

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: ‑