Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,560 A→C 16.7% intergenic (+18/+23) 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,5600AC16.7% 284.3 / 54.4 126intergenic (+18/+23)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base A (48/56);  new base C (9/12);  total (57/69)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 140 total aligned reads displayed.

TTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654461‑654659
                                                                                                   |                                                                                                   
ttGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCggagga                                                                                                     <  1:2721551/100‑1 (MQ=255)
 tGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaa                                                                                                    >  2:209685/1‑100 (MQ=255)
  gcgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaaa                                                                                                   >  2:2726275/1‑100 (MQ=255)
  gcgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaaa                                                                                                   <  1:2558632/100‑1 (MQ=255)
   cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                  <  2:1371685/100‑1 (MQ=255)
      cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                               <  2:645136/100‑1 (MQ=255)
      cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                               <  1:203593/100‑1 (MQ=255)
      cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                               <  2:105215/100‑1 (MQ=255)
        gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                             <  2:249012/100‑1 (MQ=255)
            ccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCg                                                                                         >  2:2508514/1‑100 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCgg                                                                                        <  1:2902230/99‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  2:1758736/100‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:2775493/100‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  2:75840/100‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:2246497/100‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       <  2:1929586/100‑1 (MQ=255)
              aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       <  2:1686760/100‑1 (MQ=255)
               tCTGATTATTGCTGAACTTGTTCACGGCAACGTTGGCAGACCAAACCGCATATCCCACAGTGGGTGACATGCTAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                      <  2:131683/100‑2 (MQ=255)
                cTGATTATTGCTGAAATTGTTCACGGCAAAGTGGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                      <  2:869348/99‑2 (MQ=38)
                 tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCttt                                                                                    <  2:246317/100‑1 (MQ=38)
                 tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                    >  2:447694/1‑100 (MQ=255)
                  gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                   <  1:2146407/100‑1 (MQ=255)
                  gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                   >  2:2397374/1‑100 (MQ=255)
                   attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGCTTTat                                                                                  <  2:1660722/100‑2 (MQ=37)
                   attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                  >  1:648101/1‑100 (MQ=255)
                   attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                  <  2:648104/100‑1 (MQ=255)
                     tattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                <  1:2214516/100‑1 (MQ=255)
                      attGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGCTTTatcat                                                                               <  1:1265840/100‑5 (MQ=37)
                       ttGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                              >  2:551389/1‑100 (MQ=255)
                        tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             <  1:1227949/100‑1 (MQ=255)
                        tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             <  2:356057/100‑1 (MQ=255)
                        tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             >  1:651208/1‑100 (MQ=255)
                        tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             <  2:19021/100‑1 (MQ=255)
                          cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGTCATTACTATTTAc                                                                           >  2:2241177/1‑100 (MQ=255)
                           tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                          <  2:257740/100‑1 (MQ=255)
                            gAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACtt                                                                         >  2:1744174/1‑100 (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)
                                     cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                >  2:1492749/1‑100 (MQ=255)
                                       cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCg                                                              >  2:1571886/1‑100 (MQ=255)
                                        ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                             >  1:618010/1‑100 (MQ=255)
                                        ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                             >  2:1276400/1‑100 (MQ=255)
                                        ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                             <  2:651211/100‑1 (MQ=255)
                                         gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgct                                                            >  2:2032780/1‑100 (MQ=255)
                                            aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                         >  2:1168483/1‑100 (MQ=255)
                                              gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCtgctgc                                                       <  1:1669604/100‑1 (MQ=255)
                                               ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                      <  2:48233/100‑1 (MQ=255)
                                               ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                      <  1:2742408/100‑1 (MQ=255)
                                                tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                     >  1:1511417/1‑100 (MQ=255)
                                                  gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                   <  2:1585806/100‑1 (MQ=255)
                                                  gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                   >  2:2472254/1‑100 (MQ=255)
                                                    aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGGTcc                                                 >  2:2557250/1‑100 (MQ=255)
                                                     gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                <  2:1248107/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:2675010/100‑1 (MQ=255)
                                                        cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                             <  1:1042649/100‑1 (MQ=255)
                                                        cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                             <  1:1597446/100‑1 (MQ=255)
                                                        cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                             <  2:1661189/100‑1 (MQ=255)
                                                           aaCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                          >  2:2891709/1‑100 (MQ=255)
                                                            aCGCATATACCACATTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                          <  1:2901944/99‑1 (MQ=255)
                                                            aCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                         >  1:2442518/1‑100 (MQ=255)
                                                             cGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTGCGGCATTACTATTTACTTGCCTGTAACCGCTGCTGCCGTTCCTGGCACCAt                                        >  1:385313/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                                  >  1:2100061/1‑100 (MQ=255)
                                                                   taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                  >  2:1268963/1‑100 (MQ=255)
                                                                    aCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                 <  1:2820503/100‑1 (MQ=255)
                                                                         aGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                            >  1:738046/1‑100 (MQ=255)
                                                                          gTGGGTGACATGATAAAGCCGGAGGGAACTTGCCGCCGGCATTACTATTTACTGGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAg                             <  1:885045/98‑1 (MQ=255)
                                                                           tGGGTGACATGATAAAGCCGGAGGAAACGTGCCTGCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCaa                          <  1:2819606/100‑1 (MQ=255)
                                                                             ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTTCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACg                        <  1:1102637/100‑1 (MQ=255)
                                                                                  cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAGTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                   <  1:1858630/100‑1 (MQ=255)
                                                                                  cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATCAc                   >  1:959450/1‑100 (MQ=255)
                                                                                      aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc               <  2:210005/100‑1 (MQ=255)
                                                                                       tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa              >  1:2333974/1‑100 (MQ=255)
                                                                                        aaaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  2:2138686/1‑100 (MQ=255)
                                                                                         atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa              >  2:1079106/3‑98 (MQ=39)
                                                                                         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           >  2:1417502/2‑100 (MQ=255)
                                                                                          tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca           >  1:701096/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)
                                                                                           gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  2:262547/1‑100 (MQ=255)
                                                                                            ccGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa         >  2:2787269/1‑100 (MQ=38)
                                                                                             cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac        <  1:361477/100‑1 (MQ=255)
                                                                                              ggaggcAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc       >  2:539881/1‑100 (MQ=255)
                                                                                              ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc       <  1:1097321/100‑1 (MQ=255)
                                                                                                 ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac    >  2:2818307/1‑100 (MQ=255)
                                                                                                 ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac    <  1:269016/100‑1 (MQ=255)
                                                                                                 ggaAACTTGCCTCAGGCATTAATATTTACAGGCCTGTGACCGCTGCTGCCGTTCCTGGGACCATCAGCGCAGTAGAAACGATAACAATCAACAAACcaac    >  1:1799709/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654461‑654659

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: