Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,503,519 C→T 5.7% 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.7% 307.7 / 6.5 106A323T (GCC→ACC) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (44/55);  new base T (2/4);  total (47/59)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

Only 100 of 107 total aligned reads displayed.

GAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCT  >  NC_000913/2503371‑2503663
                                                                                                                                                    |                                                                                                                                                
gAATAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGc                                                                                                                                                  <  2:1139075/149‑1 (MQ=255)
  aTAGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTAGCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCTTCATTCGGGTTACTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGATGCCGGGCAt                                                                                                                                                >  1:1064733/1‑149 (MQ=255)
    aGGCCAGCTGCTCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCGTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATaa                                                                                                                                              >  2:1347999/1‑149 (MQ=255)
    aGGCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATaa                                                                                                                                              >  2:360433/1‑149 (MQ=255)
      gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATACCg                                                                                                                                            >  1:221557/1‑149 (MQ=255)
      gCCAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACg                                                                                                                                            <  2:149544/149‑1 (MQ=255)
       ccAGCTGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCCCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACgg                                                                                                                                           >  1:796027/1‑149 (MQ=255)
            tgagcgTCAGCCTGCGCCGGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                      <  2:247357/146‑2 (MQ=255)
            tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTCTCCGCCAGCGTCTGCGCGACCTGATAA‑AACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                     <  2:519013/149‑1 (MQ=255)
            tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                      <  1:1449303/149‑2 (MQ=255)
            tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                      <  1:509328/149‑2 (MQ=255)
            tGCGCGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCaag                                                                                                                                      >  1:718821/1‑148 (MQ=255)
            cgcgcGGCAGCCTGCGCCTGTTGTTGTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGCTAATCGCCAGCACGCCGCACTGGGCGTCTAGCACTGCAGGTTGCCCGTTATGCCGGGCaac                                                                                                                                      <  1:1066686/148‑2 (MQ=255)
             gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACAGGCAATc                                                                                                                                     <  1:785460/149‑1 (MQ=255)
             gcgcGGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGCTTGCCCGGCATAACGGGCAATc                                                                                                                                     >  1:1149816/1‑149 (MQ=255)
             gcgcGGCAGCCTGCGCAGGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACAGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGAGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATc                                                                                                                                     <  1:69743/149‑1 (MQ=255)
             gcgcAGCAGCCTGCGCCTGTTGTTTTTGGCGTGTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATc                                                                                                                                     <  1:877435/149‑1 (MQ=255)
               tcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACATGATAATAACCTCTAACCGCGTCATTCGGGTTAATCGCAAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGAATAACGTGCAATctc                                                                                                                                   >  2:576656/2‑147 (MQ=255)
               gcgGCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGc                                                                                                                                   >  2:570190/1‑149 (MQ=255)
                 ggCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCAt                                                                                                                                 <  1:544393/149‑1 (MQ=255)
                 ggCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCAt                                                                                                                                 >  1:1250175/1‑149 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCCTTCGGGTTAATCGCCAGCACGCCGCACTGAGCGTCAAGCACTGCAGGTTGCCCGGCAGAACGGGCAATCGCATc                                                                                                                                >  1:1391493/1‑149 (MQ=255)
                  gCAGCCTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATc                                                                                                                                <  2:808538/149‑1 (MQ=255)
                      ccTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGc                                                                                                                            <  1:61256/149‑1 (MQ=255)
                       cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                           <  2:561106/149‑1 (MQ=255)
                       cTGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCg                                                                                                                           <  2:552123/149‑1 (MQ=255)
                        tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                          <  2:81184/149‑1 (MQ=255)
                        tGCGCCTGTTGTTTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCgg                                                                                                                          <  1:758295/149‑1 (MQ=255)
                           gcCTGTTGTTTTTGGCGTTTATCCGCCAGCGGCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCaa                                                                                                                       <  1:1565683/149‑1 (MQ=255)
                               gttgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCCTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCAGAACGGGCATTCGCATCCAGCGGCAAGCCa                                                                                                                   >  2:1014504/1‑149 (MQ=255)
                                 tgttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACt                                                                                                                 <  1:11685/149‑1 (MQ=255)
                                  gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTc                                                                                                                >  2:920887/1‑149 (MQ=255)
                                  gttTTTGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTc                                                                                                                <  2:499636/149‑1 (MQ=255)
                                   tttttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGTGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCATTCGCATCCAGCGGCAAGACCCTcc                                                                                                               >  1:863501/1‑148 (MQ=255)
                                   tttttGGCGTTTATCCGCCAGCGTCTGCGCGACATGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCCCACTGGGCGTCAAGAACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACtaa                                                                                                               >  2:1186847/1‑147 (MQ=255)
                                     tttGGCGTTTAGCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAAGCGCATCCAGCGGCAAGCCACTCAGt                                                                                                             <  2:530280/149‑1 (MQ=255)
                                      ttGGCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTa                                                                                                            <  1:1146283/149‑1 (MQ=255)
                                        ggCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTAcc                                                                                                          >  2:1475496/1‑149 (MQ=255)
                                        ggCGTTAATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTAcc                                                                                                          <  1:292065/149‑1 (MQ=255)
                                         gCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                         <  2:354465/149‑1 (MQ=255)
                                         gCGTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCg                                                                                                         <  1:1068919/149‑1 (MQ=255)
                                           gTTTATCCGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCggg                                                                                                       >  2:579761/1‑149 (MQ=255)
                                             ttATCCGCCAGCGTCTGCGCGACCTGATAATACCCGCTCCCCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGAt                                                                                                     <  2:856647/149‑1 (MQ=255)
                                                  ctccACCGTCTGCTCGACCTGATAATAACCGCTCACCGAGTCATTCGGGTTAATAGCAAGCACGCTGCACTGGTCGACAAGCACTTCACGTTGCACGCCATAACGTGAAATAGCATCCAGCCTCAAGCCACACAGTACCGGGATCGCCg                                                                                                >  2:687065/3‑149 (MQ=255)
                                                  cGCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGACCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCg                                                                                                >  2:1350803/1‑149 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                               <  1:125559/149‑1 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                               <  1:516050/149‑1 (MQ=255)
                                                   gCCAGCGTCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGa                                                                                               >  2:31343/1‑149 (MQ=255)
                                                    ccaggttCTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGTGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCagag                                                                                              <  2:126466/143‑1 (MQ=255)
                                                           cTGCGCGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACggg                                                                                       <  1:1475496/149‑1 (MQ=255)
                                                             gcgcGACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGcc                                                                                     <  1:1332147/149‑1 (MQ=255)
                                                              cgcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCa                                                                                    >  1:524360/1‑149 (MQ=255)
                                                               gcgACCTGATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCAAGCGGCAAGCCACTCAGAACCGGGACCGCCGAGGCACGGGCCAg                                                                                   >  1:1499528/1‑149 (MQ=255)
                                                                  cccTGATAAAAACCGCTCACCGCCGCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGCTTGCCCGGCATACCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTCCCGGGATCGCCGAGGCACGGGCCAGAAt                                                                                <  2:963132/148‑1 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                            >  2:1503499/1‑149 (MQ=255)
                                                                      gATAATAACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTCGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgt                                                                            <  2:390690/149‑1 (MQ=255)
                                                                       ataataACCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAgtg                                                                           >  1:425498/1‑149 (MQ=255)
                                                                          ataaCCGCTCACCGAGTCATTCGAGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGATATCCCGGCATAACGGGCAATCGCATCCAGCGGAATGCCACTAAGTACCGGGATAGCCGAGGCACGCGCCAGAATCAGAGTAt                                                                        >  2:109769/1‑149 (MQ=255)
                                                                             aCCGCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCg                                                                     >  1:378951/1‑149 (MQ=255)
                                                                               cGCTCACCGCGGCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAg                                                                   <  1:301403/149‑1 (MQ=255)
                                                                                gCTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAgg                                                                  <  1:920887/149‑1 (MQ=255)
                                                                                 cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                 <  1:553204/149‑1 (MQ=255)
                                                                                 cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                 <  1:668527/149‑1 (MQ=255)
                                                                                 cTCACCGCGTCATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGt                                                                 <  2:776812/149‑1 (MQ=255)
                                                                                           cATTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAAGGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtt                                                       <  1:1347999/149‑1 (MQ=255)
                                                                                            aTTCGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCAAGAATCAGTGGATGCGCGGTGCGCCCGGttg                                                      >  2:465851/1‑148 (MQ=255)
                                                                                             ctCGGGTCACTCGCCAGCACGCCGCACTGGGCGGCAAGCACTGCAGGTTGCCCGGCATAACCGGCAATCGCATCCAGCGGCAAGCCACTCAGTCCCGGGAGCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGtttt                                                     <  1:207446/148‑1 (MQ=255)
                                                                                              ccGGGTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTtc                                                    <  2:728999/148‑1 (MQ=255)
                                                                                                 ggTTCATCGCCAGCACGCCGCACTGGGCGTCAAGCACGGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGAACCGGGATCGCCGAGGCACGGGCCAGAATCAGGGTATGCGAGGTGCGCCCGGGTTTCTcc                                                 <  1:227043/149‑1 (MQ=255)
                                                                                                 ggTTAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTcc                                                 >  2:1108980/1‑149 (MQ=255)
                                                                                                   ttACTCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                               <  1:1014504/149‑1 (MQ=255)
                                                                                                   ttAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGGTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAg                                               <  1:1493759/149‑1 (MQ=255)
                                                                                                    tAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGa                                              <  1:1119267/149‑1 (MQ=255)
                                                                                                    tAATCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTATTCTCCAGa                                              >  2:469567/1‑149 (MQ=255)
                                                                                                       tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                           <  2:341491/149‑1 (MQ=255)
                                                                                                       tCGCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATc                                           <  1:23311/149‑1 (MQ=255)
                                                                                                         gCCAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCAt                                         >  1:1224460/1‑149 (MQ=255)
                                                                                                          ccAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGATCGCCGAGGCACGGGCCAGATTCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATg                                        >  1:1352567/1‑149 (MQ=255)
                                                                                                           cAGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                       <  1:516686/149‑1 (MQ=255)
                                                                                                           cAGCACGCCGCACTGCGCGTCAAGCACTGCAGGTGGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGc                                       <  2:1391493/149‑1 (MQ=255)
                                                                                                            aGCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGcc                                      <  2:1250175/149‑1 (MQ=255)
                                                                                                             gCACGCCGCACTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGccc                                     >  1:444243/1‑149 (MQ=255)
                                                                                                                      aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                            >  1:770096/1‑149 (MQ=255)
                                                                                                                      aCTGGGCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGAttt                            >  2:12302/1‑149 (MQ=255)
                                                                                                                         gggCGTCAAGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTc                         <  1:826272/149‑1 (MQ=255)
                                                                                                                                aaGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCaa                  >  2:1139434/1‑149 (MQ=255)
                                                                                                                                 aGCACTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCTGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCacg                 >  1:248121/1‑147 (MQ=255)
                                                                                                                                    aCTGCAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCGGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCt              >  2:919785/1‑149 (MQ=255)
                                                                                                                                        cAGGTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACGGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGGGTTCTCCAGAATCCTGCCCGCCAGCTTTTTCATATCAAGGCTCcaa          >  1:394342/1‑149 (MQ=255)
                                                                                                                                         aGGTTGCCCGGCAAAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTACAGAATCATGCCCGCAAGATTTTTCATATCAAGTCTCtaaa         <  2:1214964/149‑1 (MQ=255)
                                                                                                                                          ggTTGGCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaa        <  1:570190/149‑1 (MQ=255)
                                                                                                                                          ggTTGCCCGGCATAACGGGCAATCGCATCCATCGGCAAGCCACTCAGTAGCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGACCGCAAGATTTTTCAAATCAAGGCTCaaaaa        <  1:576656/149‑1 (MQ=255)
                                                                                                                                           gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCCCCGTACCGGGATCGCCGAGGCACGTGCCAGAATCAGTGTATGCGAGGTGCGCACGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATAAAGGCTCAAAAAc       <  2:1064733/149‑1 (MQ=255)
                                                                                                                                           gTTGCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc       >  1:1491996/1‑149 (MQ=255)
                                                                                                                                           gTTGCCAGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGACAGAATCAGTGTATGCGAGGTGCGCCAGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAAc       >  2:837441/1‑149 (MQ=255)
                                                                                                                                           attGCCCGTTATGCCGGGCAATCGCATCCACCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCACGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCaaaaaa       >  2:588456/2‑148 (MQ=255)
                                                                                                                                            ttGCCCGTTATGCCGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCCAGGCTCAAAAACt      >  1:836142/1‑149 (MQ=255)
                                                                                                                                              gCCCGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTgg    >  1:1553908/1‑149 (MQ=255)
                                                                                                                                               cccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAACGCTCAAAAACTGGc   <  2:378951/149‑1 (MQ=255)
                                                                                                                                                ccGGCATAACGGGCAATCGCATCCAGCGGCAAGCCCCTCAGTACCGGGCTCGCCGAGGCACGGGCCAGCATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCCAATCACGGCTCAAAAACTGGCt  >  1:989389/1‑149 (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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑