Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 914,634 G→C 7.8% A448P (GCC→CCC) ‡ fruA → PTS fructose transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947914,6340GC 7.8% 620.7 / 32.9 192A448P (GCC→CCC) ‡fruAPTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base G (96/80);  new base C (5/10);  total (101/91)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 200 total aligned reads displayed.

GCATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAG  >  NC_002947/914485‑914783
                                                                                                                                                      |                                                                                                                                                       
    gggCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGgccgc                                                                                                                                                     >  2:548013/1‑150 (MQ=255)
      gCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGccgccg                                                                                                                                                   >  2:4988984/1‑150 (MQ=255)
          caccaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CAt                                                                                                                                              >  1:2420441/1‑150 (MQ=255)
           acctATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTGCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGGCCATGGTCGCCGC‑CAtt                                                                                                                                             <  2:4161072/150‑2 (MQ=255)
           accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑Cagg                                                                                                                                             >  2:3199378/1‑148 (MQ=255)
           accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATg                                                                                                                                             <  2:1851204/150‑1 (MQ=255)
           accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATg                                                                                                                                             <  2:2928390/150‑1 (MQ=255)
              aaTGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGAGCCAATCGGCGGC‑ACCATGCCGGCGGC‑CATGGtc                                                                                                                                          <  1:4321250/150‑2 (MQ=255)
              aaTGCCATCCTGCTGGGACTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGCCGGCGGC‑CATGGtc                                                                                                                                          <  1:5083340/150‑2 (MQ=255)
                       ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGC‑ACCATGCCGGCGGC‑CATGGTCgccgccatt                                                                                                                                 <  1:1304262/150‑4 (MQ=37)
                       ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGAtt                                                                                                                                <  1:1671274/150‑1 (MQ=255)
                       ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGAtt                                                                                                                                <  1:3936690/150‑1 (MQ=255)
                        tgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCAGGGCCGCCGG‑CCTGGT‑GCCGCCGATTg                                                                                                                               >  2:686964/1‑150 (MQ=255)
                        tgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTAGGCGGGCCGATCAACAAGGCCGCCTATGCGTACTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTg                                                                                                                               >  1:2112148/1‑150 (MQ=255)
                         gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGC‑ACCATGCCGGCGGC‑CATGGTCGCCGCC‑ATTgg                                                                                                                              <  1:1266523/150‑1 (MQ=37)
                         gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTgg                                                                                                                              <  1:2271183/150‑1 (MQ=255)
                          ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGc                                                                                                                             <  2:1735840/150‑1 (MQ=255)
                            ggggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                           >  1:2222559/1‑150 (MQ=255)
                             gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                          <  2:3405989/150‑1 (MQ=255)
                              ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTgg                                                                                                                         <  1:98649/150‑1 (MQ=255)
                              ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTgg                                                                                                                         <  1:3668300/150‑1 (MQ=255)
                              ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTgg                                                                                                                         <  2:4757905/150‑1 (MQ=255)
                                   ttGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCACT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTCtt                                                                                                                    >  2:3736248/1‑150 (MQ=255)
                                    tGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTg                                                                                                                   <  2:2391863/150‑1 (MQ=255)
                                       tGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCa                                                                                                                >  1:5083092/1‑150 (MQ=255)
                                         ggcggcATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGGTTGGGTATTGCCAcc                                                                                                              >  1:3089150/1‑150 (MQ=255)
                                              cttgatgTGCGTCGAGCTTGGCGGGCCGATCATCAAGTCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCt                                                                                                         <  2:4895385/148‑1 (MQ=255)
                                              cATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCt                                                                                                         >  1:2642951/1‑150 (MQ=255)
                                                      gCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CCTGGT‑GCCGCCGAGTGGCGTGGGTCTTGCCACCTTCCTGGCCCGgc                                                                                                 >  2:2364996/1‑150 (MQ=255)
                                                      gCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGgc                                                                                                 <  1:1934923/150‑1 (MQ=255)
                                                        gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                          <  1:2719464/123‑1 (MQ=255)
                                                        gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                          >  2:2719464/1‑123 (MQ=255)
                                                          cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                             <  2:1018441/150‑1 (MQ=255)
                                                          cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                             <  1:59715/150‑1 (MQ=255)
                                                             ccTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGAATTGCCACCTTCCTGGCCCGGCGCAAGtt                                                                                          >  1:664661/1‑150 (MQ=255)
                                                              cTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGAGTGGCTGGGGTATTGCCACCTTCCTGGCCCGGCGCACGTTc                                                                                         >  2:859534/1‑150 (MQ=255)
                                                                 ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCCTGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                      <  1:2890887/150‑1 (MQ=255)
                                                                 ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                      <  2:1652125/150‑1 (MQ=255)
                                                                 ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                      <  1:4262032/150‑1 (MQ=255)
                                                                 ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                      <  1:3134794/150‑1 (MQ=255)
                                                                     ggCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCaga                                                                                  <  2:3619692/150‑1 (MQ=255)
                                                                       ccGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                                <  1:4459568/150‑1 (MQ=255)
                                                                       ccGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                                >  2:3349321/1‑150 (MQ=255)
                                                                         gATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcga                                                                              >  1:2347740/1‑150 (MQ=255)
                                                                         gATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCGGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcga                                                                              >  1:155035/1‑150 (MQ=255)
                                                                           tCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgagc                                                                            <  2:1264025/150‑1 (MQ=255)
                                                                            caacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGagcgagcg                                                                           >  1:1097321/1‑150 (MQ=255)
                                                                              acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                         >  1:3838861/1‑150 (MQ=255)
                                                                              acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                         >  1:123823/1‑150 (MQ=255)
                                                                              acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                         >  2:2488736/1‑150 (MQ=255)
                                                                              acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                         >  1:4434351/1‑150 (MQ=255)
                                                                               caaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGa                                                                        <  2:1310094/150‑1 (MQ=255)
                                                                                aaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAg                                                                       <  2:2194550/150‑1 (MQ=255)
                                                                                   gccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCg                                                                    <  2:3033578/150‑1 (MQ=255)
                                                                                   gccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCg                                                                    <  2:5083572/150‑1 (MQ=255)
                                                                                    ccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                   <  1:4590168/150‑1 (MQ=255)
                                                                                      gccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                 <  1:4559214/150‑1 (MQ=255)
                                                                                           tGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CAGGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAggc                                                            >  2:3945930/1‑150 (MQ=255)
                                                                                                      tGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctggcgc                                                 <  1:3349321/150‑1 (MQ=255)
                                                                                                         ggCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                              <  2:2995243/150‑1 (MQ=255)
                                                                                                          gctgctGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                             >  2:4460740/1‑150 (MQ=255)
                                                                                                          gctgctGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                             <  1:1081381/150‑1 (MQ=255)
                                                                                                          gctgctGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                             >  1:73257/1‑150 (MQ=255)
                                                                                                           ctgctgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTgggg                                            <  2:4709720/150‑1 (MQ=255)
                                                                                                              ctgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtg                                         >  1:86965/1‑150 (MQ=255)
                                                                                                               tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CTTGGT‑GCCTCCGATTGGCTTGGGTATTGCGACCTTCCTGGCCCGGCGCAAGTCCGACCTTATCGAGCGCTAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                        <  1:1678114/150‑1 (MQ=255)
                                                                                                               tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                        >  2:915199/1‑150 (MQ=255)
                                                                                                               tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                        >  2:915168/1‑150 (MQ=255)
                                                                                                               tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                        <  1:4902957/150‑1 (MQ=255)
                                                                                                               tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCAAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                        <  2:3036223/150‑1 (MQ=255)
                                                                                                                   cTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                    >  2:4381995/1‑150 (MQ=255)
                                                                                                                   cTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                    >  2:2893503/1‑150 (MQ=255)
                                                                                                                    tcgtcgAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                   >  1:3590767/1‑150 (MQ=255)
                                                                                                                       tcgAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtc                                >  1:1755685/1‑150 (MQ=255)
                                                                                                                       tcgAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtc                                >  1:3031568/1‑150 (MQ=255)
                                                                                                                         gAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                              >  2:449503/1‑150 (MQ=255)
                                                                                                                          aGCTACGCGCCAAT‑GGCGGCGACCATGGCCTGCGG‑AATGGT‑GCCGCCGATTGTCTTGGGTATTGCCTCCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTcc                             <  1:3807989/150‑1 (MQ=255)
                                                                                                                          aGCTACGCGCCAAT‑GGCGGCGACCATGACCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTcc                             <  1:5100241/150‑1 (MQ=255)
                                                                                                                             tACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAg                          <  2:4250192/150‑1 (MQ=255)
                                                                                                                              aCGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgg                         <  2:2112148/150‑1 (MQ=255)
                                                                                                                                  gcCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                              >  2:4837417/1‑141 (MQ=255)
                                                                                                                                  gcCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                              <  1:4837417/141‑1 (MQ=255)
                                                                                                                                  gcCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGc                     <  2:2861207/150‑1 (MQ=255)
                                                                                                                                      aT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAtt                 <  1:5056480/150‑1 (MQ=255)
                                                                                                                                      aT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAtt                 <  2:573807/150‑1 (MQ=255)
                                                                                                                                      aT‑GCCGGCGGCCCGGGTCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAtt                 >  2:2629153/1‑150 (MQ=255)
                                                                                                                                       t‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTc                >  1:4349204/1‑150 (MQ=255)
                                                                                                                                          gcggcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCg              <  2:674630/150‑1 (MQ=255)
                                                                                                                                           gcggcaCCATGCCGGCGGC‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGt             >  1:3116012/6‑150 (MQ=255)
                                                                                                                                             gcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt           >  1:201527/1‑150 (MQ=255)
                                                                                                                                             gcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt           <  1:635547/150‑1 (MQ=255)
                                                                                                                                             gcgACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt           <  1:1456652/150‑1 (MQ=255)
                                                                                                                                               gACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgc         >  2:1292191/1‑150 (MQ=255)
                                                                                                                                                  cATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc      >  1:3739709/1‑150 (MQ=255)
                                                                                                                                                  cATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc      >  2:2802417/1‑150 (MQ=255)
                                                                                                                                                    tGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCa    >  2:1783326/1‑150 (MQ=255)
                                                                                                                                                     ggCCGCCGG‑CATTGT‑GCCGCCGATTGTCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa   <  2:2072482/150‑1 (MQ=255)
                                                                                                                                                      gccggcGGCCATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa   <  2:3116012/150‑1 (MQ=255)
                                                                                                                                                      gccgccGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAg  <  1:4050930/150‑1 (MQ=255)
                                                                                                                                                      gccgccGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAg  <  1:557840/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GCATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAG  >  NC_002947/914485‑914783

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: ‑