Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,615,502 G→C 5.4% E160D (GAG→GAC PP_4933 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,615,5020GC 5.4% 931.2 / 30.5 275E160D (GAG→GACPP_4933hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (130/130);  new base C (3/12);  total (133/142)
Fisher's exact test for biased strand distribution p-value = 3.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 277 total aligned reads displayed.

CACACGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTCC  >  NC_002947/5615353‑5615650
                                                                                                                                                     |                                                                                                                                                     
    cGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCg                                                                                                                                                   <  1:3046877/150‑1 (MQ=255)
       ccGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTg                                                                                                                                                <  1:964207/150‑1 (MQ=255)
        cGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGc                                                                                                                                               >  1:4090212/1‑150 (MQ=255)
         gTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGcc                                                                                                                                              <  2:231583/150‑1 (MQ=255)
         gTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCGTCAACGACACCTGGCTGCCCATCCCCGACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGcc                                                                                                                                              <  2:1679925/150‑1 (MQ=255)
          tACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGACGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCCGGCACCGGACTCCGGAGcca                                                                                                                                             <  1:8528765/150‑2 (MQ=255)
          tACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCCGGCACCGGACTCCGGAGcca                                                                                                                                             <  1:6513017/150‑2 (MQ=255)
           aCAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGc                                                                                                                                           >  2:1386952/1‑150 (MQ=255)
            cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCa                                                                                                                                          >  1:2163840/1‑150 (MQ=255)
            cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCGGCACCGGACTCCGGAGccagca                                                                                                                                          <  1:6104545/150‑5 (MQ=255)
             aGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCCGGCACCGGACTCCGGAGccagca                                                                                                                                          <  2:4692604/150‑5 (MQ=255)
                aaGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGcc                                                                                                                                      <  1:4292717/150‑1 (MQ=255)
                aaGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCGGCACCGGACTCCGGAGCCAGCACgaa                                                                                                                                      <  1:6918207/150‑3 (MQ=255)
                    ggCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGt                                                                                                                                  <  2:2940426/150‑1 (MQ=255)
                    ggCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGt                                                                                                                                  <  2:2674998/150‑1 (MQ=255)
                     gCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTg                                                                                                                                 >  1:6675174/1‑150 (MQ=255)
                       cGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCa                                                                                                                               >  1:9078975/1‑150 (MQ=255)
                       cGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCa                                                                                                                               >  2:8242833/1‑150 (MQ=255)
                        gaggagTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCAt                                                                                                                              >  2:5481948/1‑150 (MQ=255)
                         aggagTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCAtt                                                                                                                             >  2:523149/1‑150 (MQ=255)
                           gagTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGc                                                                                                                           <  1:5048675/150‑1 (MQ=255)
                            agTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCg                                                                                                                          <  2:956879/150‑1 (MQ=255)
                                  tgctgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAAc                                                                                                                    >  2:5066484/1‑150 (MQ=255)
                                   gctgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACt                                                                                                                   >  2:8636344/1‑150 (MQ=255)
                                     tgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGa                                                                                                                 >  2:3392790/1‑150 (MQ=255)
                                     tgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCAGTGCGGCAAACTGa                                                                                                                 >  2:2703150/1‑150 (MQ=255)
                                       gAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATg                                                                                                               <  1:7860660/150‑1 (MQ=255)
                                        aaaCCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGc                                                                                                              <  1:5732425/150‑1 (MQ=255)
                                           ccccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTAc                                                                                                           <  2:6840558/150‑1 (MQ=255)
                                            cccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACt                                                                                                          <  1:5522416/150‑1 (MQ=255)
                                             ccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTg                                                                                                         >  2:1862723/1‑150 (MQ=255)
                                              cGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTgg                                                                                                        <  2:3690825/150‑1 (MQ=255)
                                               ggAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGa                                                                                                       <  2:7119281/150‑1 (MQ=255)
                                               ggAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGa                                                                                                       <  1:8425760/150‑1 (MQ=255)
                                                   ggTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAgcgc                                                                                                   <  1:5086484/150‑1 (MQ=255)
                                                    gTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGcc                                                                                                  <  2:6121575/150‑1 (MQ=255)
                                                         cTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACgg                                                                                             >  1:3726150/1‑150 (MQ=255)
                                                            ggcgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGcc                                                                                          >  1:443828/1‑150 (MQ=255)
                                                             gcgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCt                                                                                         <  1:2346507/150‑1 (MQ=255)
                                                              cgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACCCCCCCCTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTa                                                                                        <  1:522476/150‑1 (MQ=255)
                                                               gggcgCTTCAACGACTCCTGGCTGGCCATCCCCTACGCCGCCGCCCTTCACTACCGCGAAACTCCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACTCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAt                                                                                       <  2:3729035/150‑1 (MQ=255)
                                                                     ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCAc                                                                                 <  1:8637304/150‑1 (MQ=255)
                                                                     ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCAc                                                                                 >  1:1388333/1‑150 (MQ=255)
                                                                     ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCAc                                                                                 >  2:4449317/1‑150 (MQ=255)
                                                                     ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCAc                                                                                 >  1:1782932/1‑150 (MQ=255)
                                                                       cAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCt                                                                               >  1:3598766/1‑150 (MQ=255)
                                                                        aaCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTg                                                                              >  2:8344173/1‑150 (MQ=255)
                                                                         aCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGc                                                                             >  2:4268992/1‑150 (MQ=255)
                                                                           gACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGccc                                                                           >  1:4710377/1‑150 (MQ=255)
                                                                            acacCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCa                                                                          <  1:8802992/150‑1 (MQ=255)
                                                                                cTGGCTGCCCCTCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGc                                                                      <  1:3876904/150‑1 (MQ=255)
                                                                                cTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGc                                                                      >  1:2196658/1‑150 (MQ=255)
                                                                                     tGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGc                                                                 <  1:1460339/150‑1 (MQ=255)
                                                                                       cccATCCCCGACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTa                                                               <  2:7580721/150‑1 (MQ=255)
                                                                                        ccATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTAc                                                              >  1:2903416/1‑150 (MQ=255)
                                                                                           tCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGc                                                           >  2:5817999/1‑150 (MQ=255)
                                                                                            ccccTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCt                                                          <  2:1873983/150‑1 (MQ=255)
                                                                                             cccTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTg                                                         >  2:2351309/1‑150 (MQ=255)
                                                                                              ccTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGa                                                        >  2:7399789/1‑150 (MQ=255)
                                                                                              ccTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGa                                                        >  2:2135261/1‑150 (MQ=255)
                                                                                              ccTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCAATGCGGCAAAATGATGCTACTGGAGAGCCCACCGGCCTATGTGCACCTGCCCACCGCTACGCTACAGCTGa                                                        >  2:6410212/1‑150 (MQ=255)
                                                                                                   gccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTAt                                                   <  2:1779787/150‑1 (MQ=255)
                                                                                                   gccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTAt                                                   >  1:2744515/1‑150 (MQ=255)
                                                                                                   gccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTAt                                                   >  1:6533044/1‑150 (MQ=255)
                                                                                                   gccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTAt                                                   <  2:335399/150‑1 (MQ=255)
                                                                                                    ccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGTGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATg                                                  <  1:3375481/150‑1 (MQ=255)
                                                                                                    ccgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATg                                                  <  1:1916880/150‑1 (MQ=255)
                                                                                                     cgccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGa                                                 <  1:2951688/150‑1 (MQ=255)
                                                                                                      gccgccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGAc                                                >  1:1266818/1‑150 (MQ=255)
                                                                                                         gccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATg                                             >  2:7433053/1‑150 (MQ=255)
                                                                                                         gccCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATg                                             <  1:6871965/150‑1 (MQ=255)
                                                                                                                  tACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCaa                                    <  2:3952138/150‑1 (MQ=255)
                                                                                                                  tACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCaa                                    <  1:7150705/150‑1 (MQ=255)
                                                                                                                  tACCGCGAAACACCACACTTCGTGCTGGCTACGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCAACAGCTGATCTATGACATGCGCCTGCac                                    >  1:4946874/1‑149 (MQ=255)
                                                                                                                    ccGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGt                                  <  1:7437619/150‑1 (MQ=255)
                                                                                                                      gcgAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGAGCCTTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTCTGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGc                                >  2:8116703/1‑150 (MQ=255)
                                                                                                                        gAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGccc                              >  2:3634848/1‑150 (MQ=255)
                                                                                                                          aaCACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCg                            >  1:4220110/1‑150 (MQ=255)
                                                                                                                              ccccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGAc                        <  1:2031144/150‑1 (MQ=255)
                                                                                                                               cccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACg                       >  1:7309667/1‑150 (MQ=255)
                                                                                                                                ccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCTGGACGc                      <  2:5269917/150‑1 (MQ=255)
                                                                                                                                ccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGc                      <  1:7298485/150‑1 (MQ=255)
                                                                                                                                ccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGc                      >  1:5568083/1‑150 (MQ=255)
                                                                                                                                ccACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGc                      >  2:7956569/1‑150 (MQ=255)
                                                                                                                                   cTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCaa                   <  1:2538363/150‑1 (MQ=255)
                                                                                                                                   cTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCaa                   <  1:8202717/150‑1 (MQ=255)
                                                                                                                                     tCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCAGTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGa                 <  2:3081539/150‑1 (MQ=255)
                                                                                                                                     ggcgtgcggCACCGGACTCCGGAGCCAGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAg                  >  1:8598914/8‑150 (MQ=255)
                                                                                                                                       gTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGaaa               <  2:4501237/150‑1 (MQ=255)
                                                                                                                                       gTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGaaa               <  1:3641528/150‑1 (MQ=255)
                                                                                                                                        tGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAAc              >  1:1121316/1‑150 (MQ=255)
                                                                                                                                        tGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAAc              >  1:6571458/1‑150 (MQ=255)
                                                                                                                                        tGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAAc              <  2:4190771/150‑1 (MQ=255)
                                                                                                                                           tGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCa           >  1:586770/1‑150 (MQ=255)
                                                                                                                                           tGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCa           <  1:8304137/150‑1 (MQ=255)
                                                                                                                                           tGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCa           >  2:1025251/1‑150 (MQ=255)
                                                                                                                                              cTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCt        <  2:903441/150‑1 (MQ=255)
                                                                                                                                               tCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgc                                            >  1:4467367/1‑113 (MQ=255)
                                                                                                                                               tCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCtt       >  1:309719/1‑150 (MQ=255)
                                                                                                                                                    aGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTcc  <  2:6523803/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACACGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCCGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCTTGTTCC  >  NC_002947/5615353‑5615650

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: