Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,836,864 A→G 8.9% M291V (ATG→GTG)  ccoN‑I → cbb3‑type cytochrome c oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,836,8640AG 8.9% 124.2 / 3.1 45M291V (ATG→GTG) ccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

TTCTGGGCGCTGATCACCCTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGTCCCACTACACC  >  NC_002947/4836723‑4837010
                                                                                                                                             |                                                                                                                                                  
ttCTGGGCGCTGATCACCCTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGAc                                                                                                                                             >  1:235422/1‑149 (MQ=255)
ttCTGGGCGCTGATCACCCTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGAc                                                                                                                                             >  3:361380/1‑149 (MQ=255)
             tCACCCTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGcc                                                                                                                                <  7:197216/149‑1 (MQ=255)
                cccTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCat                                                                                                                                                   >  8:19338/1‑127 (MQ=255)
                cccTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCat                                                                                                                                                   <  7:19338/127‑1 (MQ=255)
                cccTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTgg                                                                                                                             <  2:235422/149‑1 (MQ=255)
                 ccTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTgg                                                                                                                             >  4:20569/1‑148 (MQ=255)
                        aTCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACt                                                                                                                     >  1:160656/1‑149 (MQ=255)
                         tCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTg                                                                                                                    >  7:56895/1‑149 (MQ=255)
                             ggccggccCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTgc                                                                                                                   <  3:20569/146‑1 (MQ=255)
                                      gCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACACACTGCGCAGCGACCCGa                                                                                                       >  2:40920/1‑149 (MQ=255)
                                      gCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGa                                                                                                       >  6:377095/1‑149 (MQ=255)
                                      gCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGa                                                                                                       <  7:269828/149‑1 (MQ=255)
                                        accacTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATc                                                                                                     <  5:377095/149‑1 (MQ=255)
                                        accacTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATc                                                                                                     >  6:138795/1‑149 (MQ=255)
                                        accacTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATc                                                                                                     >  1:324657/1‑149 (MQ=255)
                                            cTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTc                                                                                                 >  3:132610/1‑149 (MQ=255)
                                            cTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTc                                                                                                 <  6:54460/149‑1 (MQ=255)
                                              tGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGc                                                                                               >  5:18595/1‑149 (MQ=255)
                                               tcACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCt                                                                                              <  7:209988/148‑1 (MQ=255)
                                               gCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGc                                                                                               <  2:324657/148‑1 (MQ=255)
                                                 aCTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCt                                                                                              >  8:209988/1‑147 (MQ=255)
                                                   tACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCGTGATGACCCTGTCGGGCGCCTGGCACAACCTGCGCAGCGACCCGATCCTTCGCTTCCt                                                                                          >  3:91430/1‑149 (MQ=255)
                                                   tACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCGTGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCt                                                                                          >  3:86322/1‑149 (MQ=255)
                                                      acCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCtggt                                                                                       >  3:213974/1‑149 (MQ=255)
                                                        cGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCtggtgg                                                                                     >  3:282253/1‑149 (MQ=255)
                                                          cccTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGgtggtg                                                                                   <  2:20965/149‑1 (MQ=255)
                                                           ccTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGgtgt                                                                                  >  4:282718/1‑149 (MQ=255)
                                                                      gggCGCAGTCGCTGGGCATGGTGATNTCGCTGATCNTGCTNGCGCCAAGCTGGGGCGGCATGATCAACGGCGTGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGtt                                                                        <  4:91430/148‑1 (MQ=255)
                                                                      gggCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCGTGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGtt                                                                        <  4:86322/148‑1 (MQ=255)
                                                                                gCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGt                                                             <  6:34489/149‑1 (MQ=255)
                                                                                  tGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTcc                                                           <  4:132610/149‑1 (MQ=255)
                                                                                    ggCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCAc                                                         <  4:361380/149‑1 (MQ=255)
                                                                                                   cTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCgat                                          <  4:213974/149‑1 (MQ=39)
                                                                                                       tCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGatgatg                                      >  4:251156/1‑149 (MQ=37)
                                                                                                        ccTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATgg                                     >  7:34615/1‑149 (MQ=37)
                                                                                                          accgggCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGc                                    <  1:38297/144‑1 (MQ=37)
                                                                                                          accgggCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGc                                    >  2:38297/5‑148 (MQ=37)
                                                                                                                 gcCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGa                            <  3:251156/149‑1 (MQ=37)
                                                                                                                     aGCTGGGGCGGCATGATCAACGGCATGATGACGCTGTCGGTCGCCTGGCACAAACTGCGCATCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGt                        <  1:40920/149‑1 (MQ=255)
                                                                                                                      gCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTc                       <  8:34615/149‑1 (MQ=40)
                                                                                                                                cATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGt             <  5:138795/149‑1 (MQ=255)
                                                                                                                                     tCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGaa                                                  <  7:266146/107‑1 (MQ=40)
                                                                                                                                     tCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGaa                                                  >  8:266146/1‑107 (MQ=40)
                                                                                                                                     tCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGTCCCa         >  6:190043/1‑148 (MQ=255)
                                                                                                                                     tCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGTCCCAc        <  2:160656/149‑1 (MQ=255)
                                                                                                                                            cATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGTCCCACTACAcc  >  8:3119/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TTCTGGGCGCTGATCACCCTGTACATCTGGGCCGGCCCGCACCACTTGCACTACACCGCCCTGCCCGATTGGGCGCAGTCGCTGGGCATGGTGATGTCGCTGATCCTGCTGGCGCCAAGCTGGGGCGGCATGATCAACGGCATGATGACCCTGTCGGGCGCCTGGCACAAACTGCGCAGCGACCCGATCCTTCGCTTCCTGGTGGTGTCACTGGCGTTCTACGGCATGTCCACCTTCGAAGGCCCGATGATGGCCATCAAGACGGTCAACGCCCTGTCCCACTACACC  >  NC_002947/4836723‑4837010

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