Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,523 A→C 9.4% R321R (CGT→CGG 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,5230AC 9.4% 514.4 / 33.1 161R321R (CGT→CGGfryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (56/89);  new base C (6/10);  total (62/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAA  >  NC_000913/2503377‑2503672
                                                                                                                                                  |                                                                                                                                                     
gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                                    >  2:395979/1‑150 (MQ=255)
 ccAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACggg                                                                                                                                                   >  1:1651141/1‑150 (MQ=255)
 ccAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACggg                                                                                                                                                   >  1:1821782/1‑150 (MQ=255)
 ccAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACggg                                                                                                                                                   >  2:2032484/1‑150 (MQ=255)
   aGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCa                                                                                                                                                 >  1:2256273/1‑150 (MQ=255)
    gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaa                                                                                                                                                <  2:1481290/150‑1 (MQ=255)
    gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaa                                                                                                                                                <  2:2156403/150‑1 (MQ=255)
    gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaa                                                                                                                                                <  2:804567/150‑1 (MQ=255)
    gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaa                                                                                                                                                <  2:804705/150‑1 (MQ=255)
    gCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaa                                                                                                                                                <  2:1758070/150‑1 (MQ=255)
      tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                              <  1:1016561/150‑1 (MQ=255)
              aGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAg                                                                                                                                      <  2:124519/150‑1 (MQ=255)
                 cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCTCGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                                   >  2:1827569/1‑150 (MQ=255)
                   gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                                 >  2:2000190/1‑150 (MQ=255)
                   gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                                 >  2:61259/1‑150 (MQ=255)
                     gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAg                                                                                                                               <  1:2207184/150‑1 (MQ=255)
                     gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAg                                                                                                                               <  2:1993244/150‑1 (MQ=255)
                        tgttgtTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCa                                                                                                                            >  2:363546/1‑150 (MQ=255)
                        tgttgtTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCa                                                                                                                            <  1:253724/150‑1 (MQ=255)
                        tgttgtTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCa                                                                                                                            <  1:678947/150‑1 (MQ=255)
                                     gTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGa                                                                                                               >  2:1673561/1‑150 (MQ=255)
                                      tttATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAt                                                                                                              <  1:1350849/150‑1 (MQ=255)
                                      tttATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAt                                                                                                              <  1:1160412/150‑1 (MQ=255)
                                          tCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGcc                                                                                                          <  2:2376995/150‑1 (MQ=255)
                                            cGCCAGCGTCTGCGCGTCCTGATAATAACCGCTCACCGCGTCTTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCTGGTTGCCAGGCATACCGGGCATTCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                        <  2:69070/150‑1 (MQ=255)
                                                 gCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCAc                                                                                                   <  1:785805/150‑1 (MQ=255)
                                                 gCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCAc                                                                                                   <  1:208259/150‑1 (MQ=255)
                                                   gTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACgg                                                                                                 >  2:1166956/1‑150 (MQ=255)
                                                    tCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACggg                                                                                                <  2:1234207/150‑1 (MQ=255)
                                                              cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                      <  1:1517214/150‑1 (MQ=255)
                                                              cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                      <  2:1261583/150‑1 (MQ=255)
                                                              cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                      <  2:1341223/150‑1 (MQ=255)
                                                               tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                                     >  2:2105688/1‑150 (MQ=255)
                                                                            tCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTgc                                                                        <  1:30462/150‑1 (MQ=255)
                                                                            tCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTgc                                                                        >  1:1027208/1‑150 (MQ=255)
                                                                                 gcgTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCg                                                                   <  2:1692809/150‑1 (MQ=255)
                                                                                    tCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtt                                                                <  2:670297/150‑1 (MQ=255)
                                                                                      aTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtttt                                                              >  2:2259385/1‑150 (MQ=255)
                                                                                      aTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtttt                                                              >  1:2063166/1‑150 (MQ=255)
                                                                                        tCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCGAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtct                                                            <  2:2080047/150‑1 (MQ=255)
                                                                                         cgcgatAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTCCCCGGCATACCGGGCAATCGCAACCACCCCCAAGCCACTCAGTACCGGGAACGCCGAGGCACGGGACAGAAACAGTCTATGCGAGCTGCGCCCGGTTTtctc                                                           <  2:796467/145‑1 (MQ=255)
                                                                                         cGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtctc                                                           <  1:1031964/150‑1 (MQ=255)
                                                                                         cGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtctc                                                           <  1:2226410/150‑1 (MQ=255)
                                                                                         cGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtctc                                                           <  2:2108293/150‑1 (MQ=255)
                                                                                               aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAAt                                                     >  2:1387192/1‑150 (MQ=255)
                                                                                                  cGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCAt                                                  <  1:702930/150‑1 (MQ=255)
                                                                                                  cGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCAt                                                  <  2:1998888/150‑1 (MQ=255)
                                                                                                    ccAGCACGCCGCACTGGGCGTCCAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                                >  2:2360438/1‑150 (MQ=255)
                                                                                                    ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                                >  2:805890/1‑150 (MQ=255)
                                                                                                    ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                                >  2:1097219/1‑150 (MQ=255)
                                                                                                    ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                                >  2:409084/1‑150 (MQ=255)
                                                                                                    ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                                >  2:451745/1‑150 (MQ=255)
                                                                                                     cAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGcc                                               >  2:326163/1‑150 (MQ=255)
                                                                                                      aGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGccc                                              >  1:1105014/1‑150 (MQ=255)
                                                                                                        cACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGc                                            >  1:1723617/1‑150 (MQ=255)
                                                                                                          cgccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCaa                                          <  1:1889424/150‑1 (MQ=255)
                                                                                                          cgccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCaa                                          <  2:1608073/150‑1 (MQ=255)
                                                                                                          cgccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCaa                                          <  1:23441/150‑1 (MQ=255)
                                                                                                                  tGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                                  >  1:984707/1‑150 (MQ=255)
                                                                                                                   gggCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCa                                 <  2:286122/150‑1 (MQ=255)
                                                                                                                     gCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCaaa                               <  1:2086964/150‑1 (MQ=255)
                                                                                                                        tCAAGCACTGCAGGTTCCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCCGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCa                            >  2:585328/1‑150 (MQ=255)
                                                                                                                         cAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCaa                           <  1:2285400/150‑1 (MQ=255)
                                                                                                                          aaGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                                  <  2:774433/142‑1 (MQ=255)
                                                                                                                           aGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAgg                         >  1:1896186/1‑150 (MQ=255)
                                                                                                                               cTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCa                     <  1:1974310/150‑1 (MQ=255)
                                                                                                                                tGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaa                    <  2:2318548/150‑1 (MQ=255)
                                                                                                                                tGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaa                    <  2:909716/150‑1 (MQ=255)
                                                                                                                                   aGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa                 >  1:2158508/1‑150 (MQ=255)
                                                                                                                                     gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  2:1504831/1‑150 (MQ=255)
                                                                                                                                     gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               <  2:743548/150‑1 (MQ=255)
                                                                                                                                     gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               <  2:344262/150‑1 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCCAAAACt               >  2:1713968/2‑150 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  2:1719145/2‑150 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  2:1648212/2‑150 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  1:961536/2‑150 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  2:2196419/2‑150 (MQ=255)
                                                                                                                                     attGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACt               >  2:241972/2‑150 (MQ=255)
                                                                                                                                      ttGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTg              <  1:1304850/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:2000190/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:1600465/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:1600498/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:441071/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:870553/150‑1 (MQ=255)
                                                                                                                                       tGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg             <  1:61259/150‑1 (MQ=255)
                                                                                                                                         cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCt           >  1:12158/1‑150 (MQ=255)
                                                                                                                                         cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCt           >  1:12230/1‑150 (MQ=255)
                                                                                                                                         cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCt           >  2:1785324/1‑150 (MQ=255)
                                                                                                                                          ccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCtt          >  1:307330/1‑150 (MQ=255)
                                                                                                                                            ggCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGTATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgg        <  1:1437325/150‑1 (MQ=255)
                                                                                                                                            ggCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgg        <  1:1437337/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:1712891/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:1812053/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:636339/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:1567480/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:1317552/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  2:118668/150‑1 (MQ=255)
                                                                                                                                             gCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTggg       <  1:1272295/150‑1 (MQ=255)
                                                                                                                                                 aaCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTa   >  1:1991050/1‑150 (MQ=255)
                                                                                                                                                  aCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTaa  >  2:721419/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAA  >  NC_000913/2503377‑2503672

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