Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,105 A→T 12.0% intergenic (+78/‑210) 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,1050AT12.0% 117.6 / 4.7 46intergenic (+78/‑210)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base T (2/4);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAA  >  NC_002947/444957‑445251
                                                                                                                                                    |                                                                                                                                                  
tCGATCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGAcaacaa                                                                                                                                                    <  8:98612/149‑1 (MQ=255)
    gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                >  4:29504/1‑149 (MQ=255)
        ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                              <  2:144513/147‑1 (MQ=255)
        ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                            <  3:29504/149‑1 (MQ=255)
        ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAATTCGTCCGATCAGACAAAAATTTTGtgt                                                                                                                                            >  7:27770/1‑147 (MQ=255)
                ggTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                    <  5:194308/149‑1 (MQ=255)
                          gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                          <  6:83789/149‑1 (MQ=18)
                               cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGTATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                         <  7:25324/129‑1 (MQ=255)
                               cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGTATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                         >  8:25324/1‑129 (MQ=255)
                                 acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGATTCCCGCGAATGCGGCCGAAGAGCCAAAAAATTTTTATGGTGGGGGCTCGTAGCGGGGCAGcc                                                                                                                   >  5:226943/1‑149 (MQ=255)
                                         cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                          <  1:283880/134‑1 (MQ=17)
                                         cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                          >  2:283880/1‑134 (MQ=17)
                                                 gATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGCATGCGTCCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCGCCTAAAAggggt                                                                                                   >  3:322540/1‑147 (MQ=255)
                                                   tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATAAGACAACAATTTTTTCTGGTCGGGCCTCTTAGCGGGCAATCCCGCTCCGACAACGGTTtg                                                                                                 >  6:101027/1‑149 (MQ=255)
                                                         acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGTTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                           >  5:88123/1‑149 (MQ=255)
                                                         acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                            <  2:342555/100‑1 (MQ=255)
                                                                    tGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                >  2:78828/1‑85 (MQ=255)
                                                                       aTACTTGCGTCATCCCTTATTTCTGCGCCTGCCTTGTTGGTGGGGGCTTGCCGGCTATTGCGTCCTATCACAAAAAATTTTTGTCTGGTCGGGCCTCTTCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                             <  2:341892/149‑1 (MQ=255)
                                                                           ttGCGTCATCCACTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATGCCCGCGAATGAGTCCCATAACACAACAAAATTTTCTGGTGGGGCCTCGTACCGGGCAAGCCAGGGCCTACAAGCGTTTGAGGCATTGCATGCCAGCTTTggc                                                                          >  2:263719/1‑146 (MQ=255)
                                                                                  aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGTTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                                  <  6:88123/149‑1 (MQ=255)
                                                                                    ccccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATTCGTCCCATAAGACAACAATTTTTTATGGTCGGGCCACTTAGCGGGCCAGCCCCCTCCGAAAAGGGTTTTTGGCATGGCAAGTCTGATTTTGCATTTTTAGTg                                                                >  2:32186/1‑149 (MQ=255)
                                                                                            tCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcg                                                                                                                             <  3:278652/80‑1 (MQ=255)
                                                                                            tCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcg                                                                                                                             >  4:278652/1‑80 (MQ=255)
                                                                                                 ccctgcctTGTTGCTTCGGGCTTGCCCGCGAATGCGTCCGATCATCAAACAATTTTGGCTGGTCGGGACTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                    <  1:263719/147‑1 (MQ=255)
                                                                                                               cgCGGGCTGCCCCGCTAAGGCGTCCGATCCCACAAAAATTTTGTCTGGCGGGGCCTCTCCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                     <  6:226943/148‑1 (MQ=255)
                                                                                                                 cGGGTTTGCCCGCTAATGCGGCCGATCACAAAACAATTTTGTCGGGGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                   <  8:27770/149‑1 (MQ=255)
                                                                                                                        gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                             >  6:347114/1‑132 (MQ=255)
                                                                                                                        gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                             <  5:347114/132‑1 (MQ=255)
                                                                                                                         cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCTTAAATGTCACCATAGTCCATCCTTGACAGAAAGc                           >  2:174557/1‑149 (MQ=255)
                                                                                                                         cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                            <  1:47634/148‑1 (MQ=255)
                                                                                                                         cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                           >  5:39090/1‑149 (MQ=18)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                      <  5:327631/137‑1 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                      >  6:327631/1‑137 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                   >  8:350400/1‑140 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                   <  7:350400/140‑1 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  4:21296/1‑149 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  5:275895/1‑149 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  5:134588/1‑149 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  1:229380/1‑149 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  4:178987/1‑149 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTACTACAAGGGTTTGTGGCATGGCATGACTGATTTGCAATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAgaa                          >  4:230645/1‑147 (MQ=255)
                                                                                                                            gctaaTGCGTCCGATCAGACAACAATTTTGTCTGGGCGGGCCTATTCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                        <  1:221988/146‑1 (MQ=255)
                                                                                                                            gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCa                                          <  5:423140/131‑1 (MQ=255)
                                                                                                                            gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCa                                          >  6:423140/1‑131 (MQ=255)
                                                                                                                            gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                          >  2:106113/1‑147 (MQ=255)
                                                                                                                                           cacacaaaaTTTTTGTTTGGCCGGGCCTCTCCGGGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt         <  2:365036/146‑1 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc        >  1:403752/1‑148 (MQ=255)
                                                                                                                                               caacaaTTTTTTCTGTTCGGGCCTCTTCGCGTGCAAGGCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa     <  3:230645/149‑1 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa     >  2:33581/1‑148 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa     >  2:27288/1‑148 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGAATTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa     >  1:101088/1‑148 (MQ=255)
                                                                                                                                                  caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                     >  6:379266/1‑50 (MQ=255)
                                                                                                                                                  caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                     <  5:379266/50‑1 (MQ=255)
                                                                                                                                                  caaTTTTGTCTGGTCGGGCCGCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa  <  5:101027/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAA  >  NC_002947/444957‑445251

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: