Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 952,079 T→G 6.6% S303A (TCG→GCG)  cyoB → cytochrome bo terminal oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947952,0790TG 6.6% 602.4 / 21.4 181S303A (TCG→GCG) cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (92/77);  new base G (9/3);  total (101/80)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 186 total aligned reads displayed.

TGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGAC  >  NC_002947/951930‑952227
                                                                                                                                                     |                                                                                                                                                    
  aCTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCg                                                                                                                                                    >  2:2485391/1‑150 (MQ=255)
   cTGTTGACCGTTATCTGGACTTCCACATCTTCACCATCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                                                                                                                                                   <  1:5313316/150‑1 (MQ=255)
    tGTTGCCCGTTATCTGGACTTCCACATCTGCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                                                                                  <  1:3179408/150‑1 (MQ=255)
    tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                                                                                  >  2:3264248/1‑150 (MQ=255)
    tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                                                                                  >  1:380949/1‑150 (MQ=255)
    tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                                                                                  <  2:736376/150‑1 (MQ=255)
      ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACATCGACGt                                                                                                                                                >  1:1151190/1‑150 (MQ=255)
         aCCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtct                                                                                                                                             >  1:2212542/1‑150 (MQ=255)
           cGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTctct                                                                                                                                           <  1:5446111/150‑1 (MQ=255)
            gTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTg                                                                                                                                          >  2:3621622/1‑150 (MQ=255)
            gTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTg                                                                                                                                          <  1:4166372/150‑1 (MQ=255)
              tATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGc                                                                                                                                        <  2:4564146/150‑1 (MQ=255)
              tATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGc                                                                                                                                        <  1:1009657/150‑1 (MQ=255)
                tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                                                                                      >  2:3864144/1‑150 (MQ=255)
                tCTGGATTTCCACTTCTTCACCAATGATGGTGGAGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                                                                                      >  2:2521844/1‑150 (MQ=255)
                  tGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaagc                                                                                                                                    <  2:1109616/150‑3 (MQ=255)
                  tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                                                                                    >  2:3287074/1‑150 (MQ=255)
                      cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                                                                >  2:2440096/1‑150 (MQ=255)
                      cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                                                                <  1:4883245/150‑1 (MQ=255)
                      cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                                                                <  1:4991199/150‑1 (MQ=255)
                       ttCCACATCTGCCCCAACGAGCTTGGTGGGAACCCGAGGATGTACGTCAACCTGTTCGGGGCGCGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                                                                                               <  2:4213284/150‑1 (MQ=255)
                        tCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGt                                                                                                                              <  2:1826440/150‑1 (MQ=255)
                          cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                                                                            <  2:3087340/150‑1 (MQ=255)
                          cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                                                                            <  1:1059496/150‑1 (MQ=255)
                            caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                                                                                                                          <  2:2854200/150‑1 (MQ=255)
                             ttCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACATTTTCGGGCAAGCGTCTGTTTGGc                                                                                                                         >  2:598662/2‑150 (MQ=255)
                               cTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGcca                                                                                                                       <  2:242826/150‑1 (MQ=255)
                                ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                                                                                                                      >  2:2196022/1‑150 (MQ=255)
                                     cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                                                                                 >  2:3918321/1‑150 (MQ=255)
                                      aaCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                                                                                    <  2:2162925/146‑1 (MQ=255)
                                      aaCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                                                                                    >  1:2162925/1‑146 (MQ=255)
                                            cTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCCGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                                                                          <  2:282201/150‑1 (MQ=255)
                                                gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                                                                                      >  2:5165212/1‑150 (MQ=255)
                                                gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                                                                                      >  2:5250641/1‑150 (MQ=255)
                                                gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                                                                                      >  2:2523459/1‑150 (MQ=255)
                                                 tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                                                                                     >  2:1556058/1‑150 (MQ=255)
                                                  gggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCa                                                                                                    >  2:4893045/1‑150 (MQ=255)
                                                   ggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCAt                                                                                                   >  2:3904874/1‑150 (MQ=255)
                                                   ggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCAt                                                                                                   >  2:43499/1‑150 (MQ=255)
                                                         cGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAATAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcg                                                                                             >  1:1407906/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  2:933319/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  2:3974883/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  1:2874928/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  1:2205247/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  1:218509/1‑150 (MQ=255)
                                                          gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                            >  1:2874739/1‑150 (MQ=255)
                                                                    gTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATc                                                                                  >  2:2602451/1‑150 (MQ=255)
                                                                                 gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                                                                          <  1:3471645/113‑1 (MQ=255)
                                                                                 gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTg                                                                     >  2:2406979/1‑150 (MQ=255)
                                                                                  ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt                                                                    >  1:3064242/1‑150 (MQ=255)
                                                                                    cGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGGTTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTAt                                                                  <  2:20954/150‑1 (MQ=255)
                                                                                      tGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATgg                                                                <  1:2602451/150‑1 (MQ=255)
                                                                                       ggggTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGc                                                               <  1:3974883/150‑1 (MQ=255)
                                                                                           tCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGca                                                           >  1:754387/1‑150 (MQ=255)
                                                                                               ccTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTTCGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGcaccac                                                       <  1:2377431/150‑1 (MQ=255)
                                                                                               ccTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGcaccac                                                       <  1:2071333/150‑1 (MQ=255)
                                                                                                 tGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActt                                                     >  2:461531/1‑150 (MQ=255)
                                                                                                  gAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtc                                                    <  1:1729771/150‑1 (MQ=255)
                                                                                                  gAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtc                                                    <  1:4118558/150‑1 (MQ=255)
                                                                                                   aGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtct                                                   <  1:4641779/150‑1 (MQ=255)
                                                                                                   aGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtct                                                   <  1:1843969/150‑1 (MQ=255)
                                                                                                    ggTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtctt                                                  <  2:500153/150‑1 (MQ=255)
                                                                                                     gTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACttcttc                                                 <  2:3761555/150‑1 (MQ=255)
                                                                                                      tttACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCa                                                <  2:1668633/150‑1 (MQ=255)
                                                                                                      tttACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCa                                                >  1:2087157/1‑150 (MQ=255)
                                                                                                       ttACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCAc                                               >  1:4255369/1‑150 (MQ=255)
                                                                                                           atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGcgcg                                                                                           <  2:645478/102‑1 (MQ=255)
                                                                                                           atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATg                                           <  1:2582193/150‑1 (MQ=255)
                                                                                                           atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATg                                           >  1:3787845/1‑150 (MQ=255)
                                                                                                           atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATg                                           >  1:3564354/1‑150 (MQ=255)
                                                                                                            tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                          <  1:1256857/150‑1 (MQ=255)
                                                                                                            tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                          <  2:4642260/150‑1 (MQ=255)
                                                                                                            tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                          <  2:1011561/150‑1 (MQ=255)
                                                                                                            tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                          >  2:1516411/1‑150 (MQ=255)
                                                                                                            tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                          <  1:5098894/150‑1 (MQ=255)
                                                                                                               tgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTg                                       >  1:1880661/1‑150 (MQ=255)
                                                                                                                    ctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGgc                                  >  1:64008/1‑150 (MQ=255)
                                                                                                                     tgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGgcg                                 >  2:5201038/1‑150 (MQ=255)
                                                                                                                        cGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGcaccac                                                       <  2:5198973/125‑1 (MQ=255)
                                                                                                                        cGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGcaccac                                                       >  1:5198973/1‑125 (MQ=255)
                                                                                                                        cGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGGCCAg                             <  2:2266067/150‑1 (MQ=255)
                                                                                                                             ttCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTc                         <  2:4539813/150‑1 (MQ=255)
                                                                                                                               cGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCaa                       <  2:2424045/150‑1 (MQ=255)
                                                                                                                               cGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCaa                       <  2:1353924/150‑1 (MQ=255)
                                                                                                                                  cGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAacac                    >  1:83758/1‑150 (MQ=255)
                                                                                                                                        ctcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACcttctt              <  2:2732816/150‑1 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAAGGCCTGTTCGTCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             >  2:3505740/1‑150 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             >  1:5191852/1‑150 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             <  1:2485391/150‑1 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             >  2:5123402/1‑150 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             <  1:2523459/150‑1 (MQ=255)
                                                                                                                                         tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc             <  1:5250641/150‑1 (MQ=255)
                                                                                                                                            gAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGc          <  2:4041971/150‑1 (MQ=255)
                                                                                                                                               gTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCt                                                              >  2:2055068/1‑95 (MQ=255)
                                                                                                                                               gTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCt                                                              <  1:2055068/95‑1 (MQ=255)
                                                                                                                                                tAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTgg      >  1:1553735/1‑150 (MQ=255)
                                                                                                                                                 aaCTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGc     >  2:2021128/1‑150 (MQ=255)
                                                                                                                                                  aCTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGcg    >  1:5378436/1‑150 (MQ=255)
                                                                                                                                                    ttCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGcgac  >  1:4456504/1‑150 (MQ=255)
                                                                                                                                                    ttCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCGGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGcgac  >  2:5332993/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGAC  >  NC_002947/951930‑952227

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