Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 952,084 G→C 6.5% T304T (ACG→ACC 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,0840GC 6.5% 591.3 / 21.0 180T304T (ACG→ACCcyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (92/76);  new base C (8/3);  total (101/79)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 182 total aligned reads displayed.

TTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATGCT  >  NC_002947/951936‑952233
                                                                                                                                                    |                                                                                                                                                     
ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACATCGACGt                                                                                                                                                      >  1:1151190/1‑150 (MQ=255)
   aCCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtct                                                                                                                                                   >  1:2212542/1‑150 (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:2433166/1‑150 (MQ=255)
          tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                                                                                            >  2:2522080/1‑150 (MQ=255)
            tGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaagc                                                                                                                                          <  2:1109616/150‑3 (MQ=255)
            tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                                                                                          >  2:3287074/1‑150 (MQ=255)
              gACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGc                                                                                                                                        <  1:3621622/150‑1 (MQ=255)
                cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                                                                      <  1:4991199/150‑1 (MQ=255)
                 ttCCACATCTGCCCCAACGAGCTTGGTGGGAACCCGAGGATGTACGTCAACCTGTTCGGGGCGCGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                                                                                                     <  2:4213284/150‑1 (MQ=255)
                  tCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGTCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGt                                                                                                                                    >  1:5197280/1‑150 (MQ=255)
                  tCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGt                                                                                                                                    <  2:1826440/150‑1 (MQ=255)
                    cactTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGttt                                                                                                                                  <  1:4674435/150‑2 (MQ=255)
                    cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                                                                                  <  1:1059496/150‑1 (MQ=255)
                    cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                                                                                  <  2:3087340/150‑1 (MQ=255)
                     acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                                                                                                                                 >  2:634892/1‑150 (MQ=255)
                      caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                                                                                                                                <  2:2854187/150‑1 (MQ=255)
                          ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                                                                                                                            >  2:2196022/1‑150 (MQ=255)
                               cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                                                                                       >  2:3918321/1‑150 (MQ=255)
                                     gCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCgatgat                                                                                                                 <  1:2381920/150‑1 (MQ=255)
                                          gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                                                                                            >  2:2523459/1‑150 (MQ=255)
                                          gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                                                                                            >  2:5250641/1‑150 (MQ=255)
                                           tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                                                                                           >  2:1556058/1‑150 (MQ=255)
                                            gggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCa                                                                                                          >  2:4893045/1‑150 (MQ=255)
                                             ggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCAt                                                                                                         >  2:43499/1‑150 (MQ=255)
                                               aaCCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCg                                                                                                       <  1:2196022/150‑1 (MQ=255)
                                                aCCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCgg                                                                                                      <  1:217120/150‑1 (MQ=255)
                                                 cccGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCGAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCggg                                                                                                     >  1:3793889/1‑150 (MQ=255)
                                                  ccGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgc                                                                                                    >  2:3892218/1‑150 (MQ=255)
                                                  ccGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGGTCTCTGGCAAACGCCTGGTCGGCCACCACGCGATGAGTTACGCATCGGgc                                                                                                    >  2:4626786/1‑150 (MQ=255)
                                                   cGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAATAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcg                                                                                                   >  1:1407906/1‑150 (MQ=255)
                                                    gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                                  >  2:3974883/1‑150 (MQ=255)
                                                    gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                                  >  2:933319/1‑150 (MQ=255)
                                                    gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                                  >  1:2874739/1‑150 (MQ=255)
                                                    gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                                                                                  >  1:2874928/1‑150 (MQ=255)
                                                      tgatgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGa                                                                                                <  2:5166102/150‑1 (MQ=255)
                                                              gTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATc                                                                                        >  2:4888466/1‑150 (MQ=255)
                                                                          tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt                                                                            >  2:2424429/1‑150 (MQ=255)
                                                                          tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt                                                                            <  1:2101582/150‑1 (MQ=255)
                                                                           gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                                                                                <  1:3471645/113‑1 (MQ=255)
                                                                           gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTg                                                                           >  2:2406979/1‑150 (MQ=255)
                                                                           gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTGACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTg                                                                           >  2:2406986/1‑150 (MQ=255)
                                                                            ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt                                                                          >  1:4504641/1‑150 (MQ=255)
                                                                            ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt                                                                          >  1:3064242/1‑150 (MQ=255)
                                                                                tGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATgg                                                                      <  1:2602451/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)
                                                                                            gAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTTGCAAACGCCTGTTCGTCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtc                                                          <  1:2527954/150‑1 (MQ=255)
                                                                                            gAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCGAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtc                                                          <  2:3793889/150‑1 (MQ=255)
                                                                                            gAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtc                                                          <  1:4118558/150‑1 (MQ=255)
                                                                                             aGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtct                                                         <  1:1843666/150‑1 (MQ=255)
                                                                                              ggTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtctt                                                        <  2:5310111/150‑1 (MQ=255)
                                                                                              ggTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtctt                                                        <  2:500153/150‑1 (MQ=255)
                                                                                              ggTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCActtctt                                                        <  2:70520/150‑1 (MQ=255)
                                                                                               gTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACttcttc                                                       <  2:3761555/150‑1 (MQ=255)
                                                                                                tttACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCa                                                      <  2:2474289/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)
                                                                                                  tACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCAcc                                                    <  2:2117276/150‑1 (MQ=255)
                                                                                                     atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATg                                                 <  1:2582193/150‑1 (MQ=255)
                                                                                                     atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATg                                                 >  1:3564354/1‑150 (MQ=255)
                                                                                                     atcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTGCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGcgcg                                                                                                 >  1:645478/1‑102 (MQ=255)
                                                                                                      tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                <  1:1256857/150‑1 (MQ=255)
                                                                                                      tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                <  1:5098894/150‑1 (MQ=255)
                                                                                                      tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                <  2:1011561/150‑1 (MQ=255)
                                                                                                      tcctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                >  2:1516411/1‑150 (MQ=255)
                                                                                                        ctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTTGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACGGCTTCACCATGGGt                                              >  2:2729000/1‑150 (MQ=255)
                                                                                                        ctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                >  2:555601/1‑148 (MQ=255)
                                                                                                        ctgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATgg                                                <  1:555601/148‑1 (MQ=255)
                                                                                                         tgatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTg                                             >  1:1880661/1‑150 (MQ=255)
                                                                                                          gatcctgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGc                                            >  2:3058614/1‑150 (MQ=255)
                                                                                                               tgCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGgcg                                       >  2:5201038/1‑150 (MQ=255)
                                                                                                                 ccGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGcc                                     >  1:2519530/1‑150 (MQ=255)
                                                                                                                  cGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGGCCAg                                   <  2:2266067/150‑1 (MQ=255)
                                                                                                                   ggCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGc                                            >  2:2060140/1‑141 (MQ=255)
                                                                                                                       ttCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTc                               <  2:4539813/150‑1 (MQ=255)
                                                                                                                         cGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCaa                             <  2:1353924/150‑1 (MQ=255)
                                                                                                                            cGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAacac                          >  1:83758/1‑150 (MQ=255)
                                                                                                                                ttCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACcttc                      <  2:5157657/150‑1 (MQ=255)
                                                                                                                                   tcGGAAGTAACTTCGACGTTCTCTGGCAAAGGCCTGTTCGTCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc                   >  2:3505740/1‑150 (MQ=255)
                                                                                                                                   tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc                   <  1:5250641/150‑1 (MQ=255)
                                                                                                                                   tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc                   >  1:5191852/1‑150 (MQ=255)
                                                                                                                                   tcGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCttcttc                   >  2:2765293/1‑150 (MQ=255)
                                                                                                                                      gAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGc                <  2:4041971/150‑1 (MQ=255)
                                                                                                                                         gTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCt                                                                    <  1:2055068/95‑1 (MQ=255)
                                                                                                                                         gTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCt                                                                    >  2:2055068/1‑95 (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:1624694/1‑150 (MQ=255)
                                                                                                                                              ttCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGcgac        >  1:4456504/1‑150 (MQ=255)
                                                                                                                                                  aCGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATg    <  1:3264248/150‑1 (MQ=255)
                                                                                                                                                  aCGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATg    <  1:5201038/150‑1 (MQ=255)
                                                                                                                                                  aCGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATg    <  2:2590354/150‑1 (MQ=255)
                                                                                                                                                    gTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATGCt  <  2:2087157/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCACCATGGGTGCCGGCGCCAGCGTCAACACCTTCTTCGGCCTGGCGACGATGCT  >  NC_002947/951936‑952233

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: