Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,519 C→T 7.0% A323T (GCC→ACC) ‡ fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,503,5190CT 7.0% 356.2 / 10.8 114A323T (GCC→ACC) ‡fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (51/54);  new base T (4/4);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 114 total aligned reads displayed.

GAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCT  >  NC_000913/2503371‑2503663
                                                                                                                                                    |                                                                                                                                                
gAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCTCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTCa                                                                                                                                                 <  2:1717137/150‑1 (MQ=255)
gAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCa                                                                                                                                                 <  1:613289/150‑1 (MQ=255)
     ggCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGCCCTGATAATCACCGCTCACCGCGTCAGTCGGGTTAATCGCCCGCTCGCCGCACTGGGCGTGAAGCACTGCAGGTTGCGCGGCGTAACg                                                                                                                                            >  2:89603/1‑150 (MQ=255)
     ggCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTCATAACg                                                                                                                                            >  2:1198096/1‑150 (MQ=255)
     ggCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACg                                                                                                                                            >  1:1656504/1‑150 (MQ=255)
      gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                           >  2:1570423/1‑150 (MQ=255)
          gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCaa                                                                                                                                       <  1:1374178/150‑1 (MQ=255)
           cTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGGCCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                      <  1:424046/150‑2 (MQ=255)
           cTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATACCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCCCTGGGCGTCAAGCACTGCAGGTTGCCCGGCCTCCCGGGCact                                                                                                                                      >  2:888122/1‑148 (MQ=255)
           cTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCCCTGCAGGTTGCCCGGCATAACGGGCAAt                                                                                                                                      >  2:882168/1‑150 (MQ=255)
            ggcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaacc                                                                                                                                     <  1:1503892/149‑3 (MQ=255)
             gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCg                                                                                                                                    <  1:1174057/150‑1 (MQ=255)
              cgcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaacctg                                                                                                                                   <  2:1798951/150‑5 (MQ=255)
               gcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaacctgc                                                                                                                                  <  1:1254811/150‑6 (MQ=255)
               gcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCa                                                                                                                                  <  1:1556874/150‑1 (MQ=255)
                cgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCAt                                                                                                                                 <  1:401478/150‑1 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                               <  2:332964/150‑1 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                               <  2:1756457/150‑1 (MQ=255)
                   cAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCa                                                                                                                              <  1:1744408/150‑1 (MQ=255)
                      ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                           >  1:847036/1‑150 (MQ=255)
                       cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                          <  2:77132/150‑1 (MQ=255)
                        tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGc                                                                                                                         <  2:1436628/150‑1 (MQ=255)
                        tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGc                                                                                                                         >  2:568568/1‑150 (MQ=255)
                        tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGGCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGc                                                                                                                         <  2:734282/150‑1 (MQ=255)
                           gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAg                                                                                                                      <  1:1194850/150‑1 (MQ=255)
                             cTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGcc                                                                                                                    <  2:222033/150‑1 (MQ=255)
                             cTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGcc                                                                                                                    >  1:128405/1‑150 (MQ=255)
                              tgttgtTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTAACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCa                                                                                                                   >  1:1754739/1‑150 (MQ=255)
                                  gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCa                                                                                                               <  2:292687/150‑1 (MQ=255)
                                         gCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCgg                                                                                                        >  1:172233/1‑150 (MQ=255)
                                          cGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCggg                                                                                                       >  2:22167/1‑150 (MQ=255)
                                          cGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCggg                                                                                                       <  2:1080848/150‑1 (MQ=255)
                                           gTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGa                                                                                                      <  2:1553165/150‑1 (MQ=255)
                                            tttATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAt                                                                                                     >  1:227800/1‑150 (MQ=255)
                                               aTCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGc                                                                                                  >  2:1513896/1‑150 (MQ=255)
                                                 ccGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTCCCCGGCATAACGGGCAATCGCCTCAAGCGGCAAGCCACTCAGTACCGGGATCGCCg                                                                                                >  2:255043/1‑150 (MQ=255)
                                                  cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCCTCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                               >  2:1402151/1‑150 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAg                                                                                              <  2:1300705/150‑1 (MQ=255)
                                                    ccAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTCATAACGGGCAATCGAATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAgg                                                                                             <  2:1728167/150‑1 (MQ=255)
                                                    ccAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAgg                                                                                             <  1:669492/150‑1 (MQ=255)
                                                          tCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACggg                                                                                       <  1:444354/150‑1 (MQ=255)
                                                            tGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                     <  1:988162/150‑1 (MQ=255)
                                                            tGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                     >  2:292580/1‑150 (MQ=255)
                                                               gcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGa                                                                                  <  1:1010116/150‑1 (MQ=255)
                                                                 gACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAAt                                                                                >  2:514453/1‑150 (MQ=255)
                                                                   ccTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCa                                                                              >  1:1724224/1‑150 (MQ=255)
                                                                    cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                             <  1:625195/150‑1 (MQ=255)
                                                                     tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                            >  1:683692/1‑150 (MQ=255)
                                                                     tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                            >  1:1719026/1‑150 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCACTCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                           >  1:196500/1‑150 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                           <  1:1570423/150‑1 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                           >  2:1774007/1‑150 (MQ=255)
                                                                        taataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTa                                                                         <  2:230669/150‑1 (MQ=255)
                                                                          ataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATg                                                                       >  1:1363957/1‑150 (MQ=255)
                                                                             aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGa                                                                    >  1:1691454/1‑150 (MQ=255)
                                                                                 cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTg                                                                <  1:1826383/150‑1 (MQ=255)
                                                                                     ccGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGcc                                                            <  1:1821967/150‑1 (MQ=255)
                                                                                     ccGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCAAGAATCAGTGTATGCGAGGTGCGcc                                                            >  2:690539/1‑150 (MQ=255)
                                                                                      cgcgTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGccc                                                           >  1:1669167/1‑150 (MQ=255)
                                                                                         gTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGt                                                        <  1:1520198/150‑1 (MQ=255)
                                                                                           cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGttt                                                      <  2:1656504/150‑1 (MQ=255)
                                                                                             ttCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                    <  2:823397/150‑1 (MQ=255)
                                                                                             ttCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGCGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                    <  1:1161003/150‑1 (MQ=255)
                                                                                                gggTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTcc                                                 <  1:1465607/150‑1 (MQ=255)
                                                                                                 ggTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCa                                                >  2:891551/1‑150 (MQ=255)
                                                                                                 ggTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGAACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCa                                                <  2:1100792/150‑1 (MQ=255)
                                                                                                  gTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTTCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                               <  2:646421/150‑1 (MQ=255)
                                                                                                       tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCa                                          >  1:1310222/1‑150 (MQ=255)
                                                                                                             gCACGCCGCACTGGGCGTCAAGCACTGCAGGTTTCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCAGGTTTTCTCCAGAATCATGCCCg                                    >  2:1315168/1‑150 (MQ=255)
                                                                                                              cACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGc                                   <  2:227800/150‑1 (MQ=255)
                                                                                                               aCGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGCTCGCCGATGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCa                                  <  2:1104362/150‑1 (MQ=255)
                                                                                                                cgccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCCCGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCaa                                 >  1:951166/1‑150 (MQ=255)
                                                                                                                cgccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCaa                                 >  2:486818/1‑150 (MQ=255)
                                                                                                                  ccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGATCGCCGAGGCACGGGCAAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGa                               >  1:1351597/1‑150 (MQ=255)
                                                                                                                  ccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGa                               >  2:563143/1‑150 (MQ=255)
                                                                                                                      ccTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAtttt                           <  1:234777/149‑1 (MQ=255)
                                                                                                                      aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGTCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATTCGAGGGGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAtttt                           <  1:1420356/150‑1 (MQ=255)
                                                                                                                      aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAtttt                           <  2:1363957/150‑1 (MQ=255)
                                                                                                                        tGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                         <  1:1480444/150‑1 (MQ=255)
                                                                                                                        tGGGCGTCAAGCACTGCAGGTGGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                         <  2:297867/150‑1 (MQ=255)
                                                                                                                        ggggCGTCAAGCACGGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACGCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                         <  1:888122/149‑1 (MQ=255)
                                                                                                                        ggggCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                         <  1:882168/149‑1 (MQ=255)
                                                                                                                             gTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATc                    >  2:1823382/1‑150 (MQ=255)
                                                                                                                             gTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATc                    <  1:1795697/150‑1 (MQ=255)
                                                                                                                             gTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATc                    >  2:1294712/1‑150 (MQ=255)
                                                                                                                             gTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATGTTTCAAATc                    >  2:901126/1‑150 (MQ=255)
                                                                                                                                aaGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAg                 >  1:1536436/1‑150 (MQ=255)
                                                                                                                                aaGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACATAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGAGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAg                 <  2:1729061/150‑1 (MQ=255)
                                                                                                                                aaGCACTGCAGGTTGCCCGGCATAACGGACAATCGCATCCAGCGGCAAGCCACGCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAg                 <  2:891232/150‑1 (MQ=255)
                                                                                                                                 aGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAgg                <  1:568568/150‑1 (MQ=255)
                                                                                                                                    aCTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTc             <  2:1765480/150‑1 (MQ=255)
                                                                                                                                     cTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCa            >  1:692230/1‑150 (MQ=255)
                                                                                                                                      tGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaa           >  1:1291900/1‑150 (MQ=255)
                                                                                                                                       gCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaa          <  2:683692/150‑1 (MQ=255)
                                                                                                                                        gcgattGCCCGTTATGCCGGGCAACCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaa         >  1:1069760/5‑150 (MQ=255)
                                                                                                                                         cgattGCCCGTTATGCCGGGCAACCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTATCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa        >  1:1632355/4‑150 (MQ=255)
                                                                                                                                          ggTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc       >  1:540798/1‑150 (MQ=255)
                                                                                                                                           gNTGCCCGGCATAACGAGCAATCGCATCCAGCGCCAAGCCACTCAGTACCGGGATCGACTAGGCAGGGGCCAGAATGAGTGAACGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGAAAGATATTTCACATCAATGCTCATAAACt      >  2:331209/1‑150 (MQ=255)
                                                                                                                                           attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGTGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt      >  2:1465212/2‑150 (MQ=255)
                                                                                                                                               cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCt  >  1:1310430/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCT  >  NC_000913/2503371‑2503663

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑