Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,834,069 A→G 7.0% H157R (CAC→CGC)  ccoN1_2 → Cbb3‑type cytochrome c oxidase subunit CcoN1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,834,0690AG 7.0% 174.5 / 4.4 57H157R (CAC→CGC) ccoN1_2Cbb3‑type cytochrome c oxidase subunit CcoN1
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

GCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCGATGA  >  ACB122/4833941‑4834203
                                                                                                                                |                                                                                                                                      
gCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTa                                                                                                                                  >  3:294196/1‑135 (MQ=255)
   gATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGt                                                                                                                               >  5:125472/1‑135 (MQ=255)
       tgCTACACCACCCCCAAGGAATACGCCGAGCTGGAGTGGCCGCGGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGGCAAACGCAAGACCAAGCACATCTATGTCGGc                                                                                                                           <  6:17602/134‑1 (MQ=255)
        gCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                                                                                                                           <  2:220784/134‑1 (MQ=255)
        gCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                                                                                                                           <  3:145620/134‑1 (MQ=255)
              ccaccaccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGt                                                                                                                    >  7:182238/1‑135 (MQ=255)
                 ccaccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCt                                                                                                                 <  5:22096/135‑1 (MQ=255)
                    ccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACg                                                                                                              >  1:68860/1‑135 (MQ=255)
                    ccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACg                                                                                                              >  1:88110/1‑135 (MQ=255)
                    ccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACg                                                                                                              >  4:145543/1‑135 (MQ=255)
                    ccaAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACg                                                                                                              >  3:82730/1‑135 (MQ=255)
                         gAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACGGGTTCTACGGTGCg                                                                                                         >  6:290894/1‑135 (MQ=255)
                          aaTACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGt                                                                                                        >  1:197107/1‑135 (MQ=255)
                               gCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATc                                                                                                   >  1:16499/1‑135 (MQ=255)
                                 cGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGt                                                                                                 <  5:190671/135‑1 (MQ=255)
                                 cGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGt                                                                                                 <  2:68860/135‑1 (MQ=255)
                                  gAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTg                                                                                                >  3:233164/1‑135 (MQ=255)
                                      tGGAGTGGCCGCTGGCGATCCTGCTGGCCATCTTCTGGGTTACCTACGGTCTTGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGTTTATACGGTGCGTTCATCGTGGTCa                                                                                            <  5:290894/135‑1 (MQ=255)
                                      tGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGTCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCa                                                                                            >  2:251550/1‑135 (MQ=255)
                                      tGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCa                                                                                            <  7:23791/135‑1 (MQ=255)
                                      tGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCa                                                                                            <  5:21038/135‑1 (MQ=255)
                                      tGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCa                                                                                            >  8:237271/1‑135 (MQ=255)
                                         aGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCAcc                                                                                          <  4:82730/134‑1 (MQ=255)
                                                gCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGAtgct                                                                                  >  8:50488/1‑135 (MQ=255)
                                                gCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGAtgct                                                                                  >  4:25183/1‑135 (MQ=255)
                                                      gATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAACCGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACAt                                                                            <  1:292580/135‑1 (MQ=255)
                                                         ccTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGt                                                                         <  2:213971/135‑1 (MQ=255)
                                                              tgGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCaa                                                                      <  1:251550/133‑1 (MQ=255)
                                                                   aTCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACAtc                                                               >  4:298453/1‑135 (MQ=255)
                                                                     cGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACAtctc                                                             <  5:83595/135‑1 (MQ=255)
                                                                                ccTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTg                                                   <  3:121501/134‑1 (MQ=255)
                                                                                   aCGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGcc                                               >  4:135171/1‑135 (MQ=255)
                                                                                         tGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGAtgc                                                                                   >  7:253747/1‑93 (MQ=255)
                                                                                         tGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGAtgc                                                                                   <  8:253747/93‑1 (MQ=255)
                                                                                         tGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCa                                         <  3:298453/135‑1 (MQ=255)
                                                                                             gtTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                                                                                                                           >  7:210201/1‑49 (MQ=255)
                                                                                             gtTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGc                                                                                                                           <  8:210201/49‑1 (MQ=255)
                                                                                             gtTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTc                                     >  7:105110/1‑135 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCGCATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTg                                                   >  8:182668/1‑120 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCGCATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTg                                                   >  8:187297/1‑120 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCGCATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTg                                                   <  7:187297/120‑1 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCGCATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTg                                                   <  7:182668/120‑1 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTc                                                                                                                  <  5:221295/57‑1 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTc                                                                                                                  >  6:221295/1‑57 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCt                                                                                                                 <  6:118008/58‑1 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCt                                                                                                                 >  5:118008/1‑58 (MQ=255)
                                                                                              ttcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTc                                     >  6:222544/1‑134 (MQ=255)
                                                                                              ttcttcGGCACTAGCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCg                                    <  8:182238/135‑1 (MQ=255)
                                                                                               tcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCg                                    <  6:125472/134‑1 (MQ=255)
                                                                                               tcttcGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGt                                   <  3:25183/135‑1 (MQ=255)
                                                                                                        cTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCt                          <  7:37623/135‑1 (MQ=255)
                                                                                                                     cAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCAc             <  1:206196/135‑1 (MQ=255)
                                                                                                                        gACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGa          >  7:5319/1‑135 (MQ=255)
                                                                                                                         aCCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGAc         <  8:115456/135‑1 (MQ=255)
                                                                                                                          ccAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGAcg        <  2:85464/135‑1 (MQ=255)
                                                                                                                          ccAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGAcg        <  3:104012/135‑1 (MQ=255)
                                                                                                                           cAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGAcg        >  5:54518/1‑134 (MQ=255)
                                                                                                                               cacaTCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCgatg   <  7:237271/135‑1 (MQ=255)
                                                                                                                               cacaTCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCgatg   <  3:135171/135‑1 (MQ=255)
                                                                                                                                acaTCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCgatga  >  4:158051/1‑135 (MQ=255)
                                                                                                                                acaTCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCgatga  >  3:169243/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                                      
GCCGATGGGCTACACCACCACCAAGGAATACGCCGAGCTGGAGTGGCCGCTGGCGATCCTGCTGGCCATCGTCTGGGTTACCTACGGGCTGGTGTTCTTCGGCACTATCGTCAAACGCAAGACCAAGCACATCTATGTCGGCAACTGGTTCTACGGTGCGTTCATCGTGGTCACCGCGATGCTGCACATCGTCAACCACATCTCGCTGCCGGTGAGCCTGTTCAAGTCGTACTCGGCCTACTCCGGCGCCACCGACGCGATGA  >  ACB122/4833941‑4834203

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