Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,351,736 A→G 7.8% L5S (TTG→TCG)  PP_3823 ← cytochrome c‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,351,7360AG 7.8% 136.9 / 2.4 51L5S (TTG→TCG) PP_3823cytochrome c‑type protein
Reads supporting (aligned to +/- strand):  ref base A (26/21);  new base G (2/2);  total (28/23)
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.20e-01

GCCATTGCCGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAA  >  NC_002947/4351594‑4351883
                                                                                                                                              |                                                                                                                                                   
gCCATTGCCGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGgtg                                                                                                                                               >  1:261165/1‑149 (MQ=255)
 ccATTGCCGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCTGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGgtg                                                                                                                                               <  5:90625/148‑1 (MQ=255)
        cGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCAc                                                                                                                                        >  1:390238/1‑148 (MQ=255)
        cGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCAc                                                                                                                                        >  1:411812/1‑148 (MQ=255)
         gTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCGACGGTGTCATCACg                                                                                                                                       >  7:81536/1‑148 (MQ=255)
         gTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCGACGGTGTCATCACg                                                                                                                                       >  7:87374/1‑148 (MQ=255)
              ggCGGTGGCAAGCCAGGCAGGCCTTGGGGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAAtctc                                                                                                                                 <  7:344064/149‑1 (MQ=255)
              ggCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAAtctc                                                                                                                                 <  5:235082/149‑1 (MQ=255)
                  gTGGCAACCCAGGCAGGCCATGGCGGGGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCcttt                                                                                                                             >  1:232946/1‑149 (MQ=255)
                   tGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCctttt                                                                                                                            >  2:22434/1‑149 (MQ=255)
                       aaCCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGgtgt                                                                                                                                              >  3:41780/1‑127 (MQ=255)
                       aaCCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGgtgt                                                                                                                                              <  4:41780/127‑1 (MQ=255)
                       aaCCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCggtttgttt                                                                                                                        >  6:66769/1‑149 (MQ=255)
                          ccAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGt                                                                                                                     >  8:134368/1‑149 (MQ=255)
                                   ccATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGttt                                                                                                            <  4:395902/149‑1 (MQ=255)
                                      tGGCGGCGGGGCTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAgg                                                                                                         <  8:399775/149‑1 (MQ=255)
                                           gcggGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTaa                                                                                                    <  2:261165/149‑1 (MQ=255)
                                              gggTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAAc                                                                                                 <  7:82791/149‑1 (MQ=255)
                                              gggTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAAc                                                                                                 <  7:259530/149‑1 (MQ=255)
                                                gTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACat                                                                                               <  4:69031/149‑1 (MQ=255)
                                                   ggCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACatat                                                                                             <  2:59464/148‑1 (MQ=255)
                                                    gCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCGACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACAtata                                                                                            <  8:87374/148‑1 (MQ=255)
                                                    gCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCGACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACAtata                                                                                            <  8:81536/148‑1 (MQ=255)
                                                      cccccccTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAg                                                                                          >  6:38974/1‑148 (MQ=255)
                                                      cccccccTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAg                                                                                          >  4:419478/1‑148 (MQ=255)
                                                                gTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAagga                                                                               <  5:66769/149‑1 (MQ=255)
                                                                    aGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGAc                                                                           <  5:279692/149‑1 (MQ=255)
                                                                         ttCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTTTATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGtt                                                                      >  5:284558/1‑149 (MQ=255)
                                                                         ttCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGtt                                                                      >  2:88784/1‑149 (MQ=255)
                                                                              aCCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAACGGAGGACCCGTTGCtt                                                                  <  1:392187/148‑1 (MQ=255)
                                                                                 gTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCGGTTGCTTCAGa                                                              <  2:232946/149‑1 (MQ=255)
                                                                                   ccgcggcACATGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACatat                                                                                             <  1:329816/114‑1 (MQ=255)
                                                                                   ccgcggcACATGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACatat                                                                                             >  2:329816/3‑116 (MQ=255)
                                                                                           tgtgCCGCGGGCCAGTACATGATCAGCAAGCTGCCGATCATCACTCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGATTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCg                                                    <  2:411812/149‑1 (MQ=255)
                                                                                            gtgCCGCGGGCCAGGCCAGGTTCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCg                                                    <  2:91374/148‑1 (MQ=255)
                                                                                             tgCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATg                                                                                                               >  4:49803/1‑88 (MQ=255)
                                                                                             tgCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATg                                                                                                               <  3:49803/88‑1 (MQ=255)
                                                                                               ccGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGTTTCAGATCAACTTCCg                                                    >  7:70680/1‑145 (MQ=255)
                                                                                               ccGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTc                                                                 >  4:168279/1‑132 (MQ=255)
                                                                                               ccGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCg                                                    >  3:258742/1‑145 (MQ=255)
                                                                                               ccGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCg                                                    >  2:243253/1‑145 (MQ=255)
                                                                                                cgcgGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGAt                                               >  8:290106/1‑149 (MQ=255)
                                                                                                          gccaggATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGtt                                     >  7:233587/1‑149 (MQ=255)
                                                                                                             aggATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGTTTCAGATCAACTTCCGGAGATGCATTGCGTTGcc                                  <  8:70680/149‑1 (MQ=255)
                                                                                                               gATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCg                                 >  4:404517/1‑148 (MQ=255)
                                                                                                                  cAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGt                             >  3:338307/1‑149 (MQ=255)
                                                                                                                    gCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTgc                           >  1:2666/1‑149 (MQ=255)
                                                                                                                           gCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATc                    >  3:362907/1‑149 (MQ=255)
                                                                                                                            ccGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCg                   <  4:111561/149‑1 (MQ=255)
                                                                                                                              gATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGc                  >  1:145275/1‑148 (MQ=255)
                                                                                                                               atcatcACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGcc                 >  5:124286/1‑148 (MQ=255)
                                                                                                                                              aaCGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCggg    <  3:404517/146‑1 (MQ=255)
                                                                                                                                              aaCGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCggg    <  1:22434/146‑1 (MQ=255)
                                                                                                                                              aaCGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGaa  <  6:124286/148‑1 (MQ=255)
                                                                                                                                              aaCGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGaa  <  5:38974/148‑1 (MQ=255)
                                                                                                                                              aaCGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGaa  <  3:419478/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCCATTGCCGTCCGGGCCGTGGCAACCCAGGCAGGCCATGGCGGCGGGGTTGGCCCCCCCCTGGGTGAAGACCTTCTGACCGTCGGCGGCATGTGCCGCGGGCCAGGCCAGGATCAGCAAGCTGCCGATCATCACGCGTCTCAACGGTGTCATCACGAATCTCCTTTTCTTTTGTTATATGTTTAGGCTTAAAACATATAAGCAGAGAAAGGAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAA  >  NC_002947/4351594‑4351883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑