Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,847,822 A→G 5.2% H317R (CAC→CGC) ‡ ccoP‑II → cbb3‑type cytochrome c oxidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,847,8220AG 5.2% 621.5 / 14.1 191H317R (CAC→CGC) ‡ccoP‑IIcbb3‑type cytochrome c oxidase subunit
Reads supporting (aligned to +/- strand):  ref base A (44/136);  new base G (3/7);  total (47/144)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 196 total aligned reads displayed.

CCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAA  >  NC_002947/4847673‑4847971
                                                                                                                                                     |                                                                                                                                                      
  acccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCAcc                                                                                                                                                      >  1:4622117/1‑150 (MQ=255)
  acccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCAcc                                                                                                                                                      >  1:4462167/1‑150 (MQ=255)
    ccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAg                                                                                                                                                    <  2:923814/150‑1 (MQ=255)
     cGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAgg                                                                                                                                                   <  2:4807183/150‑1 (MQ=255)
       cAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCt                                                                                                                                                 <  1:429902/150‑1 (MQ=255)
         ggCTTTCATCTACGGCTCGAGCTTTGCCGAGTTGCATCAGACCATTAGTTATTGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGTGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCTgt                                                                                                                                               <  1:80229/150‑2 (MQ=255)
         ggCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGa                                                                                                                                               <  2:5129379/150‑1 (MQ=255)
          gCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCTgta                                                                                                                                              <  1:4339525/150‑3 (MQ=255)
            gTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGcc                                                                                                                                            >  1:1958858/1‑150 (MQ=255)
            gTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCTgtaga                                                                                                                                            <  1:1686270/150‑5 (MQ=21)
              tCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCgg                                                                                                                                          <  1:1674615/150‑1 (MQ=255)
              tCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCgg                                                                                                                                          <  2:3682301/150‑1 (MQ=255)
                 tCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCa                                                                                                                                       >  2:3713464/1‑150 (MQ=255)
                  cTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCaa                                                                                                                                      <  1:2436108/150‑1 (MQ=255)
                    aCGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAgg                                                                                                                                    <  2:89977/150‑1 (MQ=255)
                      ggCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCt                                                                                                                                  <  2:1072792/150‑1 (MQ=255)
                      ggCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCt                                                                                                                                  <  2:3991087/150‑1 (MQ=255)
                      ggCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCt                                                                                                                                  <  1:3713464/150‑1 (MQ=255)
                      ggCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCt                                                                                                                                  <  1:4751832/150‑1 (MQ=255)
                        cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                                <  1:2245577/150‑1 (MQ=255)
                        cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                                <  1:5074554/150‑1 (MQ=255)
                         tCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAg                                                                                                                               >  1:1528566/1‑150 (MQ=255)
                           gAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTa                                                                                                                             >  2:2576093/1‑150 (MQ=255)
                            aGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTaa                                                                                                                            <  2:1041639/150‑1 (MQ=255)
                            aGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTaa                                                                                                                            <  2:4491208/150‑1 (MQ=255)
                            aGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTaa                                                                                                                            >  2:2074518/1‑150 (MQ=255)
                               tttGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggc                                                                                                                         >  1:5348397/1‑150 (MQ=255)
                               tttGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggc                                                                                                                         >  1:2524800/1‑150 (MQ=255)
                                 tGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggcag                                                                                                                       <  1:756597/150‑1 (MQ=255)
                                   cccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggcaggc                                                                                                                     >  2:391/1‑150 (MQ=255)
                                   cccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggcaggc                                                                                                                     <  1:1429670/150‑1 (MQ=255)
                                    ccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcc                                                                                                                    >  2:1855023/1‑150 (MQ=255)
                                     cAGTTGCAGCAGTCCATTAGTTATGGCCGCCAGGGGCAGATGCCGGACCAGGAACATCTGCAGGGCATGGACATGGTGCACTTGGTGGCGGCTTAAGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGATGAGTAAGGCAggcct                                                                                                                   <  1:4536592/150‑1 (MQ=255)
                                      aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                                  <  2:4609965/150‑1 (MQ=255)
                                      aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                                  <  2:2015891/150‑1 (MQ=255)
                                      aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                                  <  1:5478394/150‑1 (MQ=255)
                                      aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                                  <  1:3038563/150‑1 (MQ=255)
                                        ttGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggc                                                                                                                <  1:2074518/150‑1 (MQ=255)
                                        ttGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggc                                                                                                                <  2:3128015/150‑1 (MQ=255)
                                        ttGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggc                                                                                                                <  2:5528830/150‑1 (MQ=255)
                                        ttGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggc                                                                                                                <  1:1045308/150‑1 (MQ=255)
                                           cagcagACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCctc                                                                                                             <  1:2004261/150‑1 (MQ=255)
                                              cagACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTc                                                                                                          >  2:272583/1‑150 (MQ=255)
                                                 aCCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTcgcg                                                                                                       >  2:2245550/1‑150 (MQ=255)
                                                   cATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCggg                                                                                                     >  2:118104/1‑150 (MQ=255)
                                                       cGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTaaa                                                                                                 <  1:953585/150‑1 (MQ=255)
                                                       cGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTaaa                                                                                                 >  2:2578865/1‑150 (MQ=255)
                                                        gTTATGGCCTCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGAAACGACAAGGTGCACTTGCTGGCGTCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAc                                                                                                <  2:1467007/150‑1 (MQ=255)
                                                         ttATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAcc                                                                                               <  1:3553200/150‑1 (MQ=255)
                                                           aTGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCg                                                                                             <  1:4795307/150‑1 (MQ=255)
                                                            tGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGc                                                                                            <  1:3299378/150‑1 (MQ=255)
                                                            tGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGc                                                                                            <  1:2794532/150‑1 (MQ=255)
                                                              gccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTc                                                                                          <  1:5158691/150‑1 (MQ=255)
                                                               ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                         <  2:5330933/150‑4 (MQ=255)
                                                               ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                         <  2:1662329/150‑4 (MQ=255)
                                                               ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                         <  2:4199086/150‑1 (MQ=255)
                                                                 gccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaaga                                                                                       <  2:4462167/150‑6 (MQ=255)
                                                                    aGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                    <  2:5348397/150‑1 (MQ=255)
                                                                     ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                   <  2:2970069/150‑1 (MQ=255)
                                                                     ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCTAaaca                                                                                   >  1:3490095/1‑150 (MQ=255)
                                                                           gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                             <  2:3847041/150‑1 (MQ=255)
                                                                           gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCt                                                                            <  2:678928/150‑1 (MQ=255)
                                                                             tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                          <  1:2576093/150‑1 (MQ=255)
                                                                             tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgt                                                                          <  2:795331/150‑3 (MQ=255)
                                                                              gCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtg                                                                         <  2:1056611/150‑4 (MQ=255)
                                                                                 ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                      <  2:4432874/150‑1 (MQ=255)
                                                                                   cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                    <  2:3291184/150‑1 (MQ=255)
                                                                                   cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                    <  2:1108091/150‑1 (MQ=255)
                                                                                   cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                    <  1:423396/150‑1 (MQ=255)
                                                                                    ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                    <  2:1645346/134‑1 (MQ=255)
                                                                                    ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGgaggag                                                                   >  2:1911723/1‑150 (MQ=255)
                                                                                    ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAg                                                                   <  2:3307626/150‑1 (MQ=255)
                                                                                     cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                  <  1:1826117/150‑1 (MQ=255)
                                                                                     cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCAACACAGGATCCGTGTAGGAGc                                                                  >  1:695918/1‑150 (MQ=255)
                                                                                       ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                <  1:3958773/150‑1 (MQ=255)
                                                                                       ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                <  1:3837343/150‑1 (MQ=255)
                                                                                       ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                <  1:2691616/150‑1 (MQ=255)
                                                                                       ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                <  2:1389885/150‑1 (MQ=255)
                                                                                       ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                <  2:1390454/150‑1 (MQ=255)
                                                                                        gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                               <  2:3405330/150‑1 (MQ=255)
                                                                                        gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                               <  2:593166/150‑1 (MQ=255)
                                                                                        gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                               <  1:3122076/150‑1 (MQ=255)
                                                                                         aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGt                                                              <  1:2160470/150‑1 (MQ=255)
                                                                                         aacaacGGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCGTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCCTGGCCTCCTCGCGGGTACACCCGCTCCTACCCAGGATCCTGGTAGGAGCGGGt                                                              <  2:1238402/150‑1 (MQ=255)
                                                                                            aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTa                                                           <  1:1084043/150‑1 (MQ=255)
                                                                                            aacTGCAGGGCAACGACAAGGTGCACTTGCAGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCACCAACACAAGAAACAAGTAGGAGCCGGTTTa                                                           >  1:5015574/1‑150 (MQ=255)
                                                                                                  aGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTAGCGGGGTAAAACGCGCCTACACAGGTTGCGTGTTTGAGGGGGTTTAACcgcg                                                     >  1:605801/1‑150 (MQ=255)
                                                                                                       aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAagag                                                >  2:4137522/1‑150 (MQ=255)
                                                                                                       aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTGGTAGGAGCGGGTTTACCCGCGAagag                                                >  2:3171635/1‑150 (MQ=255)
                                                                                                           aCAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                  >  2:4552865/1‑113 (MQ=255)
                                                                                                           aCAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                  <  1:4552865/113‑1 (MQ=255)
                                                                                                                         gtGGCGGCTTACGTCTACCGCCTGTCGCACCTGGCTGAGCCGGCCAAGGCTGCGTACGGCCGGCCTGGCATCTGCGCGGGTAAACCCGCTCCTACACCGGATCCTGGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                              <  2:5015574/149‑1 (MQ=37)
                                                                                                                          tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGTACTGGCACAGCAt                             >  2:2505155/1‑150 (MQ=255)
                                                                                                                          tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCGACACAGGATACGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAt                             >  2:2504864/1‑150 (MQ=255)
                                                                                                                                       cTACTGCGTGTCGCACCCGGCTGAGCCGGCCACGGCTGTGTACGGCCGGCCGGGCCTCTTCGGGGGAAACCCCGCTCCTACACTGTATCCTTGTCGGCGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa                <  1:4137522/150‑1 (MQ=25)
                                                                                                                                        ctcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAAACCGCTCCGACCCACGATACGTGTGGGAGCGGGTTTACCCGCGAAGAGGGCGGTGCTGGGACAGCATTTATCGCAGGCCAg               >  1:4020143/3‑150 (MQ=25)
                                                                                                                                         tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAAACCGCTCCTAAACAGGATGCTTGTTGGTGCGGGTTTTACCGCGAAGAGGGCGGTACTGGGACAGCCTTTATAGCAGGGAAGt              >  1:1257245/2‑150 (MQ=255)
                                                                                                                                         tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGGAAAACCCCTCCGACACAGGAAACAGGTTGGTGCGGGTTTAACCCCCCACAGGCCCGTAAAGGCACACCAATTCTCCCACGCAAGt              >  1:3859227/2‑150 (MQ=255)
                                                                                                                                                     aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCaa  >  1:1847085/1‑150 (MQ=255)
                                                                                                                                                     aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCAACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCaa  >  2:3963758/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAA  >  NC_002947/4847673‑4847971

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