Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,095 T→C 14.5% intergenic (+68/‑220) 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,0950TC14.5% 65.4 / 3.8 30intergenic (+68/‑220)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (10/14);  new base C (2/2);  total (14/16)
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.51e-01

GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACG  >  NC_002947/444949‑445237
                                                                                                                                                   |                                                                                                                                               
gACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                               <  1:71259/149‑1 (MQ=255)
gACCGAAGTCCACCTAGGCCCAAAGTTATTCAAGCCATCCGACGAAATGCTTGTCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                               <  6:235531/149‑1 (MQ=255)
     aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCTGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                <  4:73919/143‑1 (MQ=255)
     aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                >  3:73919/1‑143 (MQ=255)
        tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGacaaaat                                                                                                                                       >  5:214908/1‑146 (MQ=38)
            gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                   <  5:72009/149‑1 (MQ=255)
             cTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTg                                                                                                                                  <  1:237921/149‑1 (MQ=255)
              tCGGCGCAAAGGTCGTCAAGGCATCCAGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTg                                                                                                                                  >  1:140209/1‑149 (MQ=255)
                ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                 <  6:112385/147‑1 (MQ=255)
                ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  7:249476/149‑1 (MQ=255)
                 gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                              >  3:189963/1‑149 (MQ=255)
                                  gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                             <  4:189963/149‑1 (MQ=255)
                                  gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                             <  6:214908/149‑1 (MQ=18)
                                    aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                           <  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)
                                                                             tGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATAAGACAACAATTTTGTCTGGTCGGGCCTCTTAGCGGGCAAGCCCGCTCCTACAAGGGTTTTTGGCATGGGATGACTGa                                                                  >  1:170622/1‑149 (MQ=255)
                                                                              gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  2:150464/88‑1 (MQ=255)
                                                                              gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               >  1:150464/1‑88 (MQ=255)
                                                                               aaTACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTATGGTGGGGCCTCTTAGCGGGGAAGCCCGCTCCGAAAAGGGTTTGTGGCATGGCATTACTGAtt                                                                >  8:73320/1‑149 (MQ=39)
                                                                                     tGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAAATGGTCCGAAAACACAAAAAATTTTTATGGGGGGGGCGCTTAGCGGGGAAGCCCCCCACTAACAAGGTTTTTTGGATGTCATGACATATTTTCCAt                                                          >  4:72329/1‑149 (MQ=255)
                                                                                               ccTATTTCAGCGCCTGCCTTGTTGGCGCGGGCTTGCCCGCTAATGCGTCCGATCCCACAACATTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                <  2:27098/149‑1 (MQ=38)
                                                                                                                               ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               <  4:193116/39‑1 (MQ=255)
                                                                                                                               ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                               >  3:193116/1‑39 (MQ=255)
                                                                                                                                  cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGTGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc             >  4:215236/1‑149 (MQ=255)
                                                                                                                                  cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc             >  6:133673/1‑149 (MQ=255)
                                                                                                                                     gcgAATGCGTCCTATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa          <  2:170622/149‑1 (MQ=255)
                                                                                                                                              tCCGATCAGACAACAATATTGTCTCGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTCATTAGCCAATTTTCGTGGAGTAAATGGCACCATACTCCATCCATGACAGAAAGCTCAGCTAAAca  >  4:172979/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACG  >  NC_002947/444949‑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: