Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,837,333 C→A 6.4% 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,3330CA 6.4% 300.7 / 6.2 93P432Q (CCA→CAA) ccoN‑Icbb3‑type cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (47/40);  new base A (2/4);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCT  >  NC_002947/4837186‑4837475
                                                                                                                                                   |                                                                                                                                              
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)
                          aCCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCggcggc                                                                                                                    <  2:985622/150‑1 (MQ=255)
                            cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                  >  2:2369198/1‑150 (MQ=255)
                            cGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGgcgc                                                                                                                  >  2:1206562/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                                                                                          >  1:982667/1‑150 (MQ=255)
                                                    ggTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCt                                                                                          >  2:1766972/1‑150 (MQ=255)
                                                     gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcg                                                                                                                                        >  2:683929/1‑103 (MQ=255)
                                                     gTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTtcg                                                                                                                                        <  1:683929/103‑1 (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:1421759/150‑1 (MQ=255)
                                                                cGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                              <  1:1462171/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:456382/150‑1 (MQ=255)
                                                                        gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                      <  1:131454/150‑1 (MQ=255)
                                                                        gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                      >  2:1422914/1‑150 (MQ=255)
                                                                        gCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGgc                                                                      <  1:2085144/150‑1 (MQ=255)
                                                                           tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                     >  2:21280/1‑116 (MQ=255)
                                                                           tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGc                                                                                                     <  1:21280/116‑1 (MQ=255)
                                                                           tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                              <  2:450758/139‑1 (MQ=255)
                                                                           tGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGaa                                                                              >  1:450758/1‑139 (MQ=255)
                                                                                  gcgcgcGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                                 >  1:980194/1‑113 (MQ=255)
                                                                                  gcgcgcGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCa                                                                                                 <  2:980194/113‑1 (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    <  1:2380583/150‑1 (MQ=255)
                                                                                                                                          ccagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGc    <  2:951022/150‑1 (MQ=255)
                                                                                                                                           cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                <  2:2157191/137‑1 (MQ=255)
                                                                                                                                           cagccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCtt                >  1:2157191/1‑137 (MQ=255)
                                                                                                                                            agccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCt  >  1:251361/1‑150 (MQ=255)
                                                                                                                                            agccaCCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCt  <  2:1518539/150‑1 (MQ=255)
                                                                                                                                                  ccAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTgg                                                                       <  2:136099/75‑1 (MQ=255)
                                                                                                                                                  ccAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTgg                                                                       >  1:136099/1‑75 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGCACTTCTGGCTGGCCACCATCGGTACCGTGCTGTACATCGCCTCGATGTGGGTCAACGGCATCGCCCAGGGCCTGATGTGGCGCGCGGTCAACAGCGACGGCACGCTCACCTACTCGTTCGTGGAAACCTTGGTGGCCAGCCACCCAGGCTTCGTCGTGCGCTTCGTCGGCGGCGCGATCTTCCTCAGCGGCATGTTCCTGATGGCCTGGAACACCTGGCGCACCGTGCGTTCGCCAGCGCTTGACGTCGCCCCTGCGAACGCCCAGTTGGCTTGAGGAGAAACGCCT  >  NC_002947/4837186‑4837475

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