Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,512 G→C 7.6% L72L (CTG→CTC)  PP_4203 → ubiquinone oxidoreductase electron transfer flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,5120GC 7.6% 164.8 / 2.7 52L72L (CTG→CTC) PP_4203ubiquinone oxidoreductase electron transfer flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (29/19);  new base C (2/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTCCGCCGCCTGCCGCCTGAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTACATCATCT  >  NC_002947/4748374‑4748660
                                                                                                                                          |                                                                                                                                                    
gTCCGCCGCCTGCCGCCTGAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACt                                                                                                                                            >  3:318740/1‑149 (MQ=255)
    gccgcctgccgcctGAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGa                                                                                                                                         >  1:238497/1‑148 (MQ=255)
         ctgccgcctgAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCAGACTGGAAAGAGc                                                                                                                                   >  4:46421/1‑149 (MQ=255)
            ccgcctgAAGCAGACCGCCGCCGCAGCCGGTAGCGAAATCAGCGGCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCg                                                                                                                                <  4:256322/149‑1 (MQ=255)
            ccgcctgAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCg                                                                                                                                <  1:65669/149‑1 (MQ=255)
                 tgAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTCTTTCCAGTCGGGGAACAGCTCgtt                                                                                                                              <  7:104131/146‑3 (MQ=255)
                 tgAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCACGCGCCCTGAACGAGCTCTTTCCAGTCGGGGAACAGCTCgtt                                                                                                                              >  8:104131/1‑144 (MQ=255)
                  gAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCAc                                                                                                                          >  1:331566/1‑149 (MQ=255)
                     gCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCACCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGc                                                                                                                       >  4:166545/1‑149 (MQ=255)
                           ggCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAaca                                                                                                                 >  2:314000/1‑149 (MQ=255)
                            gccgccGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGCACGAGCTGTGCCCCGACGGGAAAGAGCTCGGTGCACCGCTGAagac                                                                                                                >  4:206175/1‑149 (MQ=255)
                               gccgAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCg                                                                                                              >  6:20807/1‑148 (MQ=255)
                                ccgAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCGGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGaa                                                                                                            <  6:50632/149‑1 (MQ=255)
                                         gTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGTTCGGTGCACCGCTGAACACCGAAGTGAAgcgc                                                                                                   <  4:318740/149‑1 (MQ=255)
                                                    aGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCt                                                                                         <  1:314000/148‑1 (MQ=255)
                                                     gCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTAc                                                                                       >  7:156293/1‑149 (MQ=255)
                                                     gCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTAc                                                                                       <  3:206175/149‑1 (MQ=255)
                                                                  ggAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATg                                                                          >  7:117010/1‑149 (MQ=255)
                                                                   gagaAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATg                                                                          <  2:238497/148‑1 (MQ=255)
                                                                     gaAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCTCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCg                                                                       >  8:240256/1‑149 (MQ=255)
                                                                     gaAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCg                                                                       >  3:16477/1‑149 (MQ=255)
                                                                          gCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAgc                                                                  <  3:46421/149‑1 (MQ=255)
                                                                                   tCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAg                                                         >  7:358075/1‑149 (MQ=255)
                                                                                     ggCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCg                                                                       >  8:214439/1‑133 (MQ=255)
                                                                                     ggCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCg                                                                       <  7:214439/133‑1 (MQ=255)
                                                                                          ccACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTg                                                  >  3:280987/1‑149 (MQ=255)
                                                                                            acaTCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGAc                                                >  2:1755/1‑149 (MQ=255)
                                                                                                 ctctctGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGtt                                           >  5:10015/1‑149 (MQ=255)
                                                                                                 ctctctGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGt                                            >  8:149403/1‑148 (MQ=255)
                                                                                                  tctctGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCAGAGTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGCAGCGCGACGACATGTACGTACTCAAGGATGCCGGCAGGGCGCCCAAGCTGCCTGACCTGTTc                                          >  8:190822/1‑149 (MQ=255)
                                                                                                      tGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGt                                        <  1:1755/147‑1 (MQ=255)
                                                                                                        gcgcAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTg                                       <  8:117010/146‑1 (MQ=255)
                                                                                                        gcgcAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGcc                                     >  7:273244/1‑148 (MQ=255)
                                                                                                        gcgcAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGcc                                     >  5:170611/1‑148 (MQ=255)
                                                                                                        gcgcAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGcc                                     >  7:41371/1‑148 (MQ=255)
                                                                                                            aGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGa                                <  6:170611/149‑1 (MQ=255)
                                                                                                                ttCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGAc                                                <  3:299998/129‑1 (MQ=255)
                                                                                                                ttCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGAc                                                >  4:299998/1‑129 (MQ=255)
                                                                                                                  cGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATg                           <  3:166545/148‑1 (MQ=255)
                                                                                                                  cGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGc                          >  6:17828/1‑149 (MQ=255)
                                                                                                                   gAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGca                         >  4:357965/1‑149 (MQ=255)
                                                                                                                        ccGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGccc                                                                       >  4:298760/1‑97 (MQ=255)
                                                                                                                         cgcgCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCa                   >  6:221281/1‑149 (MQ=255)
                                                                                                                          gcgcCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCAt                  <  5:160234/149‑1 (MQ=255)
                                                                                                                           cgcCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATg                 <  4:16477/149‑1 (MQ=255)
                                                                                                                              ccTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGc               >  2:284265/1‑148 (MQ=255)
                                                                                                                                  accGAGCTCTTTCCAGTCGGGGAACAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTa          >  4:90446/3‑149 (MQ=255)
                                                                                                                                   ccGAGCTCTTTCCAGTCGGGGAACAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTAc         <  3:90446/148‑1 (MQ=255)
                                                                                                                                    cGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTAca        <  3:261880/149‑1 (MQ=255)
                                                                                                                                     gAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTAcat       <  7:149403/149‑1 (MQ=255)
                                                                                                                                     gAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAAAACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTAcat       <  8:358075/149‑1 (MQ=255)
                                                                                                                                       gCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTAcatca     >  4:383424/1‑149 (MQ=255)
                                                                                                                                          gTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTACATCAtct  >  1:3281/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GTCCGCCGCCTGCCGCCTGAAGCAGAAGGCCGCCGAAGCCGGTAGCGAAATCAGCGTCTGCGTGGTGGAGAAAGGCTCTGAAGTCGGCGCCCACATCCTCTCTGGCGCAGTGTTCGAGCCCCGCGCCCTGAACGAGCTGTTCCCCGACTGGAAAGAGCTCGGTGCACCGCTGAACACCGAAGTGAAGCGCGACGACATCTACGTACTCAAGGATGCCGGCAGCGCGACCAAGGTGCCTGACCTGTTCGTGCCCAAGACCATGCACAACCATGGCAACTACATCATCT  >  NC_002947/4748374‑4748660

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