Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 952,033 T→C 7.1% V287V (GTT→GTC 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,0330TC 7.1% 966.8 / 38.2 310V287V (GTT→GTCcyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (163/125);  new base C (12/10);  total (175/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

TGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCAC  >  NC_002947/951885‑952182
                                                                                                                                                    |                                                                                                                                                     
  aTCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTat                                                                                                                                                    <  1:5983972/150‑2 (MQ=255)
      tcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcc                                                                                                                                                <  2:2671660/150‑1 (MQ=255)
        gCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctg                                                                                                                                              >  1:8649762/1‑150 (MQ=255)
          ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgat                                                                                                                                            >  1:2371795/1‑150 (MQ=255)
          ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACAGTTATCTGGACTTCAACATCTTCACCAACGAGCTTGGTGGGAACACTATGATGTAAGTCAACCTGTTCTGGGCGTGGGGTCACACTGAGCTTTACatcctgat                                                                                                                                            >  2:1665238/1‑150 (MQ=255)
                 ccGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGcct                                                                                                                                     >  2:8730347/1‑149 (MQ=255)
                  cGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCgg                                                                                                                                    >  2:2873146/1‑150 (MQ=255)
                  cGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCgg                                                                                                                                    >  1:8314218/1‑150 (MQ=255)
                       cTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTc                                                                                                                               >  2:4475312/1‑150 (MQ=255)
                        tGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCg                                                                                                                              <  1:315532/150‑1 (MQ=255)
                          aCTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                                                                                            <  1:2183406/150‑1 (MQ=255)
                          aCTGCTGCACTCGCTCTGCTGACTGTTGACAGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCAGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCATGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                                                                                            <  2:5733309/150‑1 (MQ=255)
                            tgctgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt                                                                                                                          >  1:4364174/1‑150 (MQ=255)
                              tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGgtgttt                                                                                                                        >  1:6955801/6‑148 (MQ=255)
                              tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGgtgttt                                                                                                                        >  1:8506216/6‑148 (MQ=255)
                               tgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGtt                                                                                                                       <  1:8124583/150‑1 (MQ=255)
                               tgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGtt                                                                                                                       <  1:3466373/150‑1 (MQ=255)
                                 cACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtct                                                                                                                     >  1:7151703/1‑150 (MQ=255)
                                  aCTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtctc                                                                                                                    <  1:3479222/150‑1 (MQ=255)
                                  aCTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtctc                                                                                                                    <  2:4287660/150‑1 (MQ=255)
                                  aCTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtctc                                                                                                                    <  2:369507/150‑1 (MQ=255)
                                   cTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCg                                                                                                                   <  2:7045349/150‑1 (MQ=255)
                                     cGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGa                                                                                                                 <  2:7441693/150‑1 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCATCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                <  2:1521936/150‑1 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                <  2:6913582/150‑1 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                >  1:2795524/1‑150 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                <  1:74868/150‑1 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                <  2:918724/150‑1 (MQ=255)
                                      gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                                <  1:306845/150‑1 (MQ=255)
                                       ctctGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAg                                                                                                               <  1:4585744/150‑1 (MQ=255)
                                          tgctgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTaa                                                                                                            <  1:3818615/150‑1 (MQ=255)
                                          tgctgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTaa                                                                                                            >  2:3850796/1‑150 (MQ=255)
                                              gACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTc                                                                                                        <  1:4904620/150‑1 (MQ=255)
                                              gACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTc                                                                                                        >  2:625610/1‑150 (MQ=255)
                                               aCTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCg                                                                                                       >  1:1748041/1‑150 (MQ=255)
                                                cTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                                                                                                      <  1:8583371/150‑1 (MQ=255)
                                                    tGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGtt                                                                                                  >  2:3736330/1‑150 (MQ=255)
                                                         gTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTg                                                                                             >  2:8451896/1‑150 (MQ=255)
                                                             tCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCaa                                                                                         >  1:5756437/1‑150 (MQ=255)
                                                              cTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCaaa                                                                                        >  1:7268620/1‑150 (MQ=255)
                                                               tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                                       >  1:8052415/1‑150 (MQ=255)
                                                               tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                                       >  1:4097250/1‑150 (MQ=255)
                                                               tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                                       >  2:4120387/1‑150 (MQ=255)
                                                                ggACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACg                                                                                      >  2:6737069/1‑150 (MQ=255)
                                                                   cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                   <  1:6086960/150‑1 (MQ=255)
                                                                   cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                                   <  1:2462774/150‑1 (MQ=255)
                                                                    ttCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                                                  >  1:5226281/1‑150 (MQ=255)
                                                                      ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt                                                                                >  1:6491334/1‑150 (MQ=255)
                                                                      ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt                                                                                >  1:5823767/1‑150 (MQ=255)
                                                                      ccACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGtt                                                                                >  1:1162908/1‑150 (MQ=255)
                                                                       cactTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGttt                                                                               <  2:2361374/150‑2 (MQ=255)
                                                                       cactTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGttt                                                                               <  1:3421972/150‑2 (MQ=255)
                                                                       cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                               >  1:1221959/1‑150 (MQ=255)
                                                                       cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                               <  1:3759914/150‑1 (MQ=255)
                                                                        acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                                                                              >  1:662962/1‑150 (MQ=255)
                                                                          aTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGc                                                                            <  2:7287500/150‑1 (MQ=255)
                                                                           tCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGcc                                                                           >  1:5471606/1‑150 (MQ=255)
                                                                             ttCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGTTTGGctat                                                                         >  1:944367/1‑147 (MQ=255)
                                                                             ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                                                                         >  1:1295883/1‑150 (MQ=255)
                                                                              tCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGTTTGGctatg                                                                        >  2:4791971/1‑146 (MQ=255)
                                                                               cACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                                       >  2:8977048/1‑150 (MQ=255)
                                                                               cACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                                       <  1:5061096/150‑1 (MQ=255)
                                                                               cACCAACGAGCTTGGTGGGAACCCGATGATGAACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCGGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGAAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                                       >  1:2038597/1‑150 (MQ=255)
                                                                                aCCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCccccac                                                                      >  2:1804768/1‑150 (MQ=255)
                                                                                aCCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCcaccac                                                                      >  2:8052975/1‑150 (MQ=255)
                                                                                 ccAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACt                                                                     >  1:7955737/1‑150 (MQ=255)
                                                                                  cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                                    >  1:6589700/1‑150 (MQ=255)
                                                                                  cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                                    <  1:3004310/150‑1 (MQ=255)
                                                                                     cGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCgat                                                                 <  2:3525837/150‑1 (MQ=255)
                                                                                     cGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCgat                                                                 >  2:861454/1‑150 (MQ=255)
                                                                                     cGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCgat                                                                 <  2:6073058/150‑1 (MQ=255)
                                                                                         cTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                             <  2:5808570/150‑1 (MQ=255)
                                                                                              tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                                        >  2:7598939/1‑150 (MQ=255)
                                                                                              tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                                        >  2:8932995/1‑150 (MQ=255)
                                                                                              tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                                        >  1:7620646/1‑150 (MQ=255)
                                                                                                           atgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATc                                           <  1:3971492/150‑1 (MQ=255)
                                                                                                           atgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATc                                           <  1:8927074/150‑1 (MQ=255)
                                                                                                              tACGTCATCCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGcc                                        <  1:2581616/150‑1 (MQ=255)
                                                                                                                cGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCAt                                      <  2:8271566/150‑1 (MQ=255)
                                                                                                                  tCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAATCGCCTGTTCGGCGACCACTCGTTGATTTACGCATCGGGCgcgatcgcgatcg                                    >  2:5629713/1‑149 (MQ=255)
                                                                                                                   cAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCt                                   >  1:370763/1‑150 (MQ=255)
                                                                                                                     aCCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCg                                 >  1:5534322/1‑150 (MQ=255)
                                                                                                                      ccTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCgg                                <  2:7268620/150‑1 (MQ=255)
                                                                                                                       cTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGc                               >  1:3603606/1‑150 (MQ=255)
                                                                                                                       cTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGc                               >  1:1835738/1‑150 (MQ=255)
                                                                                                                            cTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGc                          <  1:5236683/150‑1 (MQ=255)
                                                                                                                             tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt                         >  2:5032995/1‑150 (MQ=255)
                                                                                                                             tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt                         >  2:6586583/1‑150 (MQ=255)
                                                                                                                             tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt                         >  2:4748431/1‑150 (MQ=255)
                                                                                                                               ggCGTGGGGTCACCCTGAGGTTTACATCCTGTTCCTGCGGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAGTTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGTCTTTTCTGt                       >  2:4881497/1‑150 (MQ=255)
                                                                                                                               ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt                       >  1:2241700/1‑150 (MQ=255)
                                                                                                                                   tGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATgg                   <  1:1665238/150‑1 (MQ=255)
                                                                                                                                    ggggTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGc                  <  2:5276957/150‑1 (MQ=255)
                                                                                                                                    ggggTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGc                  <  2:2128211/150‑1 (MQ=255)
                                                                                                                                      ggTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTg                >  2:1851155/1‑150 (MQ=255)
                                                                                                                                        tCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGca              >  1:5044504/1‑150 (MQ=255)
                                                                                                                                            ccTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGcaccac          >  1:7065537/1‑150 (MQ=255)
                                                                                                                                             cTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCAct         >  1:7273345/1‑150 (MQ=255)
                                                                                                                                                  gTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGACAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACttcttc    <  1:281776/150‑1 (MQ=255)
                                                                                                                                                   tttACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCa   >  1:1506464/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCACCACTTCTTCAC  >  NC_002947/951885‑952182

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: