Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,791 T→G 5.2% S316A (TCG→GCG)  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,844,7910TG 5.2% 510.5 / 14.8 153S316A (TCG→GCG) ccoP‑IIcbb3‑type cytochrome c oxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (32/113);  new base G (2/6);  total (34/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 153 total aligned reads displayed.

GACCCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTG  >  NC_002947/4844643‑4844938
                                                                                                                                                    |                                                                                                                                                     
gACCCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTc                                                                                                                                                      <  2:3195342/150‑1 (MQ=255)
gACCCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTc                                                                                                                                                      <  1:1816542/150‑1 (MQ=255)
   ccacccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCACCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCa                                                                                                                                                   >  1:2021168/1‑150 (MQ=255)
     acccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCAcc                                                                                                                                                 <  2:2927775/150‑1 (MQ=255)
     acccGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCAcc                                                                                                                                                 >  1:4222639/1‑150 (MQ=255)
        cGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAgg                                                                                                                                              <  1:2643358/150‑1 (MQ=255)
        cGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAgg                                                                                                                                              <  1:3375918/150‑1 (MQ=255)
          cAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCt                                                                                                                                            <  1:3458473/150‑1 (MQ=255)
            ggCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGa                                                                                                                                          <  1:2961821/150‑1 (MQ=255)
             gCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCTgta                                                                                                                                         <  1:1502727/150‑3 (MQ=255)
              cGTTCATCTACGGCTCGAGCTTTGCCCTGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGTGCGACAGGCTgtag                                                                                                                                        <  1:1636353/150‑4 (MQ=21)
               gTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGcc                                                                                                                                       >  2:3325576/1‑150 (MQ=255)
                ttCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCg                                                                                                                                      >  1:4369071/1‑150 (MQ=255)
                 tCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCgg                                                                                                                                     <  2:758541/150‑1 (MQ=255)
                  cATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGc                                                                                                                                    <  1:1755432/150‑1 (MQ=255)
                     cTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCaa                                                                                                                                 <  2:973222/150‑1 (MQ=255)
                      tACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAg                                                                                                                                <  1:785376/150‑1 (MQ=255)
                      tACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAg                                                                                                                                <  2:283858/150‑1 (MQ=255)
                         ggCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCt                                                                                                                             <  1:1215201/150‑1 (MQ=255)
                           cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                           <  1:3914521/150‑1 (MQ=255)
                           cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                           <  2:4308895/150‑1 (MQ=255)
                           cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATCGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                           <  2:2217197/150‑1 (MQ=255)
                           cTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATCGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGa                                                                                                                           <  2:2217159/150‑1 (MQ=255)
                            tCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACATGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAg                                                                                                                          >  2:902225/1‑150 (MQ=255)
                                gCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAag                                                                                                                      <  2:3718378/150‑1 (MQ=255)
                                   ttGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggca                                                                                                                   <  2:778735/150‑1 (MQ=255)
                                      cccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggcaggc                                                                                                                >  1:2729199/1‑150 (MQ=255)
                                      cccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAaggcaggc                                                                                                                <  2:1579005/150‑1 (MQ=255)
                                       ccAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcc                                                                                                               >  2:3266645/1‑150 (MQ=255)
                                         aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                             <  1:4432515/150‑1 (MQ=255)
                                         aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                             <  1:3699459/150‑1 (MQ=255)
                                         aGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctg                                                                                                             <  1:3301272/150‑1 (MQ=255)
                                          gTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctgg                                                                                                            <  1:2805389/150‑1 (MQ=255)
                                           ttGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggc                                                                                                           <  2:4179729/150‑1 (MQ=255)
                                             gcagcaGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAggcctggcct                                                                                                         <  2:3107777/150‑1 (MQ=255)
                                                 cagACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTc                                                                                                     <  1:1251899/150‑1 (MQ=255)
                                                 cagACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTc                                                                                                     <  2:718021/150‑1 (MQ=255)
                                                 cagACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTc                                                                                                     >  2:692983/1‑150 (MQ=255)
                                                  agACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTcg                                                                                                    >  2:4083479/1‑150 (MQ=255)
                                                   gACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTcgc                                                                                                   <  1:2042359/150‑1 (MQ=255)
                                                   gACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTcgc                                                                                                   <  1:1425817/150‑1 (MQ=255)
                                                   gACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTcgc                                                                                                   <  2:3313379/150‑1 (MQ=255)
                                                     ccATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCgg                                                                                                 <  1:2287350/150‑1 (MQ=255)
                                                     ccATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCgg                                                                                                 <  1:2441734/150‑1 (MQ=255)
                                                     ccATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCgg                                                                                                 <  1:1384576/150‑1 (MQ=255)
                                                     ccATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCgg                                                                                                 <  1:1275312/150‑1 (MQ=255)
                                                     ccATTCGGTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCgg                                                                                                 <  2:1021776/150‑1 (MQ=255)
                                                      cATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCggg                                                                                                >  1:1533179/1‑150 (MQ=255)
                                                        ttCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTa                                                                                              >  2:1702838/1‑150 (MQ=255)
                                                        ttCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTa                                                                                              >  2:3767341/1‑150 (MQ=255)
                                                           gTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAc                                                                                          <  1:938504/150‑1 (MQ=255)
                                                           gTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAc                                                                                          <  1:4072142/150‑1 (MQ=255)
                                                               tGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGc                                                                                      <  1:1147539/150‑1 (MQ=255)
                                                                ggCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCt                                                                                     <  2:540276/150‑1 (MQ=255)
                                                                  ccgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                   <  1:3461898/150‑4 (MQ=255)
                                                                   cgccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                  >  2:2764445/1‑150 (MQ=255)
                                                                    gccAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                 <  2:3183280/150‑1 (MQ=255)
                                                                       aGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                              >  2:4503916/1‑150 (MQ=255)
                                                                        ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                             <  1:966863/150‑1 (MQ=255)
                                                                        ggggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                             <  1:2986992/150‑1 (MQ=255)
                                                                                tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                    <  2:3666037/150‑1 (MQ=255)
                                                                                tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                    <  2:3130815/150‑1 (MQ=255)
                                                                                tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                    <  1:1580252/150‑1 (MQ=255)
                                                                                tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgt                                                                    <  1:2938977/150‑3 (MQ=255)
                                                                                    cgCCCAGGAACAACTGCAGGGCAACGACCAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                <  2:1338211/149‑1 (MQ=38)
                                                                                     gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAgg                                                               <  2:4396635/150‑1 (MQ=255)
                                                                                      cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                              <  2:2602730/150‑1 (MQ=255)
                                                                                      cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                              <  1:31292/150‑1 (MQ=255)
                                                                                      cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                              <  1:230378/150‑1 (MQ=255)
                                                                                       ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAg                                                             <  1:1550523/150‑1 (MQ=255)
                                                                                       ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                             <  1:1755031/150‑1 (MQ=255)
                                                                                        cAGGAACAACTGCAGGGCATCGACACGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGTAGc                                                            <  2:3009011/150‑1 (MQ=38)
                                                                                        cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGc                                                            <  1:2604673/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                          <  1:367324/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                          <  2:3314095/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                          <  1:3813021/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                          <  2:4057239/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                          <  1:4056858/150‑1 (MQ=255)
                                                                                          ggAACAACTGCAGGGCAACGACAAGGTGCACGTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                          <  1:3332764/150‑1 (MQ=255)
                                                                                           gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                         <  1:1470720/150‑1 (MQ=255)
                                                                                           gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                         <  1:1098867/150‑1 (MQ=255)
                                                                                           gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                         <  1:1098773/150‑1 (MQ=255)
                                                                                           gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                         <  1:591126/150‑1 (MQ=255)
                                                                                            aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGt                                                        <  2:3056698/150‑1 (MQ=255)
                                                                                               aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTa                                                     >  2:3719418/1‑150 (MQ=255)
                                                                                               aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGGGTAGGAGCGGGTTTa                                                     >  2:3719461/1‑150 (MQ=255)
                                                                                                acTGCAGGGCAACGACAAGGTGCACTTGTTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAc                                                    <  1:3777421/150‑1 (MQ=255)
                                                                                                acTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAc                                                    <  1:305213/150‑1 (MQ=255)
                                                                                                  tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTAccc                                                  >  2:1864677/1‑150 (MQ=255)
                                                                                                    cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGTGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCcgc                                                <  1:2455747/150‑1 (MQ=255)
                                                                                                              ccccggTGCACTGGCGGGCGGCTTACGTCACCAGCCTGTCGCCCCCGGCTGCGCCGGCCACGGCGGCGTACGGCCGGCGGGGCCTCTTCGCGGGTAAACCCGCTCCTACCCCGGATCCTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCg                                      <  2:3167728/146‑1 (MQ=18)
                                                                                                                        cTTGGTGGCGGCTTACGTCTACCGCCTGTCGCACCTGGCTGCGCCGGCCAAGGGTGAGTACGGCTGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGcac                            <  2:3937669/150‑1 (MQ=25)
                                                                                                                            cTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                        >  2:675550/1‑150 (MQ=255)
                                                                                                                               gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCACCTACACACGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAATGGGCCCGCAttt                     >  1:2867593/1‑150 (MQ=37)
                                                                                                                                 ggcTTACTTCTACAGCCTGTCGCACCAGGCTGCGCCGGCCACGGCTGCGTACGGCAGGCCTGGCCTCTTCGCGGTTAAACCCGCTCCTCCCCCGGATCCTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtct                   <  2:3504484/150‑1 (MQ=37)
                                                                                                                                     tACGTCTACCGCCTGTCGCACCCGGCTGAGCCGGCCACGGCTTATTAAGGCCGGCCTGGCCTCTTCGGGGTTAAACCCGCTCCTACACAGGATCCTGGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCa               <  2:1653988/150‑1 (MQ=37)
                                                                                                                                        gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGc            >  2:1370749/1‑150 (MQ=255)
                                                                                                                                          cTACCGCCTGTCGCACCAGGCGGCGCGGCCCTCGGCTGATTATGGCAGGCCGGGCCTCTACGGGGTAAACCCCGCTCCTACCCCGCATCCTGGTTGGTGCGGGTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa          <  2:2021168/150‑1 (MQ=255)
                                                                                                                                              aGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCGACACAGGATCCGTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGGACAGCATTTCTCGCAGGCAAGTGa      >  2:2045028/1‑150 (MQ=255)
                                                                                                                                                  tGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGCCCTg  >  2:2641226/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCCACCCGCAGGCGTTCATCTACGGCTCGAGCTTTGCCCAGTTGCAGCAGACCATTCGTTATGGCCGCCAGGGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTG  >  NC_002947/4844643‑4844938

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