Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,280 C→T 5.3% A429A (GCC→GCT 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,2800CT 5.3% 292.1 / 2.6 93A429A (GCC→GCTccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (42/46);  new base T (3/2);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACG  >  NC_002947/4837132‑4837427
                                                                                                                                                    |                                                                                                                                                   
tGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcc                                                                                                                                                     >  8:269927/1‑149 (MQ=255)
    cAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgccagcc                                                                                                                                                 >  5:134627/1‑149 (MQ=255)
       tGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCAcc                                                                                                                                              >  2:196171/1‑149 (MQ=255)
         cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCAccc                                                                                                                                             <  1:40615/148‑1 (MQ=255)
         cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCAcc                                                                                                                                              >  5:170996/1‑147 (MQ=255)
         cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCa                                                                                                                                            <  7:199290/149‑1 (MQ=255)
          gcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAg                                                                                                                                           <  3:51209/149‑1 (MQ=255)
                ctggctgtCCAGCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGCGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtc                                                                                                                                      <  2:230156/148‑1 (MQ=255)
                ctggctggCCACCATCGTTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGCGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtc                                                                                                                                      <  2:224027/148‑1 (MQ=255)
                   gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAGCCTGGGTGGCTGGCCACCAAGGTTtccac                                                                                                                                   <  5:208170/148‑4 (MQ=18)
                   gctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgt                                                                                                                                    <  2:138288/147‑1 (MQ=255)
                       gcCACCATCGGTAGCGTTCGGTAGATAGCCTCGCAGTGGGTCAACGGCCTCGCCCAGGGCCTGATGTGTCGCGCGGGCAACAGCGACGGCACGCTAACCAACTCCTTCGTGGAAACCTTTGTGGCCAGCCAACCACGCTTCGTCGTgc                                                                                                                               <  7:169269/148‑1 (MQ=255)
                       gcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgc                                                                                                                               <  8:279462/148‑1 (MQ=255)
                         caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCt                                                                                                                            >  7:34476/1‑149 (MQ=255)
                         caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCt                                                                                                                            >  5:333007/1‑149 (MQ=255)
                                 gTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCCACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtc                                                                                                                                   <  1:392129/134‑1 (MQ=255)
                                 gTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCCACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcgtc                                                                                                                                   >  2:392129/1‑134 (MQ=255)
                                     cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtcgtc                                                                                                                       >  5:254661/1‑142 (MQ=255)
                                     cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtcgtc                                                                                                                       <  6:254661/142‑1 (MQ=255)
                                     cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGgcggcg                                                                                                                >  4:163131/1‑149 (MQ=255)
                                      gTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                               <  6:134627/149‑1 (MQ=255)
                                          tGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATc                                                                                                           >  4:95319/1‑149 (MQ=255)
                                          tGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATc                                                                                                           >  4:170820/1‑149 (MQ=255)
                                            tACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCtt                                                                                                         <  1:196171/149‑1 (MQ=255)
                                             aCATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                        >  2:22627/1‑149 (MQ=255)
                                              cATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                        <  3:163131/148‑1 (MQ=255)
                                               aTCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                      >  7:163282/1‑149 (MQ=255)
                                                tCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                      >  6:158572/1‑148 (MQ=255)
                                                  gCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAg                                                                                                   >  1:181140/1‑149 (MQ=255)
                                                   ccTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                  >  8:47927/1‑149 (MQ=255)
                                                     tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgc                                                                                                                               <  3:215253/118‑1 (MQ=255)
                                                     tCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgc                                                                                                                               >  4:215253/1‑118 (MQ=255)
                                                      cGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgcgc                                                                                                                             <  2:205267/119‑1 (MQ=255)
                                                      cGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgcgc                                                                                                                             >  1:205267/1‑119 (MQ=255)
                                                       gATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCgg                                                                                                <  3:95319/147‑1 (MQ=255)
                                                             ggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTcc                                                                                        >  6:1173/1‑149 (MQ=255)
                                                              gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                       <  6:333007/149‑1 (MQ=255)
                                                              gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                       <  3:170820/149‑1 (MQ=255)
                                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                       >  1:128685/1‑148 (MQ=255)
                                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTg                                                                                      <  7:47927/149‑1 (MQ=255)
                                                                cAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGa                                                                                     >  1:293062/1‑149 (MQ=255)
                                                                  aCGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATg                                                                                   >  7:31194/1‑149 (MQ=255)
                                                                  aCGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATg                                                                                   >  1:215692/1‑149 (MQ=255)
                                                                   cGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATgg                                                                                  <  8:31194/149‑1 (MQ=255)
                                                                   cGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATgg                                                                                  >  7:149510/1‑149 (MQ=255)
                                                                    ggCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGc                                                                                 <  4:92293/149‑1 (MQ=255)
                                                                       aTCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAACGGCATGTTCCTGATGGCCTg                                                                              >  2:179761/1‑149 (MQ=255)
                                                                        tCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTgg                                                                             <  2:293062/149‑1 (MQ=255)
                                                                          gCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACTCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTTATGGCCTGGaa                                                                           <  4:225887/149‑1 (MQ=255)
                                                                          gCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                           >  7:268963/1‑149 (MQ=255)
                                                                           cccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                           >  6:292808/1‑148 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtc                                                                                                                          <  2:30089/100‑1 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtc                                                                                                                          >  1:30089/1‑100 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAaca                                                                         >  6:361850/1‑149 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAaca                                                                         >  1:25555/1‑149 (MQ=255)
                                                                             cAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAacac                                                                        >  3:161858/1‑149 (MQ=255)
                                                                             cAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAacac                                                                        <  5:361850/149‑1 (MQ=255)
                                                                               gggCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCt                                                                      <  2:181140/149‑1 (MQ=255)
                                                                                ggCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTg                                                                     >  7:229539/1‑149 (MQ=255)
                                                                                  ccTGATGTGGCGAGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                   >  4:198799/1‑149 (MQ=255)
                                                                                     gATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCa                                                                >  2:45225/1‑149 (MQ=255)
                                                                                         tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                      <  1:274821/107‑1 (MQ=255)
                                                                                         tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                      >  2:274821/1‑107 (MQ=255)
                                                                                         tgGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCg                                                             <  2:215692/148‑1 (MQ=255)
                                                                                             gcgcgGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGt                                                        <  6:170996/149‑1 (MQ=255)
                                                                                                   tCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGcag                                                  >  4:348895/1‑147 (MQ=255)
                                                                                                   tCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGc                                                    <  5:292808/147‑1 (MQ=255)
                                                                                                     aaCAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAgc                                                <  4:161858/149‑1 (MQ=255)
                                                                                                           gACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGa                                          >  7:231803/1‑149 (MQ=255)
                                                                                                             cGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACg                                        <  7:269927/149‑1 (MQ=255)
                                                                                                                 aCGCTCACCTACTCGTTCGTGGAAACGTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGt                                                                                           <  1:148685/94‑1 (MQ=255)
                                                                                                                 aCGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGt                                                                                           >  2:148685/1‑94 (MQ=255)
                                                                                                                     tCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCt                                >  3:212073/1‑149 (MQ=255)
                                                                                                                       aCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCTCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCt                                <  4:370906/147‑1 (MQ=255)
                                                                                                                       aCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCt                                <  4:376890/147‑1 (MQ=255)
                                                                                                                       aCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGc                              <  8:268963/149‑1 (MQ=255)
                                                                                                                            cTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAAc                          <  8:34476/148‑1 (MQ=255)
                                                                                                                               gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGc                                                    <  8:186566/119‑1 (MQ=255)
                                                                                                                               gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGa                                          >  4:215794/1‑129 (MQ=255)
                                                                                                                               gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGa                                          <  3:215794/129‑1 (MQ=255)
                                                                                                                                 acgacgAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAg                    >  3:319586/6‑149 (MQ=38)
                                                                                                                                   gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                  <  8:229539/149‑1 (MQ=255)
                                                                                                                                   gacgAAGCCTGGGTGGCTGGCCACCAAGGTTTCGTCGTGCGCTTCGTGGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACATGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                  >  6:341683/4‑149 (MQ=255)
                                                                                                                                     ggAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                >  6:204252/1‑149 (MQ=255)
                                                                                                                                      gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                <  5:1173/148‑1 (MQ=255)
                                                                                                                                        aaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt             >  3:233338/1‑149 (MQ=255)
                                                                                                                                           ctgggtggCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag          <  3:78770/146‑1 (MQ=255)
                                                                                                                                           ctgggtggCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag          >  4:78770/4‑149 (MQ=255)
                                                                                                                                           cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                <  3:349672/143‑1 (MQ=255)
                                                                                                                                           cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTgg                >  4:349672/1‑143 (MQ=255)
                                                                                                                                           cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag          <  2:128685/149‑1 (MQ=255)
                                                                                                                                           cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag          <  1:170052/149‑1 (MQ=255)
                                                                                                                                             ttgtggCCAGCGACCTTGGTTTGGTCGTTCTGTTCGTCGGCGGCGCGATGTTCCTCAGCGGGATGTTCCTGTTGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTga           <  8:231803/144‑1 (MQ=255)
                                                                                                                                              ggtggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagg         <  2:25555/147‑1 (MQ=255)
                                                                                                                                                tggCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGgaga      >  2:388146/1‑148 (MQ=255)
                                                                                                                                                  gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAAc   <  5:204252/149‑1 (MQ=255)
                                                                                                                                                   ccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACg  >  2:147440/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACG  >  NC_002947/4837132‑4837427

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: