Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,481 T→C 6.9% M17T (ATG→ACG)  ccoO‑I → cbb3‑type cytochrome c oxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,837,4810TC 6.9% 164.7 / 2.6 58M17T (ATG→ACG) ccoO‑Icbb3‑type cytochrome c oxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base C (2/2);  total (32/26)
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.12e-01

CATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGAAGGTTGCGTGGG  >  NC_002947/4837334‑4837628
                                                                                                                                                   |                                                                                                                                                   
cATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATg                                                                                                                                                    <  4:54674/149‑1 (MQ=255)
        tGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGACggcaaata                                                                                                                                             >  5:265786/1‑142 (MQ=255)
             gCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGt                                                                                                                                         >  7:336248/1‑147 (MQ=255)
                  gAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGACGGCAAATACCATCAGCagg                                                                                                                                  <  6:265786/149‑3 (MQ=255)
                    acacCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATcgg                                                                                                                                >  6:375155/1‑149 (MQ=255)
                        cTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCggcggc                                                                                                                            >  8:274733/1‑149 (MQ=255)
                        cTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCggcggc                                                                                                                            <  8:336248/149‑1 (MQ=255)
                               aCCGTGCGTTCGCCAGCGCTTGACGTCGGCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTg                                                                                                                                             >  1:232536/1‑125 (MQ=255)
                               aCCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTg                                                                                                                                             <  2:232536/125‑1 (MQ=255)
                                      gTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGt                                                                                                              <  5:296628/149‑1 (MQ=255)
                                      gTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGt                                                                                                              >  4:68549/1‑149 (MQ=255)
                                            cAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCt                                                                                                        >  7:166298/1‑149 (MQ=255)
                                                  ttGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGttctt                                                                                                  >  6:66197/1‑149 (MQ=255)
                                                    gACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGttcttc                                                                                                 <  5:66197/148‑1 (MQ=255)
                                                    gACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGttcttc                                                                                                 >  4:188260/1‑148 (MQ=255)
                                                       gTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAg                                                                                              >  5:280419/1‑148 (MQ=255)
                                                           ccccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCa                                                                                                                                       >  3:123832/1‑103 (MQ=255)
                                                           ccccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCa                                                                                                                                       <  4:123832/103‑1 (MQ=255)
                                                           ccccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTACCGTCAGCATCGGCGGCCTGACCCAGATc                                                                                                                <  3:7253/126‑1 (MQ=255)
                                                           ccccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTACCGTCAGCATCGGCGGCCTGACCCAGATc                                                                                                                >  4:7253/1‑126 (MQ=255)
                                                            cccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCAACGGCGGCCTGACCCAGATCGTCCCGCTCTTCTTCCAGGa                                                                                            >  4:376813/1‑145 (MQ=255)
                                                            cccTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGACAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGa                                                                                            <  3:376813/145‑1 (MQ=255)
                                                              cTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGACGGCAAATACCATCAGCa                                                                                                                                    >  8:373333/1‑103 (MQ=255)
                                                              cTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGACGGCAAATACCATCAGCa                                                                                                                                    <  7:373333/103‑1 (MQ=255)
                                                                        ccAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTg                                                                                                       >  4:216689/1‑122 (MQ=255)
                                                                        ccAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTg                                                                                                       <  3:216689/122‑1 (MQ=255)
                                                                        ccAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCg                                                                            <  6:280419/149‑1 (MQ=255)
                                                                           gTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTc                                                                         <  6:282484/149‑1 (MQ=255)
                                                                             tGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCg                                                                        >  2:227192/1‑148 (MQ=255)
                                                                                 ttGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCa                                                                                               <  2:266141/121‑1 (MQ=255)
                                                                                 ttGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCa                                                                                               >  1:266141/1‑121 (MQ=255)
                                                                                 ttGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAgg                                                                    >  5:177888/1‑148 (MQ=255)
                                                                                        gaAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGc                                                            <  7:274733/149‑1 (MQ=255)
                                                                                                    gaAACATGAAGTCATCGAGAAAAACCTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCgcgc                                                >  3:2795/1‑149 (MQ=255)
                                                                                                     aaaCATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCa                                                                                       <  5:157289/109‑1 (MQ=255)
                                                                                                     aaaCATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCa                                                                                       >  6:157289/1‑109 (MQ=255)
                                                                                                       aCATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTg                                              >  4:210342/1‑148 (MQ=255)
                                                                                                        cATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGc                                             >  3:129324/1‑148 (MQ=255)
                                                                                                              gTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTgg                                       >  6:320885/1‑148 (MQ=255)
                                                                                                                 aTCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTgg                                       >  4:279657/1‑145 (MQ=255)
                                                                                                                 aTCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTgg                                       <  3:279657/145‑1 (MQ=255)
                                                                                                                    gagaAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGcc                                 <  5:320885/148‑1 (MQ=255)
                                                                                                                    gagaAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCcg                                >  6:184376/1‑149 (MQ=255)
                                                                                                                    gagaAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCcg                                >  8:244117/1‑149 (MQ=255)
                                                                                                                    gagaAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTAAGCATCCGCGCCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGACGAAGCCATGAAGCCCTACAACTCGCTGCAACTGGAACGCcg                                >  7:65235/1‑149 (MQ=25)
                                                                                                                        aaaaCGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGAt                            <  6:177888/149‑1 (MQ=255)
                                                                                                                        aaaaCGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGAt                            <  7:244117/149‑1 (MQ=255)
                                                                                                                         aaaCGTCGGCCTGATGGCCCTGCTGATGGTATTAGCCGTCAGCATCCGCGGCCTGACACAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACAGCGCTGCAACTGGAAGGCCGAGATg                           >  1:263373/1‑149 (MQ=25)
                                                                                                                            cGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCTGCCTGACCCAGATCTTCCCGCTGTTCTTCCAGGACGTCACCAACAATCCGGTCGAAGGCATGAAGCCCTACACGGCGCTGCAACTGGTAGGCCGGGATGTCt                        <  2:263373/149‑1 (MQ=25)
                                                                                                                            cGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTc                         >  6:368191/1‑148 (MQ=255)
                                                                                                                                 gCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGCCCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATc                   >  7:82817/1‑149 (MQ=255)
                                                                                                                                       tGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGaa             <  5:184376/149‑1 (MQ=255)
                                                                                                                                       tGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGaa             <  3:210342/149‑1 (MQ=255)
                                                                                                                                       tGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCGGTTCTTCCATTACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGaa             <  8:166298/149‑1 (MQ=255)
                                                                                                                                           ccTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACAt                    <  4:117829/138‑1 (MQ=255)
                                                                                                                                           ccTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACAt                    >  3:117829/1‑138 (MQ=255)
                                                                                                                                            ctgctgATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGAAGGTTg        >  5:397193/1‑149 (MQ=255)
                                                                                                                                                  aTGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGAAGGTTGCGTggg  <  3:188260/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGAAACATGAAGTCATCGAGAAAAACGTCGGCCTGATGGCCCTGCTGATGGTATTTGCCGTCAGCATCGGCGGCCTGACCCAGATCGTCCCGCTGTTCTTCCAGGACGTCACCAACAAACCGGTCGAAGGCATGAAGCCCTACACCGCGCTGCAACTGGAAGGCCGCGATGTCTACATCCGCGAAGGTTGCGTGGG  >  NC_002947/4837334‑4837628

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