Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,325 C→T 8.0% 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,3250CT 8.0% 242.8 / 9.3 87A429A (GCC→GCTccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (45/35);  new base T (2/5);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGC  >  NC_002947/4837177‑4837473
                                                                                                                                                    |                                                                                                                                                    
tGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                                     <  2:821387/150‑1 (MQ=255)
 gATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACACCGCCTCGGTGTGGGTCAACGGCATCTCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTGGTTCGTGGAAACCTTGGTGgccag                                                                                                                                                    >  1:553184/1‑150 (MQ=255)
     aaTGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCAc                                                                                                                                                >  1:1386130/1‑150 (MQ=255)
         cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAg                                                                                                                                            <  2:201046/150‑1 (MQ=255)
         cgcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAg                                                                                                                                            <  2:893985/150‑1 (MQ=255)
          gcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAgg                                                                                                                                           <  2:1959198/150‑1 (MQ=255)
          gcACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAgg                                                                                                                                           >  1:672203/1‑150 (MQ=255)
           cACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCTGGCCACCAAggt                                                                                                                                          <  1:2286235/150‑2 (MQ=255)
                ctggctggCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAGCCTGGGTGGCTGGCCACCAAGGTTtcca                                                                                                                                     <  2:362439/150‑3 (MQ=18)
                    ctggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTg                                                                                                                                 <  2:209362/150‑1 (MQ=255)
                     tggcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTgc                                                                                                                                >  1:1077271/1‑150 (MQ=255)
                       gcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTcgtcgt                                                                                                                                  >  2:1875645/1‑146 (MQ=255)
                       gcCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTcgtcgt                                                                                                                                  <  1:1875645/146‑1 (MQ=255)
                         caccaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCtt                                                                                                                            <  2:465515/150‑1 (MQ=255)
                          accaTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTCCTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtc                                                                                                                           <  2:757827/150‑1 (MQ=255)
                             aTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtcgtc                                                                                                                        >  1:951022/1‑150 (MQ=255)
                             aTCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTtcgtc                                                                                                                        >  1:688853/1‑150 (MQ=255)
                               cGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                                     >  1:1374888/1‑119 (MQ=255)
                               cGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                                     <  2:1374888/119‑1 (MQ=255)
                                  tACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                                     >  2:727224/1‑116 (MQ=255)
                                  tACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGgcca                                                                                                                                                     <  1:727224/116‑1 (MQ=255)
                                   aCCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCggcggc                                                                                                                  <  2:985622/150‑1 (MQ=255)
                                     cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                >  2:1206562/1‑150 (MQ=255)
                                     cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                >  2:2369198/1‑150 (MQ=255)
                                     cGTGCTGTACATCGCCAAGATGTGGGAGCACAGCAACGCACAGGGTCTGTGATGGTGCGCCGTCAGCAGCGACGTCACGCTCACCTGCTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                >  2:1481464/1‑150 (MQ=255)
                                      gTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGcgcg                                                                                                               >  2:1519125/1‑150 (MQ=255)
                                          tGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCt                                                                                                           >  2:2085144/1‑150 (MQ=255)
                                            tACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTc                                                                                                         >  2:1421759/1‑150 (MQ=255)
                                              cATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCt                                                                                                       <  1:1705282/150‑1 (MQ=255)
                                                tCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGGGGCGCGATCTTCCTCa                                                                                                     >  2:131454/1‑150 (MQ=255)
                                                  gCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTctga                                                                                                   <  2:1077271/150‑4 (MQ=255)
                                                  gCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                   >  1:456382/1‑150 (MQ=255)
                                                   ccTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCg                                                                                                  >  2:1633740/1‑150 (MQ=255)
                                                         tgtgGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGt                                                                                            >  2:2380583/1‑150 (MQ=255)
                                                             ggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                        >  2:1766972/1‑150 (MQ=255)
                                                             ggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                        >  1:982667/1‑150 (MQ=255)
                                                              gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcg                                                                                                                                      <  1:683929/103‑1 (MQ=255)
                                                              gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcg                                                                                                                                      >  2:683929/1‑103 (MQ=255)
                                                               tCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACACAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGa                                                                                      >  2:574919/1‑150 (MQ=255)
                                                               tCAACGGCATCACCCAGATCCTGCTGTGGCGCGCGGTCAACAGCTACGGCACGCTCACCTACTCGTTCGTTGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCATGCGCTTCGTAGGCGGCTCGATCTTGCTCAGCGGCATGCTCCTgc                                                                                      >  2:859094/1‑149 (MQ=255)
                                                                 aaCGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATg                                                                                    >  1:1451134/1‑150 (MQ=255)
                                                                   cGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCCGCGGCATGTTCCTGATGGc                                                                                  >  2:1952275/1‑150 (MQ=255)
                                                                   cGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGc                                                                                  >  2:2179631/1‑150 (MQ=255)
                                                                         cGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                            <  1:1462171/150‑1 (MQ=255)
                                                                         cGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                            <  1:1421759/150‑1 (MQ=255)
                                                                           cccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTgg                                                                              <  2:1929116/145‑1 (MQ=255)
                                                                           cccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTgg                                                                              >  1:1929116/1‑145 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                <  1:967561/111‑1 (MQ=255)
                                                                            ccAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                >  2:967561/1‑111 (MQ=255)
                                                                                 gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                    >  2:1422914/1‑150 (MQ=255)
                                                                                 gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                    <  1:131454/150‑1 (MQ=255)
                                                                                 gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                    <  2:456382/150‑1 (MQ=255)
                                                                                 gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                    <  1:2085144/150‑1 (MQ=255)
                                                                                    tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                   <  1:21280/116‑1 (MQ=255)
                                                                                    tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                   >  2:21280/1‑116 (MQ=255)
                                                                                    tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                            >  1:450758/1‑139 (MQ=255)
                                                                                    tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                            <  2:450758/139‑1 (MQ=255)
                                                                                           gcgcgcGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                               <  2:980194/113‑1 (MQ=255)
                                                                                           gcgcgcGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                               >  1:980194/1‑113 (MQ=255)
                                                                                             gcgcgGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGtt                                                        >  2:610133/1‑150 (MQ=255)
                                                                                                  gTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCa                                                   <  2:925537/150‑1 (MQ=255)
                                                                                                   tCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGc                                                     <  2:1965202/147‑1 (MQ=255)
                                                                                                   tCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGc                                                     >  1:1965202/1‑147 (MQ=255)
                                                                                                   tCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAg                                                  >  2:586452/1‑150 (MQ=255)
                                                                                                     aaCAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCCCCGTGCGTTCGCCAgcg                                                >  2:418684/1‑150 (MQ=255)
                                                                                                         gCGACGGCACGCTCACCTACTCGTTCGTGTAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTg                                            >  2:1354013/1‑150 (MQ=255)
                                                                                                             cGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGt                                        <  2:1386130/150‑1 (MQ=255)
                                                                                                                  cGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGccc                                   <  1:934951/150‑1 (MQ=255)
                                                                                                                   gCTCCCCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGcccc                                  >  2:2931/1‑150 (MQ=255)
                                                                                                                     tCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTg                                >  2:479885/1‑150 (MQ=255)
                                                                                                                     tCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTg                                >  1:218598/1‑150 (MQ=255)
                                                                                                                         cTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGAt                                                                                     >  1:1111344/1‑93 (MQ=255)
                                                                                                                         cTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGAt                                                                                     <  2:1111344/93‑1 (MQ=255)
                                                                                                                         cTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGaa                            >  2:228744/1‑150 (MQ=255)
                                                                                                                          tACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAAc                           <  2:982667/150‑1 (MQ=255)
                                                                                                                          tACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAAc                           <  1:1766972/150‑1 (MQ=255)
                                                                                                                           aCTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACg                          >  2:2235990/1‑150 (MQ=255)
                                                                                                                             tcgttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGcc                        <  1:586452/150‑1 (MQ=255)
                                                                                                                              cgttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGccc                       <  1:1206562/150‑1 (MQ=255)
                                                                                                                              cgttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGccc                       <  1:1481464/150‑1 (MQ=255)
                                                                                                                               gttcgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCa                      <  2:672203/150‑1 (MQ=255)
                                                                                                                                  cgtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGtt                   <  2:688853/150‑1 (MQ=255)
                                                                                                                                   gtGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTg                  >  2:2370235/1‑150 (MQ=255)
                                                                                                                                      gAAGCCTGGGTGGCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCt               >  1:1068110/1‑150 (MQ=255)
                                                                                                                                      gAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTGCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCt               <  1:866467/150‑1 (MQ=255)
                                                                                                                                       aaaCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt              <  2:1451134/150‑1 (MQ=255)
                                                                                                                                         gccTGGGTGGCTGGCCACCAAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTga            >  2:742096/2‑150 (MQ=255)
                                                                                                                                         aCCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTga            <  1:610133/150‑1 (MQ=255)
                                                                                                                                           ctgggtggCTGGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgag           <  2:170884/146‑1 (MQ=255)
                                                                                                                                           cTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagg          >  1:2003082/1‑150 (MQ=255)
                                                                                                                                            ttGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagga         >  1:1702419/1‑150 (MQ=255)
                                                                                                                                            ttGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTgagga         >  2:2354902/1‑150 (MQ=255)
                                                                                                                                                 ggCCAGCCACCCAGGCTTCGTAGTGCGCTTCGTCGGCGGCGCGATTTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAAc    <  2:553184/150‑1 (MQ=255)
                                                                                                                                                  gccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACg   >  1:18037/1‑150 (MQ=255)
                                                                                                                                                   ccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGc  <  2:951022/150‑1 (MQ=255)
                                                                                                                                                   ccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGc  <  1:2380583/150‑1 (MQ=255)
                                                                                                                                                    cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt              <  2:2157191/137‑1 (MQ=255)
                                                                                                                                                    cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt              >  1:2157191/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGATCAATGCGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGC  >  NC_002947/4837177‑4837473

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: