Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,241 G→A 9.2% G16E (GGG→GAG)  cyoC → cytochrome bo terminal oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947953,2410GA 9.2% 118.5 / 2.0 43G16E (GGG→GAG) cyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base A (2/2);  total (20/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 = 1.00e+00

CGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCA  >  NC_002947/953092‑953385
                                                                                                                                                     |                                                                                                                                                
cGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTATGCCTTCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTTAAGTAATGCACGGCGCTGCTCATGGTCACGACCATgg                                                                                                                                                  <  2:347326/150‑1 (MQ=255)
 gaagaaGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATggg                                                                                                                                                 >  1:418036/1‑150 (MQ=255)
 gaagaaGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATggg                                                                                                                                                 <  2:194431/150‑1 (MQ=255)
                aTCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccaccac                                                                                                                                  >  1:555288/1‑150 (MQ=255)
                  cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGAccaccacca                                                                                                                                <  1:94233/150‑1 (MQ=255)
                     aGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGTACGTGTCGAATCGTTTGAGCGGGTGTAATCAATTTCCAGTCAAGTAATGCCCCGCCCTGCTCATGGTCACTACCATGGGCATGACCACCACCACCACGa                                                                                                                             <  2:313770/150‑1 (MQ=255)
                     aGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGa                                                                                                                             >  2:693860/1‑150 (MQ=255)
                     aGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGa                                                                                                                             <  2:658940/150‑1 (MQ=255)
                         gAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATgacgac                                                                                                                                        >  2:529086/1‑135 (MQ=255)
                         gAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATgacgac                                                                                                                                        <  1:529086/135‑1 (MQ=255)
                          aGCGTATGAAAGCGCTGGCAAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                        <  2:582713/150‑1 (MQ=255)
                                  aaaGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATg                                                                                                                >  2:673617/1‑150 (MQ=255)
                                           gCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTg                                                                                                       <  1:312471/150‑1 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCagc                                                                                                                                             >  1:128354/1‑108 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCagc                                                                                                                                             <  2:128354/110‑3 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCagc                                                                                                                                             <  2:128292/110‑3 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCagc                                                                                                                                             >  1:128292/1‑108 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCCGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTggg                                                                                                     >  2:657451/1‑150 (MQ=255)
                                                  cAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCt                                                                                                <  2:576696/150‑1 (MQ=255)
                                                  cAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCt                                                                                                <  2:555288/150‑1 (MQ=255)
                                                     gTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGc                                                                                             <  2:499203/150‑1 (MQ=255)
                                                         tctGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGc                                                                                                                                                >  1:507235/1‑95 (MQ=255)
                                                         tctGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGc                                                                                                                                                <  2:507235/95‑1 (MQ=255)
                                                               tgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtgat                                                                                   <  1:693860/150‑1 (MQ=255)
                                                                     cgcACGTGTCGCATCGTTTGCGCGGGTGTAATCAACGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGa                                                                             <  1:713479/150‑1 (MQ=255)
                                                                       cACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACt                                                                           <  1:611050/150‑1 (MQ=255)
                                                                              cGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCt                                                                    <  1:688142/150‑1 (MQ=255)
                                                                                aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                                  <  2:418036/150‑1 (MQ=255)
                                                                                 aTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                        <  1:305190/95‑1 (MQ=255)
                                                                                 aTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                        >  2:305190/1‑95 (MQ=255)
                                                                                  tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTc                                                                                                 <  1:693204/117‑1 (MQ=255)
                                                                                  tCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTc                                                                                                 >  2:693204/1‑117 (MQ=255)
                                                                                                    tCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                                  >  1:330523/1‑130 (MQ=255)
                                                                                                    tCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                                  <  2:330523/130‑1 (MQ=255)
                                                                                                             aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGccacccaccacgtg                                     >  2:399191/1‑144 (MQ=255)
                                                                                                               tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                    >  2:215904/1‑133 (MQ=255)
                                                                                                               tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                    <  1:215904/133‑1 (MQ=255)
                                                                                                                      aTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGAGTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGc                            >  2:535639/1‑150 (MQ=255)
                                                                                                                      aTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGc                            >  2:683992/1‑150 (MQ=255)
                                                                                                                              cgcTGCTCATGGTCACGTTGCCTGTCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTCCCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGc                    <  1:243049/150‑1 (MQ=255)
                                                                                                                                       tcgTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGccc           >  1:468866/3‑150 (MQ=255)
                                                                                                                                             cGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGAGCCTCTGGCTGTACCTGATGACCGACTGCCTCCTGTTTGCCGCGCTCCTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGCCGGCCCGtggcg     >  2:74941/1‑146 (MQ=255)
                                                                                                                                              gACCATGGGCATGACGACCACTACCACGACTCTGGCCAGAGGCCCTCACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGATGCTGTCCGGCCGTTTTTCCGTCGGCCCGTCgggc    >  2:122249/1‑149 (MQ=255)
                                                                                                                                                ccATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCa  >  1:135099/1‑150 (MQ=255)
                                                                                                                                                 cATGGGCATGACGACCACCACCACGACTCGGGACAGATGACCGTACTGGGCTTCTgg                                                                                              <  2:250805/57‑1 (MQ=255)
                                                                                                                                                 cATGGGCATGACGACCACCACCACGACTCGGGACAGATGACCGTACTGGGCTTCTgg                                                                                              >  1:250805/1‑57 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCA  >  NC_002947/953092‑953385

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