Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,095 T→C 17.2% 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,0950TC17.2% 99.5 / 12.4 49intergenic (+68/‑220)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base C (2/6);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTT  >  NC_002947/444949‑445243
                                                                                                                                                   |                                                                                                                                                    
gACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCa                                                                                                                                                    <  1:341528/149‑1 (MQ=255)
          gAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                          >  4:189648/1‑149 (MQ=255)
            gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  3:189648/149‑1 (MQ=255)
              tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTg                                                                                                                                       >  8:251460/1‑148 (MQ=255)
              tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAAGAGGCCCGACCAGACAAAATtgttgt                                                                                                                                      >  7:174226/1‑149 (MQ=37)
                ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                     <  7:261096/148‑1 (MQ=255)
                      aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                              <  7:26976/149‑1 (MQ=255)
                      aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCgg                                                                                                                              <  3:172705/149‑1 (MQ=38)
                       aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCgga                                                                                                                             <  5:69949/149‑2 (MQ=38)
                       aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCgga                                                                                                                             <  4:202749/149‑2 (MQ=38)
                        ggTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                            <  3:79293/149‑1 (MQ=255)
                               aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                      <  3:366860/148‑1 (MQ=255)
                                  gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCGGAATG‑GTGCGATCAGACAACAAATTTTTGGGGTGGGGGCGCTTcggcg                                                                                                                  >  4:33868/1‑146 (MQ=255)
                                   cATCCGACGACATGCTTGGCCGGATCATCGACTCCTACGCCTGATTACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                  <  8:174226/148‑1 (MQ=17)
                                       cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAg                                                                                                             <  4:177834/149‑1 (MQ=17)
                                         acgaTATGCTTGGCAGGATCATCGTCACCTACGCCTGAATACTTGCGTCATGCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCTATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGTCAAGcc                                                                                                           <  2:210303/149‑1 (MQ=39)
                                                    ggCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                    >  5:22090/1‑113 (MQ=255)
                                                    ggCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                    <  6:22090/113‑1 (MQ=255)
                                                           tcatcGACACCTCCGCCTGAATACTTGCGCCATCCCCTATTTCCGCCCCTGCCTTGTCGGCGCGGGC‑TTGCCCGCTATTGCGCCCGACCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGCCAAGCCCGCTCCTACAAGGGTTtg                                                                                         <  2:123212/149‑1 (MQ=255)
                                                                 acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATAAGACAACAATTTTTTCTGGTGGGGCCTCTTAGCGGGCCATCCCGCGCCAACAAGGGTTTGttgcat                                                                                   >  5:126923/1‑149 (MQ=255)
                                                                     cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                      >  8:200494/1‑110 (MQ=255)
                                                                     cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                      <  7:200494/110‑1 (MQ=255)
                                                                        cGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAAAACCTCCGATAACACAACAAATTTTTATTGTAGGGCCACTTAGCGGGGAAGCCCGCAGCAAAAAAGGTTTGTGgcatggcatg                                                                             >  1:118506/1‑148 (MQ=255)
                                                                         gCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                           <  7:251460/149‑1 (MQ=255)
                                                                             gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGAGAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                              <  6:225343/110‑1 (MQ=255)
                                                                             gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                              >  5:225343/1‑110 (MQ=255)
                                                                                           tCCCCTATTCCAGCGCCGCCCTGGGAGGAGCGGGC‑TGCCCCGCGAAGGCGTCCGCGCCTCCAACAATTTTGGCTGGGCGGGCCGCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                         <  6:285565/149‑1 (MQ=255)
                                                                                                      cccggctgcctTGTAGGAGCGGTG‑AGGCCCGGTAATGGTGCCGATCCGCCAAAAATTTTGTCGGGGCGGGCGTCTTCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATg                                               <  4:264385/143‑1 (MQ=255)
                                                                                                        cgcATGCCTTGTAGGAGCGGGC‑TTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                            <  2:118506/149‑1 (MQ=255)
                                                                                                           ctgcctTGTAGGAGCGGGC‑TTGCCCGCGATTGGGTCCTATCAGACAACAATTTTGTCTGGTGGGGCCTCTTCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCa                                         <  6:126923/149‑1 (MQ=38)
                                                                                                             gcattGTAGGAGCGGGC‑TTGCCCGCGAATGCGCCCGATCAGACAACAATTTTGTCGGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                       <  1:88468/146‑1 (MQ=255)
                                                                                                                         cGGGC‑TTCCCCGCTATTGCGCCCGAACACAAAAAAATTTTGTCTGTGCGGCCCTCTTCCCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                           <  8:336129/149‑1 (MQ=255)
                                                                                                                           ggC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                            >  5:19242/1‑130 (MQ=255)
                                                                                                                           ggC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                           >  1:176777/1‑147 (MQ=255)
                                                                                                                           ggC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGCCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                           >  1:106321/1‑147 (MQ=255)
                                                                                                                            gC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                           >  4:348018/1‑146 (MQ=255)
                                                                                                                                  cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                                <  1:363410/136‑1 (MQ=255)
                                                                                                                                  cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                                >  2:363410/1‑136 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTTCAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATgg                                              <  5:66801/121‑1 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTTCAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATgg                                              >  6:66801/1‑121 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                           >  5:363920/1‑92 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                           <  6:363920/92‑1 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                       >  3:360991/1‑128 (MQ=255)
                                                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                       <  4:360991/128‑1 (MQ=255)
                                                                                                                                     gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                  >  6:282056/1‑147 (MQ=255)
                                                                                                                                     gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                <  3:348018/149‑1 (MQ=255)
                                                                                                                                       gAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGc              >  2:338758/1‑149 (MQ=255)
                                                                                                                                              tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg        >  2:92127/1‑148 (MQ=255)
                                                                                                                                               ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt       >  1:352383/1‑148 (MQ=38)
                                                                                                                                                   tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt  >  3:384929/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTT  >  NC_002947/444949‑445243

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