Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,586,030 2 bp→TC intergenic (+113/+101) 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,0300CT92.9% 37.2 / ‑4.6 14intergenic (+113/+102)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (1/12);  total (2/12)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029474,586,0310TC92.9% 38.7 / ‑5.5 14intergenic (+114/+101)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (1/12);  minor base G (1/0);  total (2/12)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCC  >  NC_002947/4585888‑4586105
                                                                                                                                               ||                                                                            
ccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gg                                                                         <  4:74713/148‑1 (MQ=255)
              aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCt                                                            <  5:176879/148‑1 (MQ=255)
              aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCt                                                            >  7:39071/1‑148 (MQ=255)
                  gTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCa                                                        <  4:163388/148‑1 (MQ=255)
                   tGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCATCGCG‑GGTAAAACCCCGCCTACCac                                                       >  1:32255/1‑147 (MQ=255)
                    gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAgg                                                      <  4:152258/148‑1 (MQ=255)
                      aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGa                                                    <  3:76643/148‑1 (MQ=255)
                                  aGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCt                                        <  5:27615/148‑1 (MQ=39)
                                    aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                       <  1:19535/147‑1 (MQ=39)
                                         ggCACCGGCGTGAACAGCTGAGCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgc                                 <  2:82589/148‑1 (MQ=38)
                                                  gTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGcc                        <  1:127616/148‑1 (MQ=39)
                                                       cAGCTGATCTGCGCTGCTCCACTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCTGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTg                    <  6:136324/147‑1 (MQ=38)
                                                         gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                 <  8:159835/148‑1 (MQ=38)
                                                                        tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAccc  <  8:39071/148‑1 (MQ=38)
                                                                                                                                               ||                                                                            
CCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCC  >  NC_002947/4585888‑4586105

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