Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 109,962 G→T 8.3% intergenic (+25/‑19) PP_0103 → / → ctaD cytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947109,9620GT 8.3% 145.5 / 2.7 49intergenic (+25/‑19)PP_0103/ctaDcytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (17/27);  new base T (2/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

AGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109813‑110107
                                                                                                                                                     |                                                                                                                                                 
agcagcTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATg                                                                                                                                                   <  1:109322/150‑1 (MQ=255)
       gTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGCAGCCGTCACAAGCTAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCa                                                                                                                                            <  2:710777/150‑1 (MQ=255)
               tCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAg                                                                                                                                    <  2:323507/150‑1 (MQ=255)
                 gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGa                                                                                                                                  <  2:540705/150‑1 (MQ=255)
                 gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGa                                                                                                                                  <  2:821713/150‑1 (MQ=255)
                 gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGACAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCCGAAGCCAAGTGAACCGGGTTACTCCCAGAAGACCAGAGATTGCAGGAGACAGGa                                                                                                                                  >  2:257178/1‑150 (MQ=255)
                         tgctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTg                                                                                                                          >  1:559381/1‑150 (MQ=255)
                           ctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCa                                                                                                                        <  1:878179/150‑1 (MQ=255)
                             gTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGt                                                                                                                      <  2:1070975/150‑1 (MQ=255)
                             gTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGt                                                                                                                      <  1:685207/150‑1 (MQ=255)
                              tAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTg                                                                                                                     >  1:541375/1‑150 (MQ=255)
                                 tCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATc                                                                                                                  <  2:373287/150‑1 (MQ=255)
                                   aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAg                                                                                                                                           <  2:694631/123‑1 (MQ=255)
                                   aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAg                                                                                                                                           >  1:694631/1‑123 (MQ=255)
                                    ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  1:538028/150‑1 (MQ=255)
                                    ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  2:559381/150‑1 (MQ=255)
                                    ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  1:305568/150‑1 (MQ=255)
                                     cTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgacg                                                                                                              >  2:36534/1‑150 (MQ=255)
                                          aGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCAc                                                                                                         <  2:603601/150‑1 (MQ=255)
                                                  gCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGACCACgg                                                                                                 <  1:535161/150‑1 (MQ=255)
                                                              aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACa                                                                                                                                     >  2:234976/1‑102 (MQ=255)
                                                              aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACa                                                                                                                                     <  1:234976/102‑1 (MQ=255)
                                                              aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGt                                                                                                                           >  1:788142/1‑112 (MQ=255)
                                                              aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGt                                                                                                                           <  2:788142/112‑1 (MQ=255)
                                                               acaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCAAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATgcg                                                                                    >  2:942564/1‑150 (MQ=255)
                                                                  aaGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCAc                                                                                 <  1:692542/150‑1 (MQ=255)
                                                                            ggTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCa                                                                       >  1:727968/1‑150 (MQ=255)
                                                                                             tGCTGGCGCAGAAGCAGGCAGAAGACAAGTGAACCGGGTTACTCGCGACAATAAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCATGATGCGCTgg                                                      <  2:541375/150‑1 (MQ=255)
                                                                                              gCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTggg                                                     <  2:727968/150‑1 (MQ=255)
                                                                                                      tGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGacc                                             >  1:537199/1‑150 (MQ=255)
                                                                                                          gcaggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGAccacca                                         >  2:153871/1‑150 (MQ=255)
                                                                                                          gcaggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGAccacca                                         >  2:155781/1‑150 (MQ=255)
                                                                                                            aggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaac                                       <  2:895590/150‑1 (MQ=255)
                                                                                                              gcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaa                                        >  1:1108096/1‑147 (MQ=255)
                                                                                                              gcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaa                                        <  2:1108096/147‑1 (MQ=255)
                                                                                                                  aaGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATaaa                                 <  2:537199/150‑1 (MQ=255)
                                                                                                                            aaCCGGGTTACTCGCGGCAATCAATGGATTGAAGGAGACAGGACATAAGTACAGTGAGCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCa                       >  1:415148/1‑150 (MQ=255)
                                                                                                                                gggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      >  1:988428/1‑99 (MQ=255)
                                                                                                                                gggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      <  2:988428/99‑1 (MQ=255)
                                                                                                                                 ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGAt                                                             <  2:1129226/107‑1 (MQ=255)
                                                                                                                                 ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGAt                                                             >  1:1129226/1‑107 (MQ=255)
                                                                                                                                       tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtg            >  2:684590/1‑150 (MQ=255)
                                                                                                                                       tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtg            <  1:36534/150‑1 (MQ=255)
                                                                                                                                        tcctgCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               <  1:850032/110‑1 (MQ=255)
                                                                                                                                        tcctgCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               >  2:850032/5‑114 (MQ=255)
                                                                                                                                         gcgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtgtg          >  2:917891/1‑150 (MQ=255)
                                                                                                                                          ctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGac                                              >  2:213109/3‑113 (MQ=255)
                                                                                                                                          ctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGac                                              <  1:213109/111‑1 (MQ=255)
                                                                                                                                               atcaatGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAg    <  1:942564/150‑1 (MQ=255)
                                                                                                                                                 caatGGATTGCAGGAGACAGGTCATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCCGCCACGCGCCTGATGCCCTAGCTGCTGACCACCAACCATAAAGACATCTGCACGATGTACCTGTGGTTCAGtt  <  2:49654/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109813‑110107

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