Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,579,157 G→T 15.6% R164R (CGG→CGT)  talA → transaldolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,579,1570GT15.6% 467.4 / 72.3 181R164R (CGG→CGT) talAtransaldolase A
Reads supporting (aligned to +/- strand):  ref base G (78/73);  new base T (14/14);  total (93/88)
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.19e-01

Only 100 of 184 total aligned reads displayed.

CTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGA  >  NC_000913/2579059‑2579255
                                                                                                  |                                                                                                  
  ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACgggac                                                                                                 >  1:2672086/1‑98 (MQ=255)
  ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGcc                                                                                                 <  2:293752/100‑1 (MQ=255)
  ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGcc                                                                                                 >  1:4244947/1‑100 (MQ=255)
  ggCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGc                                                                                                  >  2:1904407/1‑99 (MQ=255)
      tCGACCTGGGAAGGAATTAGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCtgtg                                                                                             >  2:818163/1‑100 (MQ=255)
        gACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCg                                                                                           <  2:1641219/100‑1 (MQ=255)
        gACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCg                                                                                           <  2:3999509/100‑1 (MQ=255)
          ccTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAACAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                         <  2:1968140/100‑4 (MQ=255)
          ccTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaa                                                                                         <  2:1123361/100‑4 (MQ=38)
          ccTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGa                                                                                         <  2:1286793/100‑1 (MQ=255)
           cTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGaa                                                                                        >  1:125638/1‑100 (MQ=255)
           cTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGaa                                                                                        <  2:4136568/100‑1 (MQ=255)
            tGGGAAGGAATTCGCGCGGCAGAAGAGCTGGGTAAAGAAGGTATTAACTGCAACCTGGCGCTGCTGTTGTCTTTTGCACAGGCACGGGCCTGTGCGGAAg                                                                                       <  1:3981831/100‑1 (MQ=255)
            tGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaag                                                                                       <  2:144970/100‑6 (MQ=38)
             gggAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaaga                                                                                      <  2:3356774/100‑7 (MQ=38)
             gggAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGc                                                                                      >  2:1023455/1‑100 (MQ=255)
             gggAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGc                                                                                      <  1:140240/100‑1 (MQ=255)
              ggaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCa                                                                                     <  1:3023966/100‑1 (MQ=255)
               gaaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCCCGTGCCTGTGcaaaagaaa                                                                                    <  1:3376034/100‑9 (MQ=37)
                aaggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAgg                                                                                   <  1:3074626/100‑1 (MQ=255)
                  ggaaTTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCg                                                                                 >  2:34653/1‑100 (MQ=255)
                     nttCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGttt                                                                              >  1:836507/2‑100 (MQ=255)
                      ttCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGtttt                                                                             <  2:2543257/100‑1 (MQ=255)
                      ttCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGtttt                                                                             <  2:1991012/100‑1 (MQ=255)
                        cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTc                                                                           <  1:1071519/100‑1 (MQ=255)
                        cgcgcgGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTc                                                                           <  1:580182/100‑1 (MQ=255)
                             ggCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAtt                                                                      >  1:2636158/1‑100 (MQ=255)
                             ggCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGAtt                                                                      <  2:2590274/100‑1 (MQ=255)
                               cAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTc                                                                    >  2:3817238/1‑100 (MQ=255)
                                agaagaGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCg                                                                   <  2:462787/100‑1 (MQ=255)
                                 gaagaGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCg                                                                   <  1:3364956/99‑1 (MQ=255)
                                  aagaGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGcc                                                                 >  2:3908749/1‑100 (MQ=255)
                                   agagCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCg                                                                <  2:1909239/100‑1 (MQ=255)
                                       cTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTg                                                            >  1:1804908/1‑100 (MQ=255)
                                        tGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCTTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGt                                                           >  1:1526222/1‑100 (MQ=255)
                                        tGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGt                                                           >  1:3842883/1‑100 (MQ=255)
                                        tGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGt                                                           >  2:206972/1‑100 (MQ=255)
                                         ggAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTc                                                          >  2:205738/1‑100 (MQ=255)
                                         ggAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTc                                                          >  1:1203032/1‑100 (MQ=255)
                                           aaaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCgg                                                        <  2:2573499/100‑1 (MQ=255)
                                           aaaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCg                                                         >  1:1037581/1‑99 (MQ=255)
                                            aaaaGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCggg                                                       >  2:3370087/1‑100 (MQ=255)
                                                  agGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTAtt                                                 >  1:1504740/1‑100 (MQ=255)
                                                     tATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTAt                                              >  1:3634992/1‑100 (MQ=255)
                                                      aTTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATg                                             >  2:1522633/1‑100 (MQ=255)
                                                      aTTAACTGCAACCGGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGATGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATg                                             >  1:2727698/1‑100 (MQ=255)
                                                        tAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAc                                           >  2:513709/1‑100 (MQ=255)
                                                        tAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAc                                           >  2:156762/1‑100 (MQ=255)
                                                         aaCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACt                                          >  2:698119/1‑100 (MQ=255)
                                                         aaCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACt                                          >  2:2959532/1‑100 (MQ=255)
                                                          aCTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTg                                         <  1:4068417/100‑1 (MQ=255)
                                                           cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTgg                                        <  1:1622914/100‑1 (MQ=255)
                                                           cTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCATGTGCGGAAGCAGGCGTTTTTCGGATTTCGACGTTTGTCGGGGGTATTGATGACTgg                                        >  2:3741366/1‑100 (MQ=255)
                                                             gCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGGTTTCGCCGTTTGTCGGGCGTATTTATGACTGGt                                       >  1:2311908/1‑99 (MQ=255)
                                                               aaCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCCGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATc                                    <  1:3507320/100‑1 (MQ=255)
                                                               aaCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATc                                    >  2:3425202/1‑100 (MQ=255)
                                                                aCCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCa                                   <  1:3359268/100‑1 (MQ=255)
                                                                  cTGACGCTGCTGTTTTCTTTTGTTTATTTACGAGCCTGAGCGGAAGTAGGCGTTTTTCTGGTTTCTCCGTTTGTCGGGCGTATTTAttcct                                          >  1:1860740/1‑87 (MQ=255)
                                                                  cTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAgg                                 <  1:3109681/100‑1 (MQ=255)
                                                                   tggcGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATc                                    <  1:3693960/93‑1 (MQ=255)
                                                                   tGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGc                                >  1:662410/1‑100 (MQ=255)
                                                                   tGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGc                                >  2:2225714/1‑100 (MQ=255)
                                                                        ctgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGAGTGGTATCAGGCACGCa                           >  1:692172/1‑100 (MQ=255)
                                                                        ctgctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGGTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCa                           <  2:3537635/100‑1 (MQ=255)
                                                                          gctgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGCCTGGTATCAGGCACGCa                           >  1:1941140/1‑98 (MQ=255)
                                                                            tgTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGcc                       <  2:286173/100‑1 (MQ=255)
                                                                             gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGc                        >  1:3679757/1‑98 (MQ=255)
                                                                             gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCg                      >  2:97348/1‑100 (MQ=255)
                                                                             gTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCg                      <  1:713846/100‑1 (MQ=255)
                                                                               tttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGGTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGAt                    >  2:1742253/1‑100 (MQ=255)
                                                                                ttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATg                   <  1:433363/100‑1 (MQ=255)
                                                                                ttCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGAATTTATGCCTGGTATCAGGCCCGCAAGCCGagg                   >  2:4259186/1‑98 (MQ=255)
                                                                                 tCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATgg                  <  2:2193657/100‑1 (MQ=255)
                                                                                 gcttccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATgg                  >  2:3040065/7‑100 (MQ=39)
                                                                                  cTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGa                 >  1:2453167/1‑100 (MQ=255)
                                                                                   ttttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAc                >  2:2311900/1‑100 (MQ=255)
                                                                                   ttccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAc                >  2:1377504/5‑100 (MQ=255)
                                                                                    tttGCAAAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTAGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAcc               >  1:1760368/1‑100 (MQ=255)
                                                                                    tccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAcc               >  2:4013526/4‑100 (MQ=255)
                                                                                     ttGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc              <  2:160754/100‑1 (MQ=255)
                                                                                     ccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc              >  1:2922594/3‑100 (MQ=255)
                                                                                     ccgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAccc              >  2:3880645/3‑100 (MQ=255)
                                                                                      cgCACAGGCCCGTGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGAcc               >  2:3724241/2‑98 (MQ=255)
                                                                                       gCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGt            <  1:1017609/100‑1 (MQ=255)
                                                                                       gCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGt            >  1:3579590/1‑100 (MQ=255)
                                                                                       gCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGt            <  2:3842908/100‑1 (MQ=255)
                                                                                        cacaGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa           <  2:614105/100‑1 (MQ=255)
                                                                                        cacaGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTa           >  1:227591/1‑100 (MQ=255)
                                                                                         acaGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAt          <  2:3063198/100‑1 (MQ=255)
                                                                                          caGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtg         >  2:794296/1‑100 (MQ=255)
                                                                                           aGGCACGGGCCTGTGCGAAAGCAGGCGTTTTTCTGATTTCGCCGTATGTCGGGCGTATTCATGACTGGTATCAGTCACGCCAGCCGCTGGACCCg             <  1:4283971/95‑1 (MQ=255)
                                                                                            ggCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTCTGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgtg       >  2:1563236/1‑100 (MQ=255)
                                                                                            ggCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgtg       >  2:366602/1‑100 (MQ=255)
                                                                                            ggCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgtg       <  2:1942628/100‑1 (MQ=255)
                                                                                            ggCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTAtgtg       >  1:2547238/1‑100 (MQ=255)
                                                                                              cACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggt     >  1:1950133/1‑100 (MQ=255)
                                                                                               aCGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATgtggtg    <  1:786942/100‑1 (MQ=255)
                                                                                                cGGGCCTGTGCGGAAGCAGGCGTTTTTCTTATTCCGCCGTTTGTCGGGCGTCTTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg   <  2:2016213/100‑1 (MQ=255)
                                                                                                cGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg   <  2:2833914/100‑1 (MQ=255)
                                                                                                cGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGtggtgg   <  2:227592/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACAGGCACGGGCCTGTGCGGAAGCAGGCGTTTTTCTGATTTCGCCGTTTGTCGGGCGTATTTATGACTGGTATCAGGCACGCAAGCCGATGGACCCGTATGTGGTGGA  >  NC_000913/2579059‑2579255

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