Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,090 T→C 12.5% intergenic (+63/‑225) PP_0367 → / → PP_0368 hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947445,0900TC12.5% 79.5 / 3.3 32intergenic (+63/‑225)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (13/15);  new base C (2/2);  total (15/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 = 7.25e-01

GGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACG  >  NC_002947/444944‑445237
                                                                                                                                                   |                                                                                                                                                    
ggCAAGACCGATGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                    >  5:84191/1‑149 (MQ=255)
ggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGAGTCATCCCCTATTTCAGCGCCTGCATTGTAGGAGCGGGCTTGCCCGCGAATTCGTcc                                                                                                                                                    >  5:88091/1‑149 (MQ=255)
ggCAAGACCGAAGTCGAGCTCGGAGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTc                                                                                                                                                     <  6:191792/148‑1 (MQ=255)
     gACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                               <  1:71259/149‑1 (MQ=255)
     gACCGAAGTCCACCTAGGCCCAAAGTTATTCAAGCCATCC‑GACGAAATGCTTGTCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                               <  6:235531/149‑1 (MQ=255)
          aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCTGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                <  4:73919/143‑1 (MQ=255)
          aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                >  3:73919/1‑143 (MQ=255)
             tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGacaaaat                                                                                                                                       >  5:214908/1‑146 (MQ=38)
                 gCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                   <  5:72009/149‑1 (MQ=255)
                  cTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTg                                                                                                                                  <  1:237921/149‑1 (MQ=255)
                   tCGGCGCAAAGGTCGTCAAGGCATCCAGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTg                                                                                                                                  >  1:140209/1‑149 (MQ=255)
                     ggCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                 <  6:112385/147‑1 (MQ=255)
                     ggCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  7:249476/149‑1 (MQ=255)
                      gcgcAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                              >  3:189963/1‑149 (MQ=255)
                                       gCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                             <  4:189963/149‑1 (MQ=255)
                                       gCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                             <  6:214908/149‑1 (MQ=18)
                                         aTCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                           <  2:140209/149‑1 (MQ=255)
                                                    aggCTTGGCCGGATCATCCACACCTACGCCTTAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTGGTAGGAGCGGGTTTGCCCTCTAATGCGTCCGATCAGACAAAATTTTTGTCTGGTGGGGCGTATTCGGGGGAAAGCCCGCTc                                                                                                 <  6:88091/147‑1 (MQ=255)
                                                         tGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGccc                                                                                                     >  2:58636/1‑140 (MQ=14)
                                                         tGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGccc                                                                                                     <  1:58636/140‑1 (MQ=14)
                                                                          ccTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCGCGAAAGCGTCCGAACAGACAACAAATTTGTCAGGTCGGGCCACTTAGCGGGGAACCCCCCTc                                                                                                 >  3:62328/1‑127 (MQ=255)
                                                                                  tGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATAAGACAACAATTTTGTCTGGTCGGGCCTCTTAGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTTTGGCATGGGATGACTGa                                                                  >  1:170622/1‑149 (MQ=255)
                                                                                   gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               >  1:150464/1‑88 (MQ=255)
                                                                                   gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  2:150464/88‑1 (MQ=255)
                                                                                    aaTACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTATGGTGGGGCCTCTTAGCGGGGAAGCCCGCTCCG‑AAAAGGGTTTGTGGCATGGCATTACTGAtt                                                                >  8:73320/1‑149 (MQ=39)
                                                                                          tGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAAATGGTCCGAAAACACAAAAAATTTTTATGGGGGGGGCGCTTAGCGGGGAAGCCCCCCACTAACAAGGTTTTTTGG‑ATGTCATGACATATTTTCCAt                                                          >  4:72329/1‑149 (MQ=255)
                                                                                                    ccTATTTCAGCGCCTGCCTTGTTGGCGCGGGCTTGCCCGCTAATGCGTCCGATCCCACAACATTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                <  2:27098/149‑1 (MQ=38)
                                                                                                                                    ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  4:193116/39‑1 (MQ=255)
                                                                                                                                    ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               >  3:193116/1‑39 (MQ=255)
                                                                                                                                       cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGTGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc             >  4:215236/1‑149 (MQ=255)
                                                                                                                                       cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc             >  6:133673/1‑149 (MQ=255)
                                                                                                                                          gcgAATGCGTCCTATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa          <  2:170622/149‑1 (MQ=255)
                                                                                                                                                   tCCGATCAGACAACAATATTGTCTCGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTCATTAGCCAATTTTCGTGGAGTAAATGGCACCATACTCCATCCATGACAGAAAGCTCAGCTAAAca  >  4:172979/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCC‑GACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCT‑ACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACG  >  NC_002947/444944‑445237

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