Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,560 A→C 16.6% 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.6% 439.9 / 86.1 187intergenic (+18/+23)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base A (77/79);  new base C (15/16);  total (92/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 201 total aligned reads displayed.

TGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654462‑654659
                                                                                                  |                                                                                                   
tGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaa                                                                                                    <  2:1272977/100‑1 (MQ=255)
  cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                  <  2:4147063/100‑1 (MQ=255)
    ccAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGGATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                  >  1:1652792/1‑98 (MQ=255)
    ccAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                >  2:1858480/1‑100 (MQ=255)
          gccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTtcctcc                                                                                          <  1:1780480/100‑1 (MQ=255)
             aTCTGCTTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:727993/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTTACGGTAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:1699587/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:1715313/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:792516/100‑1 (MQ=255)
              tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                      <  1:2639623/100‑2 (MQ=255)
              tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       <  2:243509/99‑1 (MQ=255)
              tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                      <  1:4260110/100‑1 (MQ=255)
              tCTGAGTATTGCTGAAAGTGTTGACGGCAAAGTTGGCAGAGCACAACGCATATACCACGGGGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:742393/99‑1 (MQ=255)
               cTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGctt                                                                                     <  2:112760/100‑3 (MQ=39)
                tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGctt                                                                                     <  2:3981806/99‑3 (MQ=39)
                 gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                   >  2:910465/1‑100 (MQ=255)
                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                  <  2:242158/100‑1 (MQ=255)
                    tattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                >  1:2541015/1‑100 (MQ=255)
                    tattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                >  1:331951/1‑100 (MQ=255)
                       tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             >  1:2238136/1‑100 (MQ=255)
                       tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             >  2:4152696/1‑100 (MQ=255)
                       tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                             <  1:2829555/100‑1 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           <  2:1527851/100‑1 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           >  2:3288295/1‑100 (MQ=255)
                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                          >  1:1019286/1‑100 (MQ=255)
                            aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                        <  1:2989141/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGCCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      <  1:3509612/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGTTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      <  1:691902/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      <  2:2238146/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      <  2:1642975/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      <  2:2031259/100‑1 (MQ=255)
                              aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                      >  1:3912783/1‑100 (MQ=255)
                                tGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtg                                                                    <  1:3237088/100‑1 (MQ=255)
                                 gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                   >  2:373726/1‑100 (MQ=255)
                                 gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                   <  2:28674/100‑1 (MQ=255)
                                   tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                 <  2:3088750/100‑1 (MQ=255)
                                   tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                 >  1:1904793/1‑100 (MQ=255)
                                       ggCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                             <  1:563167/100‑1 (MQ=255)
                                         cAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctg                                                           >  1:1122656/1‑100 (MQ=255)
                                            aGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGctgctg                                                        <  1:3564447/100‑1 (MQ=255)
                                             gTTGGCAGACCACCACGCATCTACCCCAGTGGGTGACATGTTCAAGCCGGAGGAATCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgctc                                                        <  1:2752312/99‑2 (MQ=255)
                                             gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCtgctgc                                                       <  2:781396/100‑1 (MQ=255)
                                               tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                     <  1:2888566/100‑1 (MQ=255)
                                               tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                     >  2:465638/1‑100 (MQ=255)
                                               tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                     <  2:3229536/100‑1 (MQ=255)
                                                 gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                   >  2:2704808/1‑100 (MQ=255)
                                                 gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                   >  2:597961/1‑100 (MQ=255)
                                                 gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTCCCGtt                                                   <  1:2935678/100‑1 (MQ=255)
                                                 gCAGACCACAACGCATAGACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                   <  1:4339455/100‑1 (MQ=255)
                                                  cAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTc                                                  <  1:3125806/100‑1 (MQ=255)
                                                   aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                 <  2:1789131/100‑1 (MQ=255)
                                                   aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                 <  2:1198643/100‑1 (MQ=255)
                                                    gACCACAACGCATATACGACAGTGAATGACATCATAAAGCCGGGGGAAACTCGCCTCCGGCCTTACTGATTCCTTGCCTGTGCACgctgc                                                          >  1:2623991/1‑90 (MQ=255)
                                                     aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                               <  1:1858473/100‑1 (MQ=255)
                                                      ccACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                              >  2:2227864/1‑100 (MQ=255)
                                                      ccACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                              >  1:3407857/1‑100 (MQ=255)
                                                      agacaacaCATATACCACAGTGGGCGACACGCTAACCCGGGATTAAACTTGCTTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                              >  1:990960/3‑100 (MQ=255)
                                                       cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTAACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                             <  1:3985372/100‑1 (MQ=255)
                                                        acaACGCATATACCACAGTGGTTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                            <  1:2175105/100‑1 (MQ=255)
                                                        acaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                            <  1:247739/100‑1 (MQ=255)
                                                          aaCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                          <  2:2281869/100‑1 (MQ=255)
                                                           aCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                         >  2:3650376/1‑100 (MQ=255)
                                                            cGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCAt                                        >  2:2261900/1‑100 (MQ=255)
                                                             gCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATc                                       <  1:2491881/100‑1 (MQ=255)
                                                                tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGc                                    <  2:2112564/100‑1 (MQ=255)
                                                                tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGc                                    >  1:3411801/1‑100 (MQ=255)
                                                                  taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                  <  1:316435/100‑1 (MQ=255)
                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                               >  2:3274321/1‑100 (MQ=255)
                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                               >  2:3690466/1‑100 (MQ=255)
                                                                      acaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                               <  1:1354165/99‑1 (MQ=255)
                                                                       caGTGGGTGACATGATAAAGCCGGAGGAGACTTGCCGCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAg                             <  1:859492/100‑1 (MQ=255)
                                                                        aGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                            <  1:2913395/100‑1 (MQ=255)
                                                                         gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACAACTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCa                           >  2:556783/1‑100 (MQ=255)
                                                                           gggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                         >  2:1455649/1‑100 (MQ=255)
                                                                                aaatagtAATGCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGTCACCATCCGCTCTGTAGCAACGATaa                    >  1:4155250/7‑100 (MQ=21)
                                                                                 cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                   <  2:2512855/100‑1 (MQ=255)
                                                                                 cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                   >  1:2685270/1‑100 (MQ=255)
                                                                                   tGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACaa                 >  1:4011989/1‑100 (MQ=255)
                                                                                     gtAATGCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc               >  2:3596949/2‑100 (MQ=37)
                                                                                     gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc               >  2:342800/2‑100 (MQ=38)
                                                                                     aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc               <  1:204067/100‑1 (MQ=255)
                                                                                      tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa              >  1:2002383/1‑100 (MQ=255)
                                                                                       aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCCCCATCACCTCTGTAGCAACGATAACAATCac             >  1:2061184/4‑100 (MQ=38)
                                                                                       aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  1:1944511/4‑100 (MQ=38)
                                                                                        atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTACCGTTCCTGGCACCATCAGCTCTGTAGCATCGATAACACTCacc            >  1:3599819/3‑100 (MQ=255)
                                                                                         tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACAGCTGATGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  1:3730109/2‑98 (MQ=38)
                                                                                         aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             <  1:4177233/98‑1 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  2:2861138/1‑100 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:4083445/1‑100 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:2383257/1‑100 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          <  2:2880841/100‑1 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  2:3807383/1‑100 (MQ=255)
                                                                                            cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac        <  1:1455643/100‑1 (MQ=255)
                                                                                            cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac        >  2:43951/1‑100 (MQ=255)
                                                                                            cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac        <  2:3021966/100‑1 (MQ=255)
                                                                                             ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc       >  2:3827150/1‑100 (MQ=255)
                                                                                             ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc       <  2:2812522/100‑1 (MQ=255)
                                                                                             ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTCCTGCCGCTCCTTGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:1080525/1‑97 (MQ=255)
                                                                                                ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaaccaa     >  1:3800258/1‑99 (MQ=255)
                                                                                                  aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaacaa  <  2:1567041/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654462‑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: