Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,224 T→C 8.9% A10A (GCT→GCC 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,2240TC 8.9% 121.3 / 2.3 45A10A (GCT→GCCcyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base C (2/2);  total (27/18)
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.33e-01

CATGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGC  >  NC_002947/953080‑953371
                                                                                                                                                |                                                                                                                                                   
cATGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCCAGTAATGCACGGCGCTGCTCATgg                                                                                                                                                <  2:752471/150‑1 (MQ=255)
 aTGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCa                                                                                                                                                   >  1:548424/1‑146 (MQ=255)
 aTGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCa                                                                                                                                                   <  2:548424/146‑1 (MQ=255)
 aTGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGt                                                                                                                                               >  1:417911/1‑150 (MQ=255)
     ttCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACg                                                                                                                                           >  2:1080531/1‑150 (MQ=255)
            cGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATgg                                                                                                                                    <  2:417911/150‑1 (MQ=255)
                  agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATga                                                                                                                              <  2:954559/150‑1 (MQ=255)
                       cgcgcATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGacc                                                                                                                         <  2:309968/150‑1 (MQ=255)
                            aTCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGAGGaccaccac                                                                                                                    >  2:400249/1‑150 (MQ=255)
                              cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATgag                                                                                                                                   <  1:313628/133‑3 (MQ=255)
                              cGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATgag                                                                                                                                   >  2:313628/1‑131 (MQ=255)
                               gAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCACGAGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACcaccaccac                                                                                                                 >  1:261863/1‑150 (MQ=255)
                                             gAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGAt                                                                                                   >  1:747610/1‑150 (MQ=255)
                                                 gcgcTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGAcc                                                                                               <  1:1080531/150‑1 (MQ=255)
                                                        ccAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGCGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCAc                                                                                                                                            >  1:324401/1‑98 (MQ=255)
                                                        ccAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCAc                                                                                                                                            <  2:324401/98‑1 (MQ=255)
                                                         cAAGGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTggg                                                                                       >  2:1128244/1‑150 (MQ=255)
                                                         cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATg                                                                                                                                                 <  2:92262/92‑1 (MQ=255)
                                                         cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATg                                                                                                                                                 >  1:92262/1‑92 (MQ=255)
                                                             gCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCCGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCCCCTCCACGACGCGGGCCAGATGACCGTACTGGGCTTc                                                                                   <  2:407572/150‑1 (MQ=255)
                                                                        gcctgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCt                                                                        >  1:614026/1‑150 (MQ=255)
                                                                          ctgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccac                                                                                                                       <  1:972409/101‑1 (MQ=255)
                                                                          ctgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccac                                                                                                                       >  2:972409/1‑101 (MQ=255)
                                                                          ctgccGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCtga                                                                      <  1:152642/150‑1 (MQ=255)
                                                                                 cgcACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGTGCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGa                                                               <  1:353043/150‑1 (MQ=255)
                                                                                  gcACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCt                                                                        >  2:708381/1‑140 (MQ=255)
                                                                                  gcACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCt                                                                        <  1:708381/140‑1 (MQ=255)
                                                                                    aCGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGTGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTg                                                            >  2:195001/1‑150 (MQ=255)
                                                                                           gAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTg                                                     >  2:671675/1‑150 (MQ=255)
                                                                                            aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTgg                                                                                        <  1:440710/114‑1 (MQ=255)
                                                                                            aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTgg                                                                                        >  2:440710/1‑114 (MQ=255)
                                                                                            aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                    >  1:1103793/1‑150 (MQ=255)
                                                                                            aaTCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGt                                                    >  1:1103819/1‑150 (MQ=255)
                                                                                                gTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGc                                                >  1:444304/1‑150 (MQ=255)
                                                                                                 tttGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCg                                               >  1:142185/1‑150 (MQ=255)
                                                                                                         ggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCa                                                                                                                                 <  2:471542/60‑1 (MQ=25)
                                                                                                         ggTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCCCATGGTCGTGACCATGAGCa                                                                                                                                 >  1:471542/1‑60 (MQ=21)
                                                                                                             tAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGcc                                   >  1:616504/1‑150 (MQ=255)
                                                                                                                tCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAACGGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAcc                                >  1:560835/1‑150 (MQ=255)
                                                                                                                      tCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTa                              <  2:66693/146‑1 (MQ=255)
                                                                                                                      tCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTa                              >  1:66693/1‑146 (MQ=255)
                                                                                                                           tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGa                                                                                                               >  2:501453/1‑60 (MQ=255)
                                                                                                                           tCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGa                                                                                                               <  1:501453/60‑1 (MQ=255)
                                                                                                                            cAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTg                    >  2:621138/1‑150 (MQ=255)
                                                                                                                                   tGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCa             >  2:664475/1‑150 (MQ=255)
                                                                                                                                              gctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGGTGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCcggc  >  1:194992/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATGGTTCCGGCCGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGC  >  NC_002947/953080‑953371

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