Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 914,636 C→G 7.9% A448A (GCC→GCG) ‡ fruA → PTS fructose transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947914,6360CG 7.9% 621.0 / 32.3 191A448A (GCC→GCG) ‡fruAPTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base C (94/82);  new base G (5/10);  total (99/92)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 201 total aligned reads displayed.

ATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGA  >  NC_002947/914487‑914785
                                                                                                                                                        |                                                                                                                                                       
        caccaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CAt                                                                                                                                                >  1:2420441/1‑150 (MQ=255)
         accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑Cagg                                                                                                                                               >  2:3199378/1‑148 (MQ=255)
         accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATg                                                                                                                                               <  2:1851204/150‑1 (MQ=255)
         accaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATg                                                                                                                                               <  2:2928390/150‑1 (MQ=255)
            aaTGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGAGCCAATCGGCGGC‑‑‑ACCATGCCGGCGGC‑CATGGtc                                                                                                                                            <  1:4321250/150‑2 (MQ=255)
            aaTGCCATCCTGCTGGGACTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGCCGGCGGC‑CATGGtc                                                                                                                                            <  1:5083382/150‑2 (MQ=255)
            aaTGCCATCCTGCTGGGACTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGCCGGCGGC‑CATGGtc                                                                                                                                            <  1:5083338/150‑2 (MQ=255)
                   tCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGa                                                                                                                                    <  1:5006938/150‑1 (MQ=255)
                     ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGC‑‑‑ACCATGCCGGCGGC‑CATGGTCgccgccatt                                                                                                                                   <  1:1304262/150‑4 (MQ=37)
                     ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGAtt                                                                                                                                  <  1:1671274/150‑1 (MQ=255)
                     ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGAtt                                                                                                                                  <  1:3936690/150‑1 (MQ=255)
                      tgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTAGGCGGGCCGATCAACAAGGCCGCCTATGCGTACTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTg                                                                                                                                 >  1:2112148/1‑150 (MQ=255)
                       gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGC‑‑‑ACCATGCCGGCGGC‑CATGGTCGCCGCC‑ATTgg                                                                                                                                <  1:1266523/150‑1 (MQ=37)
                       gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCGGCACCATGCCGGCGGC‑CATGGTCgccgccatt                                                                                                                                   <  2:678960/150‑4 (MQ=38)
                       gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTgg                                                                                                                                <  1:2271183/150‑1 (MQ=255)
                          ggggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                             >  1:2222559/1‑150 (MQ=255)
                          ggggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCtt                                                                                                                             <  2:471254/150‑1 (MQ=255)
                           gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                            <  2:3405989/150‑1 (MQ=255)
                            ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTgg                                                                                                                           <  1:3668300/150‑1 (MQ=255)
                            ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTgg                                                                                                                           <  1:98649/150‑1 (MQ=255)
                               tGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTa                                                                                                                        >  1:3866729/1‑150 (MQ=255)
                                  tGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTg                                                                                                                     <  2:2391863/150‑1 (MQ=255)
                                     tGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCgg                                                                                                                                                    <  2:4305511/118‑1 (MQ=255)
                                       ggcggcATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGGTTGGGTATTGCCAcc                                                                                                                >  1:3089150/1‑150 (MQ=255)
                                            cATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCt                                                                                                           <  2:4895319/150‑1 (MQ=255)
                                            cATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCt                                                                                                           >  1:2642951/1‑150 (MQ=255)
                                               gatgTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CAGGGT‑GCCGCCGAGGGGCTGGGGTCTGGCCACCTTCCTGGc                                                                                                        >  2:4137975/1‑150 (MQ=255)
                                                    gCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CCTGGT‑GCCGCCGAGTGGCGTGGGTCTTGCCACCTTCCTGGCCCGgc                                                                                                   >  2:2364996/1‑150 (MQ=255)
                                                    gCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGgc                                                                                                   <  1:1934923/150‑1 (MQ=255)
                                                      gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                            >  2:2719464/1‑123 (MQ=255)
                                                      gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTg                                                                                                                            <  1:2719464/123‑1 (MQ=255)
                                                        cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                               <  2:1018441/150‑1 (MQ=255)
                                                        cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                               <  1:59715/150‑1 (MQ=255)
                                                           ccTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGAATTGCCACCTTCCTGGCCCGGCGCAAGtt                                                                                            >  1:664661/1‑150 (MQ=255)
                                                             ttGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCg                                                                                          >  2:5056480/1‑150 (MQ=255)
                                                               ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCCTGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                        <  1:2890887/150‑1 (MQ=255)
                                                               ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                        <  1:4262032/150‑1 (MQ=255)
                                                               ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                        <  2:1652125/150‑1 (MQ=255)
                                                               ggCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGcc                                                                                        <  1:3134794/150‑1 (MQ=255)
                                                                     ccGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                                  <  1:4459568/150‑1 (MQ=255)
                                                                     ccGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                                  <  2:1774723/150‑1 (MQ=255)
                                                                     ccGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                                  >  2:3349321/1‑150 (MQ=255)
                                                                       gATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcga                                                                                >  1:2347740/1‑150 (MQ=255)
                                                                       gATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCGGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcga                                                                                >  1:155035/1‑150 (MQ=255)
                                                                        aTCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                               <  1:548013/150‑1 (MQ=255)
                                                                         tCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgagc                                                                              <  2:1264025/150‑1 (MQ=255)
                                                                          caacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGagcgagcg                                                                             >  1:1097321/1‑150 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                           >  2:2488736/1‑150 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                           >  1:123823/1‑150 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                           >  1:4434351/1‑150 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                           >  2:665299/1‑150 (MQ=255)
                                                                             caaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGa                                                                          <  2:1310094/150‑1 (MQ=255)
                                                                              aaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAg                                                                         <  2:2194550/150‑1 (MQ=255)
                                                                                 gccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCg                                                                      <  2:3033578/150‑1 (MQ=255)
                                                                                  ccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                     <  1:4590168/150‑1 (MQ=255)
                                                                                    gccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                   <  1:4559214/150‑1 (MQ=255)
                                                                                            gTTCTCGGTGGGGCTGCTGGCCTCGTAGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGAGTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCGCggcggc                                                           >  2:3658645/1‑150 (MQ=255)
                                                                                              tctcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTGGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgc                                                         >  2:5091796/1‑150 (MQ=255)
                                                                                                tcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctg                                                       <  2:1788/150‑1 (MQ=255)
                                                                                                tcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctg                                                       >  2:848340/1‑150 (MQ=255)
                                                                                                    tGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctggcgc                                                   <  1:3349321/150‑1 (MQ=255)
                                                                                                       ggCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                                <  2:2995243/150‑1 (MQ=255)
                                                                                                       ggCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                                <  2:439272/150‑1 (MQ=255)
                                                                                                        gctgctGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                               <  1:1081381/150‑1 (MQ=255)
                                                                                                        gctgctGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                               >  1:73257/1‑150 (MQ=255)
                                                                                                            ctgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtg                                           >  1:86965/1‑150 (MQ=255)
                                                                                                             tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CTTGGT‑GCCTCCGATTGGCTTGGGTATTGCGACCTTCCTGGCCCGGCGCAAGTCCGACCTTATCGAGCGCTAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                          <  1:1678114/150‑1 (MQ=255)
                                                                                                             tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTTGATTGGCTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                          >  2:915188/1‑150 (MQ=255)
                                                                                                             tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                          >  2:915052/1‑150 (MQ=255)
                                                                                                             tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                          <  1:4902957/150‑1 (MQ=255)
                                                                                                             tgGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCAAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                          <  2:3036223/150‑1 (MQ=255)
                                                                                                               gCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                        >  2:3932818/1‑150 (MQ=255)
                                                                                                                 cTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                      >  2:550029/1‑150 (MQ=255)
                                                                                                                 cTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                      >  2:4760317/1‑150 (MQ=255)
                                                                                                                 cTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                      >  2:2893503/1‑150 (MQ=255)
                                                                                                                  tcgtcgAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                     >  1:3590767/1‑150 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAt                                   >  2:3544364/1‑150 (MQ=255)
                                                                                                                     tcgAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtc                                  >  1:1755685/1‑150 (MQ=255)
                                                                                                                     tcgAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtc                                  >  1:3031568/1‑150 (MQ=255)
                                                                                                                       gAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                                >  2:449477/1‑150 (MQ=255)
                                                                                                                        aGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCTGCGG‑AATGGT‑GCCGCCGATTGTCTTGGGTATTGCCTCCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTcc                               <  1:3807989/150‑1 (MQ=255)
                                                                                                                            ccaagcCAATCGGCGGC‑‑‑ACCATGCCGGCGGC‑CATGGTCGCCGCC‑ATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGCGCGCGAGGCGGGCACGGCGGCGCGGGCGCTGGGGCGGTGCTTCATCTCCGAgg                           >  2:4954395/5‑150 (MQ=21)
                                                                                                                            aCGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgg                           <  2:2112148/150‑1 (MQ=255)
                                                                                                                               cgcCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAggggg                        <  2:4179135/150‑1 (MQ=255)
                                                                                                                                gcCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                                <  1:4837417/141‑1 (MQ=255)
                                                                                                                                gcCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGc                       <  2:2861207/150‑1 (MQ=255)
                                                                                                                                  cAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa                     <  2:4610423/150‑1 (MQ=255)
                                                                                                                                    aT‑GCCGGCG‑‑GCCCGGGTCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAtt                   >  2:2629153/1‑150 (MQ=255)
                                                                                                                                     t‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTc                  >  1:4349204/1‑150 (MQ=255)
                                                                                                                                           gcggcaCCATGCCGGCGGC‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGt               >  1:3116012/6‑150 (MQ=255)
                                                                                                                                           gcg‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt             <  1:635547/150‑1 (MQ=255)
                                                                                                                                           gcg‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt             >  1:201527/1‑150 (MQ=255)
                                                                                                                                           gcg‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt             <  1:1456652/150‑1 (MQ=255)
                                                                                                                                             g‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgc           >  2:1292191/1‑150 (MQ=255)
                                                                                                                                                  cATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc        >  1:3739709/1‑150 (MQ=255)
                                                                                                                                                  cATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc        >  1:5060825/1‑150 (MQ=255)
                                                                                                                                                  cATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc        >  2:2802417/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)
                                                                                                                                                        |                                                                                                                                                       
ATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAAT‑GGCGGCG‑‑ACCATGGCCGCCGG‑CATGGT‑GCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGA  >  NC_002947/914487‑914785

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