Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,519 C→T 5.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 5.0% 390.6 / 5.7 119A323T (GCC→ACC) ‡fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (54/59);  new base T (3/3);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 123 total aligned reads displayed.

GAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGG  >  NC_000913/2503371‑2503668
                                                                                                                                                    |                                                                                                                                                     
gAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCa                                                                                                                                                      <  2:1119423/150‑1 (MQ=255)
gAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCa                                                                                                                                                      <  2:63595/150‑1 (MQ=255)
     ggCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACg                                                                                                                                                 >  1:1083208/1‑150 (MQ=255)
     ggCCAGCTGCGAGGCAGACTGCGCCTGTTGTTTTTGTCGTTTATCCGCCAGCGACTGCGCGAACTGATAATAACCGCTCACCGCGTCAGTCGTGTTAATCGACAGCACGCCGCACTGTGCGTCAAGGACTGCAGTTTGCCCGGCAGAACg                                                                                                                                                 >  1:290259/1‑150 (MQ=255)
      gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                                <  1:1698707/150‑1 (MQ=255)
      gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTCTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                                >  2:2087143/1‑150 (MQ=255)
        cAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGc                                                                                                                                              <  2:1151132/150‑1 (MQ=255)
         aGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCa                                                                                                                                             <  1:808506/150‑1 (MQ=255)
             gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaacct                                                                                                                                         <  2:2101249/150‑4 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATACCGGGCAATCGCATcc                                                                                                                                    >  1:1640690/1‑150 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                    <  2:454330/150‑1 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                    <  1:1012319/150‑1 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATcc                                                                                                                                    >  1:456391/1‑150 (MQ=255)
                    aGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAg                                                                                                                                  >  1:1739546/1‑150 (MQ=255)
                      ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCAAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                                >  2:2065655/1‑150 (MQ=255)
                       cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                               <  1:1547284/150‑1 (MQ=255)
                        ggctccTGTTGTTTTTGGCGGTTATCCGCCAGCGTCTGCGCGACCTGATACTAACCGCGCACCGCGTCATTCGGGTTAATCGCCCGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGc                                                                                                                              <  1:1048639/146‑1 (MQ=255)
                         gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                             >  2:1590011/1‑150 (MQ=255)
                         gcgcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCa                                                                                                                             >  2:1742404/1‑150 (MQ=255)
                           gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAg                                                                                                                           >  1:1477409/1‑150 (MQ=255)
                                ttgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACt                                                                                                                      <  1:1065847/150‑1 (MQ=255)
                                 tgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTc                                                                                                                     <  2:864130/150‑1 (MQ=255)
                                  gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTcc                                                                                                                    >  2:967367/1‑149 (MQ=255)
                                  gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCa                                                                                                                    >  2:677706/1‑150 (MQ=255)
                                     tttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTa                                                                                                                 >  1:568727/1‑150 (MQ=255)
                                        ggCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                              >  2:834088/1‑150 (MQ=255)
                                          cGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCggg                                                                                                            >  1:1693502/1‑150 (MQ=255)
                                             ttATCCGCCTGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATc                                                                                                         <  2:1941063/150‑1 (MQ=255)
                                              tATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCg                                                                                                        >  2:689365/1‑150 (MQ=255)
                                                tCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGcc                                                                                                      >  2:726588/1‑150 (MQ=255)
                                                  cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                                    >  1:987435/1‑150 (MQ=255)
                                                  cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGACAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATACAGCGGCAAGCCACTCAGTACCGGGATCGCCgc                                                                                                    >  2:1126037/1‑149 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAg                                                                                                   <  1:1264239/150‑1 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAg                                                                                                   <  2:539285/150‑1 (MQ=255)
                                                     cAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGc                                                                                                 <  2:158485/150‑1 (MQ=255)
                                                      aGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCa                                                                                                >  2:824373/1‑150 (MQ=255)
                                                        cGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACg                                                                                              <  2:703622/150‑1 (MQ=255)
                                                        cGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACg                                                                                              >  2:891648/1‑150 (MQ=255)
                                                         gTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACgg                                                                                             <  2:836221/150‑1 (MQ=255)
                                                            tGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                          <  2:1871400/150‑1 (MQ=255)
                                                             gcgcGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACGGGGCGTCAAGCACTGCAGGTGGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCa                                                                                         <  1:1185477/150‑1 (MQ=255)
                                                              cgcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAg                                                                                        <  2:316583/150‑1 (MQ=255)
                                                               gcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGa                                                                                       <  2:1453133/150‑1 (MQ=255)
                                                                cgACCTGATAATATCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGaa                                                                                      <  2:1392173/150‑1 (MQ=255)
                                                                 gACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAAt                                                                                     <  1:1835904/150‑1 (MQ=255)
                                                                    cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                  >  2:254645/1‑150 (MQ=255)
                                                                    cTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAg                                                                                  <  1:615188/150‑1 (MQ=255)
                                                                     tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCCCGGGCCAGCATCTgt                                                                                 >  2:348564/1‑150 (MQ=255)
                                                                     tGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                                 >  1:1250480/1‑150 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                                <  2:817032/150‑1 (MQ=255)
                                                                          ataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGCTCGCCGCGGCCCGGGCCAGAATCAGTGTATg                                                                            >  2:233515/1‑150 (MQ=255)
                                                                          ataaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATg                                                                            <  1:1904306/150‑1 (MQ=255)
                                                                           taaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCACTCGCCTCCAGCGGCACGCCACTCAGTACCGGGATCGCCGAGGCGCGGGCCAGAATCAGTGTATGc                                                                           >  2:1144419/1‑150 (MQ=255)
                                                                           taaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGc                                                                           >  2:959359/1‑150 (MQ=255)
                                                                            aaCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCCGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGCATCAGTCTATGCg                                                                          >  1:1386810/1‑150 (MQ=255)
                                                                             aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGa                                                                         <  1:553970/150‑1 (MQ=255)
                                                                             aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGa                                                                         <  1:1788422/150‑1 (MQ=255)
                                                                             aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGa                                                                         >  1:1454624/1‑150 (MQ=255)
                                                                                 cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCACTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTg                                                                     >  1:2165152/1‑150 (MQ=255)
                                                                                    aCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGACAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTgcgc                                                                  >  2:1470227/1‑150 (MQ=255)
                                                                                     ccGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATTGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGcc                                                                 >  2:1508369/1‑150 (MQ=255)
                                                                                     ccGCGTCAATCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCCTCCAGCGGCAAGCCACTCCGTACCTGGATCGCAGCGGCAAGGGGCAGAATCAGAGTATCCGCGGTGCGcc                                                                 >  1:1621821/1‑150 (MQ=255)
                                                                                       gcgTCATTCGGGTTAATCGCCAGCACGCCGGACTGTGCGGCAAGCACTGCAGTATGCCCGTCATAACGGGCAATCGCATCAAGCGGCAAGGCACTCAGTAGCGGGATCGCCGAGGCACGTGCCAGAATCAGTGTAAGCGAGGTGCGCCCg                                                               <  2:332669/150‑1 (MQ=255)
                                                                                        cgTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCgg                                                              <  1:2029696/150‑1 (MQ=255)
                                                                                        cgTCATTCGGCGCACTCGCCAGCCCGCCGCACTGGGCGGCAAGCACAGCAGGTTGCCCGGCATAACGGGCAAGCGCACCCCGCGGCAAGCCACTCCGTACCGGGCGCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCgg                                                              <  1:2027350/150‑1 (MQ=255)
                                                                                         gTCATTCGGGTTAATCGCCATCACGACGCACTGGGCGTCAAGCACTGAATGTTGCACGGCATAAAGGGCAATCGCATCCAGCGTCAAGACACTCATTACATCGATCGCCGCGGCACGGGCCAGAATCAGTTTATGCGATGTTCTACCGGt                                                             >  2:2151592/1‑150 (MQ=255)
                                                                                         gTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGt                                                             <  2:1739546/150‑1 (MQ=255)
                                                                                           cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGttt                                                           <  1:601332/150‑1 (MQ=255)
                                                                                             ttCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATGACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                         <  2:538291/150‑1 (MQ=255)
                                                                                                  gctAATCGCCAGCCCGCCGCACTGGCCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                                    <  1:292359/148‑1 (MQ=255)
                                                                                                     aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCCCGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGGTTTCTCCAGAAt                                                 >  2:1667466/1‑150 (MQ=255)
                                                                                                     aaTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAAt                                                 <  2:1420133/150‑1 (MQ=255)
                                                                                                      aTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                                >  2:1197799/1‑150 (MQ=255)
                                                                                                       tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCa                                               <  2:2154242/150‑1 (MQ=255)
                                                                                                         gCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATg                                             >  2:1868178/1‑150 (MQ=255)
                                                                                                         gCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTGTTCTCCAGAATCATg                                             >  1:1113957/1‑150 (MQ=255)
                                                                                                           cAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGcc                                           <  1:456816/150‑1 (MQ=255)
                                                                                                                 gccgcACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAg                                     <  2:456391/150‑1 (MQ=255)
                                                                                                                     cACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                                 >  2:1234691/1‑150 (MQ=255)
                                                                                                                       cTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttttt                               >  1:777552/1‑150 (MQ=255)
                                                                                                                        tGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                              <  1:1559885/150‑1 (MQ=255)
                                                                                                                         gggCGTCAAGCACTGCAGGTTGCCCGGCATACCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCa                             <  2:655526/150‑1 (MQ=255)
                                                                                                                         gggCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCa                             <  2:1783532/150‑1 (MQ=255)
                                                                                                                           gCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCaaa                           <  1:891648/150‑1 (MQ=255)
                                                                                                                               cAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCaa                       >  1:176477/1‑150 (MQ=255)
                                                                                                                                  gCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGc                    <  1:2087143/150‑1 (MQ=255)
                                                                                                                                  gCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGc                    <  2:182817/150‑1 (MQ=255)
                                                                                                                                     cTGCAGGTTGCCCGGCATAACGGGCAATCGCGTCCAGCGGCAAGCCACTCAGTACCGGTATCGCCGAGGCACGTGCCAGAATCAGTGTATGCGAGGTGCGCCCGTTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCa                 <  2:290259/150‑1 (MQ=255)
                                                                                                                                       gCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCATTGCTCaaa               >  2:623766/1‑150 (MQ=255)
                                                                                                                                       gCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaa               >  2:1759198/1‑150 (MQ=255)
                                                                                                                                        cAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaa              <  2:616733/150‑1 (MQ=255)
                                                                                                                                         cgattGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa             >  2:652368/4‑150 (MQ=255)
                                                                                                                                         cgattGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa             >  2:2169878/4‑150 (MQ=255)
                                                                                                                                          ggTTGCCCGGCATAACGGGCAATAGCATCCAGAGGCAAGACACTCATTTCCGGCATAGCCGCGTCACGGGCCTGAATCAGTGTATGCGAGGTGCGACCGGATTTCTACCGAATCATGCCCGCAACAATTTTCAAATCAAGGCTCAAAAAc            >  2:1186025/1‑150 (MQ=255)
                                                                                                                                            ttGCCCGGCATAACGGGAAATCGAATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTg          <  2:2056401/150‑1 (MQ=255)
                                                                                                                                               cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCt       >  2:1912027/1‑150 (MQ=255)
                                                                                                                                                   gCATAACGGGCAATCGCTTCCAGGGGCAATCCACTCAGTACTGGGATAGCCGAGGCACGGGCCATAATCAGTGTATGCGAGGTTCGCCCGGTTTTCTCCATAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAACACTGGCTTggg   <  2:126934/150‑1 (MQ=255)
                                                                                                                                                    cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg  <  2:1693502/150‑1 (MQ=255)
                                                                                                                                                    cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg  <  1:1217734/150‑1 (MQ=255)
                                                                                                                                                    cATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTgggg  <  1:824373/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGG  >  NC_000913/2503371‑2503668

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