Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,566 G→T 16.4% intergenic (+24/+17) citB → / ← dcuC response regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913654,5660GT16.4% 352.7 / 53.4 133intergenic (+24/+17)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base G (53/58);  new base T (10/12);  total (63/70)
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.08e-01

Only 100 of 136 total aligned reads displayed.

CAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCA  >  NC_000913/654467‑654662
                                                                                                   |                                                                                                
cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                                  <  2:105215/100‑1 (MQ=255)
cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                                  <  2:645136/100‑1 (MQ=255)
  gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCCCAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTcc                                                                                                <  2:1301986/100‑1 (MQ=255)
  gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                                <  2:249012/100‑1 (MQ=255)
  gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                                <  2:280874/100‑1 (MQ=255)
      ccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCg                                                                                            >  2:2508514/1‑100 (MQ=255)
        aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  2:1758736/100‑1 (MQ=255)
        aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  2:75840/100‑1 (MQ=255)
        aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  1:2775493/100‑1 (MQ=255)
        aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  1:2246497/100‑1 (MQ=255)
        aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                          <  2:1686760/100‑1 (MQ=255)
          cTGATTATTGCTGAAATTGTTCACGGCAAAGTGGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                         <  2:869348/99‑2 (MQ=38)
           tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCttt                                                                                       <  2:246598/100‑1 (MQ=38)
           tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCttt                                                                                       <  2:246317/100‑1 (MQ=38)
           tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                       >  2:447694/1‑100 (MQ=255)
            gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                      >  2:2397374/1‑100 (MQ=255)
            gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                      <  1:2146407/100‑1 (MQ=255)
             attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGCTTTat                                                                                     <  2:1660722/100‑2 (MQ=37)
             attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                     <  2:648104/100‑1 (MQ=255)
             attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                     >  1:648101/1‑100 (MQ=255)
                attGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGCTTTatcat                                                                                  <  1:1265840/100‑5 (MQ=37)
                 ttGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                 >  2:551389/1‑100 (MQ=255)
                  tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                <  2:356057/100‑1 (MQ=255)
                  tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                <  2:19021/100‑1 (MQ=255)
                  tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                <  1:1227949/100‑1 (MQ=255)
                  tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                >  1:651208/1‑100 (MQ=255)
                   gCTGAAATTGTTCACGGCAAAGTTGGCAGACCAGAACGCATATACCACAGTGGGTGACATGATAAAGCCGAAGGAAACTTGCCTCCGGCATTACTATTTa                                                                               >  1:1753855/1‑100 (MQ=255)
                     tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                             <  2:257740/100‑1 (MQ=255)
                      gAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACtt                                                                            >  2:1744174/1‑100 (MQ=255)
                       aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                           <  1:2166605/100‑1 (MQ=255)
                           tGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtg                                                                       <  2:2247862/100‑1 (MQ=255)
                            gTTCACGGGAAAGTTGGCAGACCACAGCGGATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtg                                                                       <  1:20850/99‑1 (MQ=255)
                             ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgtg                                                                     <  1:2631825/100‑1 (MQ=255)
                               cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                   <  1:551386/100‑1 (MQ=255)
                                aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAcc                                                                  >  2:2123080/1‑100 (MQ=255)
                                 cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCg                                                                 >  2:1571886/1‑100 (MQ=255)
                                  ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                                <  2:651211/100‑1 (MQ=255)
                                  ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                                >  1:618010/1‑100 (MQ=255)
                                  ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                                >  2:1276400/1‑100 (MQ=255)
                                   gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgct                                                               >  2:2032780/1‑100 (MQ=255)
                                      aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                            >  1:1661182/1‑100 (MQ=255)
                                      aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                            >  2:1168483/1‑100 (MQ=255)
                                         ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                         <  1:869409/100‑1 (MQ=255)
                                         ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                         <  1:2742408/100‑1 (MQ=255)
                                           ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGt                                                       >  1:2754244/1‑100 (MQ=255)
                                            gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                      <  1:1360432/100‑1 (MQ=255)
                                            gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                      <  2:1585806/100‑1 (MQ=255)
                                                aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                  <  2:2586434/100‑1 (MQ=255)
                                                aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                  <  1:2180268/100‑1 (MQ=255)
                                                 ccACAACACACATACCACAGTGGGTGACATGAGAAAACCCGAGGAAACTTGCCCCCGGCATTACCATTTACTTGCCTGTGACCGCTGCTGCCGTACCTgg                                                 <  1:168424/100‑1 (MQ=255)
                                                  cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTc                                                     <  2:1999666/95‑1 (MQ=255)
                                                  cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                                <  1:1042649/100‑1 (MQ=255)
                                                  cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                                <  1:2675010/100‑1 (MQ=255)
                                                  cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                                <  2:1661189/100‑1 (MQ=255)
                                                  cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                                <  1:1597446/100‑1 (MQ=255)
                                                      aCGCATATACCACATTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                             <  1:2901944/99‑1 (MQ=255)
                                                      aCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                            >  1:2442518/1‑100 (MQ=255)
                                                       cGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTGCGGCATTACTATTTACTTGCCTGTAACCGCTGCTGCCGTTCCTGGCACCAt                                           >  1:385313/1‑100 (MQ=255)
                                                       cGCATATACCACAGTGGGGGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCAt                                           >  1:2065067/1‑100 (MQ=255)
                                                        gCATATACCACAGTGGGTGACATGATAAAGCCGGCGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATc                                          <  1:1536268/100‑1 (MQ=255)
                                                         natatACCACAGTGGGTGCCATGATAAAGCCGCAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                            >  1:1914487/2‑97 (MQ=255)
                                                          atatACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAg                                        <  1:395316/100‑1 (MQ=255)
                                                            ataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGct                                      >  2:1885985/1‑100 (MQ=255)
                                                             taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                     >  2:1268963/1‑100 (MQ=255)
                                                              aCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                    <  1:2820503/100‑1 (MQ=255)
                                                               ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                   >  2:515214/1‑100 (MQ=255)
                                                                 acaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTa                                 <  2:1569518/100‑1 (MQ=255)
                                                                   aGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                               >  1:738046/1‑100 (MQ=255)
                                                                       ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTTCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACg                           <  1:1102637/100‑1 (MQ=255)
                                                                            cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATCAc                      >  1:959450/1‑100 (MQ=255)
                                                                                gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  >  1:2225713/2‑100 (MQ=38)
                                                                                aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAGCATCAGGTCTGTAGCAACGATAACAAt                   >  1:2430773/1‑99 (MQ=255)
                                                                                aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  <  2:210005/100‑1 (MQ=255)
                                                                                   atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc               >  2:2202814/3‑100 (MQ=255)
                                                                                   aaGCCGGAGGAAACTTGCCTCCGGCATTACTTTTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc               >  1:2259119/1‑100 (MQ=255)
                                                                                   aaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc               <  2:2260560/100‑1 (MQ=255)
                                                                                   aaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAACATCAGCTCTGTAGCAACGATAACAATCacc               >  2:1417500/1‑100 (MQ=255)
                                                                                    tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  1:204269/2‑100 (MQ=255)
                                                                                    tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  1:701096/2‑100 (MQ=255)
                                                                                    tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  2:1417502/2‑100 (MQ=255)
                                                                                    tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACATGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc               >  2:2055652/2‑99 (MQ=38)
                                                                                    aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  2:2191270/1‑100 (MQ=255)
                                                                                    aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  1:528204/1‑100 (MQ=255)
                                                                                    aGCCGGAGGAAACTTGCCTCCGGCATTACCATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  2:332881/1‑100 (MQ=255)
                                                                                     gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:1173643/1‑100 (MQ=255)
                                                                                     gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:1405022/1‑100 (MQ=255)
                                                                                     gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:262547/1‑100 (MQ=255)
                                                                                       cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac           <  1:361477/100‑1 (MQ=255)
                                                                                        ggaggcAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          >  2:539881/1‑100 (MQ=255)
                                                                                        ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          <  1:1097321/100‑1 (MQ=255)
                                                                                        ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          >  2:1078370/1‑100 (MQ=255)
                                                                                        gaaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          >  2:1976823/3‑100 (MQ=255)
                                                                                           ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac       >  2:2818307/1‑100 (MQ=255)
                                                                                           ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac       <  1:269016/100‑1 (MQ=255)
                                                                                           ggaAACTTGCCTCAGGCATTAATATTTACAGGCCTGTGACCGCTGCTGCCGTTCCTGGGACCATCAGCGCAGTAGAAACGATAACAATCAACAAACcaac       >  1:1799709/1‑100 (MQ=255)
                                                                                            gaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaaca      <  2:2520321/100‑1 (MQ=255)
                                                                                            gaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac       <  1:2828900/99‑1 (MQ=255)
                                                                                             aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaacaa     >  1:2023006/1‑100 (MQ=255)
                                                                                               aCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGc   <  2:2241976/100‑1 (MQ=255)
                                                                                                cTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCa  >  2:1891754/1‑100 (MQ=255)
                                                                                                   |                                                                                                
CAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCA  >  NC_000913/654467‑654662

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: