Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,288 C→A 10.3% P432Q (CCA→CAA)  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,837,2880CA10.3% 202.5 / 10.9 68P432Q (CCA→CAA) ccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (36/25);  new base A (4/3);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGA  >  NC_002947/4837142‑4837435
                                                                                                                                                  |                                                                                                                                                   
gcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAg                                                                                                                                                   >  4:365586/1‑149 (MQ=255)
      ctggctggCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtc                                                                                                                                              <  4:340685/148‑1 (MQ=255)
      ctggctggCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtc                                                                                                                                              <  6:343986/148‑1 (MQ=255)
             gcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTg                                                                                                                                        <  1:150354/147‑1 (MQ=255)
               caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAg                                                                                                                                                   >  2:160465/1‑134 (MQ=255)
               caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAg                                                                                                                                                   <  1:160465/134‑1 (MQ=255)
               caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtcg                                                                                                                                          >  8:153529/1‑143 (MQ=255)
               caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtcg                                                                                                                                          <  7:153529/143‑1 (MQ=255)
               caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCt                                                                                                                                    >  8:130068/1‑149 (MQ=255)
                   aTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtcgt                                                                                                                                >  2:199039/1‑149 (MQ=255)
                        tACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGc                                                                                                                                                 <  4:94410/127‑1 (MQ=255)
                        tACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGc                                                                                                                                                 >  3:94410/1‑127 (MQ=255)
                           cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAGGTTtccac                                                                                                                                           >  1:210844/1‑127 (MQ=21)
                           cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAGGTTtccac                                                                                                                                           <  2:210844/130‑4 (MQ=21)
                            gTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                       >  6:59836/1‑149 (MQ=255)
                                    cATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                                <  5:59836/148‑1 (MQ=255)
                                         ccTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                          >  2:83515/1‑149 (MQ=255)
                                               tgtgGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATg                                                                                                    >  6:210772/1‑149 (MQ=255)
                                                   ggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTcc                                                                                                >  3:1005/1‑149 (MQ=255)
                                                    gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                               >  7:321310/1‑149 (MQ=255)
                                                     tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                               >  2:368769/1‑148 (MQ=255)
                                                        aCGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATg                                                                                           <  5:210772/149‑1 (MQ=255)
                                                         cGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATGTTCCTCAGCGGCATGTTCCTGATgg                                                                                          >  2:72245/1‑149 (MQ=255)
                                                            cATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCt                                                                                       <  1:368769/149‑1 (MQ=255)
                                                                  ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCTTAGGCGGCGCGATCTTCCTCAGCGGCATTTTCCTGATGGCAGGGAaca                                                                                 >  6:244047/1‑149 (MQ=37)
                                                                    aGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACAcc                                                                               >  4:178977/1‑149 (MQ=255)
                                                                           gATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCa                                                                        >  2:208544/1‑149 (MQ=255)
                                                                              gtgGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACATGGCGCACCg                                                                     >  5:288492/1‑149 (MQ=255)
                                                                               tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGt                                                                    >  3:18397/1‑149 (MQ=255)
                                                                                ggCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTg                                                                   >  3:233990/1‑149 (MQ=255)
                                                                                     gcgGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTc                                                              <  7:130068/149‑1 (MQ=255)
                                                                                     gcgGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTc                                                              <  2:126912/149‑1 (MQ=255)
                                                                                            aCAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                                   >  8:81939/1‑121 (MQ=255)
                                                                                            aCAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                                   <  7:81939/121‑1 (MQ=255)
                                                                                            aCAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAgcg                                                       <  1:72245/149‑1 (MQ=255)
                                                                                             cAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAgcgc                                                      >  4:257326/1‑149 (MQ=255)
                                                                                              aGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCt                                                     >  7:287062/1‑149 (MQ=255)
                                                                                                  aCGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCTGCGCGATCTTCGTCACCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCCTGCGTTCGCCAGCGCTTGAc                                                 >  7:109622/1‑149 (MQ=255)
                                                                                                  aCGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGAc                                                 <  4:1005/149‑1 (MQ=255)
                                                                                                       aCGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGc                                            <  8:321310/149‑1 (MQ=255)
                                                                                                         gCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGATCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGccc                                          >  1:92982/1‑149 (MQ=255)
                                                                                                           tCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCt                                        >  4:313910/1‑149 (MQ=255)
                                                                                                             aCCTACTCGTTCGTTGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTg                                       <  5:244047/148‑1 (MQ=255)
                                                                                                              ccTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGc                                      >  5:38538/1‑148 (MQ=255)
                                                                                                               cTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCg                                     >  1:260114/1‑148 (MQ=255)
                                                                                                                  cTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACg                                 <  3:257326/149‑1 (MQ=255)
                                                                                                                  cTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACg                                 <  4:18397/149‑1 (MQ=255)
                                                                                                                   tcgttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGc                                >  5:52575/1‑149 (MQ=255)
                                                                                                                     gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGc                                <  3:49962/147‑1 (MQ=255)
                                                                                                                      ttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                               >  6:103188/1‑147 (MQ=255)
                                                                                                                         gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                                      <  6:180848/73‑1 (MQ=255)
                                                                                                                         gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                                      >  5:180848/1‑73 (MQ=255)
                                                                                                                         gacgAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTGCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCa                             >  1:125495/4‑146 (MQ=255)
                                                                                                                         gacgAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCa                             <  2:125495/143‑1 (MQ=255)
                                                                                                                           ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCTGCATGTGCCTGATGGCCTGGAACAACTGGCGCACCGGGCGTTCGACAGCGCTTGACGTCGCCCCTGCGACCGCCCAGTTgt                        >  6:255137/1‑148 (MQ=255)
                                                                                                                           ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                        >  8:27989/1‑149 (MQ=255)
                                                                                                                            gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                        <  3:365586/148‑1 (MQ=255)
                                                                                                                            gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGc                       >  6:356752/1‑149 (MQ=255)
                                                                                                                            gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGc                       >  7:88399/1‑149 (MQ=255)
                                                                                                                            gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGc                       <  8:287062/149‑1 (MQ=255)
                                                                                                                               gccTGGGTGGCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACATCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                     >  8:207131/2‑148 (MQ=255)
                                                                                                                                 ctgggtggCTGGCCACCCAGGCTTCGGCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag                  >  1:203764/4‑149 (MQ=255)
                                                                                                                                    ggtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagga                <  1:380117/148‑1 (MQ=255)
                                                                                                                                    ggtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagg                 <  1:199039/147‑1 (MQ=255)
                                                                                                                                      tggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGaa             <  4:233990/149‑1 (MQ=255)
                                                                                                                                             ccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCtt                                                    <  2:156267/103‑1 (MQ=255)
                                                                                                                                             ccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCtt                                                    >  1:156267/1‑103 (MQ=255)
                                                                                                                                             ccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCt       >  8:295722/1‑148 (MQ=255)
                                                                                                                                                  cAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCtgatga  >  5:249243/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCTGATGA  >  NC_002947/4837142‑4837435

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: