Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,778,364 C→T 18.6% intergenic (+17/+26) ydiO → / ← ydiP putative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,778,3640CT18.6% 258.4 / 44.5 113intergenic (+17/+26)ydiO/ydiPputative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (54/38);  new base T (13/8);  total (67/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 123 total aligned reads displayed.

TGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTC  >  NC_000913/1778266‑1778463
                                                                                                  |                                                                                                   
tgAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGct                                                                                                    >  2:840414/1‑100 (MQ=255)
tgAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGct                                                                                                    >  2:1486014/1‑100 (MQ=255)
   aCGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttct                                                                                                 >  2:4399598/1‑100 (MQ=255)
    cGTATGGGCGGCGGTACAGACGAAATGATGATTTACGTAGCAGGACGGCAGATCCTGGGAGACTCTCAGAACAAATAATCTGCAGGCGGCGCAGcttc                                                                                                  <  1:2865001/98‑1 (MQ=255)
    cGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttctt                                                                                                >  2:1943499/1‑100 (MQ=255)
      tATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                              <  2:1191421/100‑1 (MQ=255)
       aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                              >  1:4111789/1‑99 (MQ=255)
       aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                             <  2:547755/100‑1 (MQ=255)
       aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                             <  1:2234035/100‑1 (MQ=255)
       aTCGGCGGCGGTACAGACGAAATTATGATTTACATAGCAGGTCGGCAGAGCATGAAAGACTATCAGAACAGATAATCAGCAGGCGGCGCAGcttc                                                                                                  >  1:2647673/1‑95 (MQ=255)
        tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACa                                                                                            >  1:2262216/1‑100 (MQ=255)
          ggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACaaa                                                                                          >  2:3861571/1‑100 (MQ=255)
          ggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGCCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACaaa                                                                                          >  1:2397604/1‑100 (MQ=255)
            cggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACt                                                                                        <  2:85346/100‑1 (MQ=255)
                gggaCAGACGAAATCATGATTTACCTAGCAGGTCGTCAGATTCTGAAACACTATCAGAACAAATAAGCCCCAGACGACACAACTTCTTAACAAACTg                                                                                       >  1:3778970/4‑97 (MQ=255)
                ggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgcgc                                                                                    >  2:361866/1‑100 (MQ=255)
                    cAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcc                                                                                <  2:4040644/100‑1 (MQ=38)
                    cAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCgccgcc                                                                                >  2:4283899/1‑100 (MQ=255)
                    cAGACGAAATTATGATTTACGTAGCAGGTCAGCAGGTCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcc                                                                                <  1:3051266/100‑1 (MQ=37)
                      gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCa                                                                               >  1:3338136/1‑99 (MQ=255)
                      gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAg                                                                              >  2:3402503/1‑100 (MQ=255)
                      gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAg                                                                              >  2:3703704/1‑100 (MQ=255)
                       aCGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgc                                                                             <  1:2556532/100‑4 (MQ=37)
                        cGAAAGTGTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgc                                                                             <  1:2795227/99‑4 (MQ=25)
                         gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAt                                                                            >  1:260172/1‑99 (MQ=255)
                          aaaTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGCCTCTCAGACCAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgcag                                                                           <  1:443084/99‑6 (MQ=21)
                           aaTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgcaga                                                                          <  1:2618998/99‑7 (MQ=25)
                            aTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgcaga                                                                          <  2:4114616/98‑7 (MQ=25)
                            aTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTAt                                                                        >  1:4023335/1‑100 (MQ=255)
                             ttATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATc                                                                       <  1:1172009/100‑1 (MQ=255)
                              tATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATcc                                                                      <  2:2282022/100‑1 (MQ=255)
                              tATGATTTACGTAGCAGGTCGACAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATcc                                                                      <  2:1082856/100‑1 (MQ=255)
                               aTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCca                                                                     >  2:3246556/1‑100 (MQ=255)
                                tGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                    >  2:3924937/1‑100 (MQ=255)
                                 gATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                   >  2:3020103/1‑100 (MQ=255)
                                  aTTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaaca                                                                  <  1:4283868/100‑1 (MQ=255)
                                  aTTTACGCAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCGTAGCAAACTGCGCCGCCAGATTTATCcaaca                                                                  <  1:104002/100‑1 (MQ=255)
                                   tttACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaacaa                                                                 >  2:1928062/1‑100 (MQ=255)
                                    ttACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATACTCGGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAg                                                                <  1:1930240/100‑1 (MQ=255)
                                     tACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGa                                                               >  1:623147/1‑100 (MQ=255)
                                      aCGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGAc                                                              <  1:3026381/100‑1 (MQ=255)
                                       cGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACt                                                             >  1:3103422/1‑100 (MQ=255)
                                              ggTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCgg                                                       >  2:3846435/1‑99 (MQ=255)
                                              ggTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGt                                                      >  2:3837025/1‑100 (MQ=255)
                                              ggTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGt                                                      >  2:380147/1‑100 (MQ=255)
                                                 cGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGa                                                   >  2:577464/1‑100 (MQ=255)
                                                  ggCAGAGCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCCGATTTATCCAACAAGACTTACCGGTTGAg                                                  >  1:4058310/1‑100 (MQ=255)
                                                   gCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAgg                                                 <  2:3153929/100‑1 (MQ=255)
                                                    cAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGa                                                >  2:3243586/1‑100 (MQ=255)
                                                      gATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGCTTTACCCAACAAGACTTACCGGTCGAgcca                                               >  1:1000033/1‑96 (MQ=255)
                                                      gATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGaaa                                              >  1:2920511/1‑100 (MQ=255)
                                                      gATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAACTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGaaa                                              >  1:318908/1‑100 (MQ=255)
                                                         ccTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTc                                           <  2:2363933/100‑1 (MQ=255)
                                                          cTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCa                                          >  1:342315/1‑100 (MQ=255)
                                                           tGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAg                                         >  1:3353931/1‑100 (MQ=255)
                                                           tGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAg                                         >  2:3706384/1‑100 (MQ=255)
                                                            gAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGa                                        <  1:2606108/100‑1 (MQ=255)
                                                             aaaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGAt                                       <  1:1113209/100‑1 (MQ=255)
                                                              aaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATa                                      <  1:3124245/100‑1 (MQ=255)
                                                               aGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAg                                     >  1:2482362/1‑100 (MQ=255)
                                                                  cTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAt                                  >  2:751642/1‑100 (MQ=255)
                                                                   tATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAtt                                 <  2:2108767/100‑1 (MQ=255)
                                                                     tCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTc                               >  2:3182590/1‑100 (MQ=255)
                                                                      cAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCt                              <  2:1777233/100‑1 (MQ=255)
                                                                      cAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCt                              >  2:3854895/1‑100 (MQ=255)
                                                                       aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                             >  2:169261/1‑100 (MQ=255)
                                                                       aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCAAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                             <  1:2405676/100‑1 (MQ=255)
                                                                        gAACAAATAATCTGCAGGCGGTGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTTAGGAAATTCAGATAGGATTTCTTg                            <  1:2257125/100‑1 (MQ=255)
                                                                        gAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTg                            <  2:1239588/100‑1 (MQ=255)
                                                                         aaCAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGt                           <  2:2864765/100‑1 (MQ=255)
                                                                            aaaTAATCTGCAGGCGGCGCAGCTTCTTAACAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTcc                        <  1:4399567/100‑1 (MQ=255)
                                                                            aaaTAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTcc                        <  2:4023362/100‑1 (MQ=255)
                                                                            aaaTAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTcc                        <  1:90392/100‑1 (MQ=255)
                                                                             aataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCg                       >  1:453758/1‑100 (MQ=255)
                                                                               taatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCa                     <  1:616164/100‑1 (MQ=255)
                                                                                  tCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaaa                  >  1:3990801/1‑100 (MQ=255)
                                                                                   ccgatctcggcgCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                 <  2:1262789/93‑1 (MQ=255)
                                                                                   cTGCAGGCGGCGCAGCTTCTGAACACACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                 <  1:3936136/100‑1 (MQ=255)
                                                                                     gCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAg               <  1:474508/100‑1 (MQ=255)
                                                                                      ttggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCTGTTGAGGCAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt              >  1:644197/3‑100 (MQ=37)
                                                                                      cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt              >  1:3235184/1‑100 (MQ=255)
                                                                                      cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt              >  1:2464510/1‑100 (MQ=255)
                                                                                       tggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt             >  2:3745007/2‑100 (MQ=37)
                                                                                       tggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt             >  1:2813043/2‑100 (MQ=37)
                                                                                        ggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAACCAGtt             >  2:721026/1‑99 (MQ=37)
                                                                                        ggcggcGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttt            >  2:2802187/1‑100 (MQ=255)
                                                                                        ggccgcGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGCTTTCTTGTTCCGCAAAACAGttt            >  1:4159407/1‑100 (MQ=255)
                                                                                         gcggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt           >  2:2021232/1‑100 (MQ=38)
                                                                                          cggcgCAGTTTGTTAAGAAGGTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAGTCAGATACAAGTTCTTCTATAGCAAAACAGctt            >  1:4003813/1‑98 (MQ=16)
                                                                                          cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt          >  2:3056295/1‑100 (MQ=38)
                                                                                          cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt          >  1:334118/1‑100 (MQ=38)
                                                                                           ggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa         >  2:2279969/1‑100 (MQ=38)
                                                                                           ggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa         >  1:3985477/1‑100 (MQ=38)
                                                                                           ggcgCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa         <  1:2762435/100‑1 (MQ=255)
                                                                                            gcgcAGTTTGTTAAGAAGCTGCGCAGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTATTGTTACACAACACAGTTTTTaa        >  1:4185497/1‑100 (MQ=255)
                                                                                            gcgcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaa        >  2:2570068/1‑100 (MQ=255)
                                                                                             cgcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaaa       >  1:124749/1‑100 (MQ=255)
                                                                                               cAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAAtt     <  2:434071/100‑1 (MQ=255)
                                                                                                  cttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc  <  2:3221693/100‑1 (MQ=255)
                                                                                                  cttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc  <  2:3307764/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTC  >  NC_000913/1778266‑1778463

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: