Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,763 A→T 15.6% intergenic (+29/+19) yeaD → / ← yeaE D‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,864,7630AT15.6% 249.1 / 35.3 90intergenic (+29/+19)yeaD/yeaED‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base T (7/7);  total (45/45)
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.70e-01

CAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864664‑1864862
                                                                                                   |                                                                                                   
cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                     >  1:1329197/1‑100 (MQ=255)
    aaCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCa                                                                                                 >  2:2023525/1‑100 (MQ=255)
      cGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATc                                                                                               >  1:2706042/1‑100 (MQ=255)
       gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              >  1:1295199/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCGTCCg                                                                                             <  2:1293200/100‑1 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              >  1:2297972/1‑99 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             <  2:3454664/100‑1 (MQ=255)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                            <  1:732590/100‑1 (MQ=39)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:2647553/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:1355565/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:1713126/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:2627038/100‑1 (MQ=38)
           aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                          <  2:1684541/100‑1 (MQ=255)
            aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  1:1934135/100‑1 (MQ=38)
             aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcaca                                                                                        >  2:397423/1‑100 (MQ=255)
             aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcaca                                                                                        <  2:2573591/100‑1 (MQ=255)
             aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcaca                                                                                        <  2:1965535/100‑1 (MQ=255)
              gTGACCAAAGAGAAACCTGTACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                       <  2:2614964/100‑1 (MQ=255)
               tGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacaca                                                                                      <  1:1879207/100‑1 (MQ=255)
                 aCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaac                                                                                    >  2:122599/1‑100 (MQ=255)
                  ncAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  1:2859476/2‑100 (MQ=255)
                     aaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATca                                                                                >  2:2681957/1‑100 (MQ=255)
                      agagAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcac                                                                               <  1:1477396/100‑1 (MQ=255)
                      agagAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcac                                                                               <  2:1691898/100‑1 (MQ=255)
                         naaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            >  2:2760739/2‑100 (MQ=255)
                                 cacaTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAg                                                                    >  2:2813562/1‑100 (MQ=255)
                                   caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  >  1:2040515/1‑100 (MQ=255)
                                   caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  >  2:3455665/1‑100 (MQ=255)
                                    aTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcgn                                                                 <  2:2859499/100‑2 (MQ=255)
                                      cTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAg                                                               >  2:1661909/1‑100 (MQ=255)
                                       tGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGt                                                              <  2:2467561/100‑1 (MQ=255)
                                         gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttt                                                            >  1:531482/1‑100 (MQ=255)
                                         gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttt                                                            <  1:344279/100‑1 (MQ=255)
                                         gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttt                                                            <  1:2800572/100‑1 (MQ=255)
                                          cgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGtttt                                                           >  2:1088486/1‑100 (MQ=255)
                                            cAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTc                                                         >  1:503252/1‑100 (MQ=255)
                                             aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        >  1:1632323/1‑100 (MQ=255)
                                             aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        <  2:443969/100‑1 (MQ=255)
                                              aTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCt                                                       <  1:2743987/100‑1 (MQ=255)
                                                ccATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCttt                                                     >  1:1812262/1‑100 (MQ=255)
                                                 cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtttt                                                    <  2:1593695/100‑1 (MQ=255)
                                                  aTTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                   >  1:1619361/1‑100 (MQ=255)
                                                  aTTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                   <  1:151954/100‑1 (MQ=255)
                                                   ttCGCGTTGCGAAACGTTAATTTACGTTCATGTTGTGTGCCGGGTCCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtt                                                      <  1:2476734/96‑1 (MQ=255)
                                                   ttCGCGTCGCGAAACGTTAATTTACTTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTATTCCTTTTgg                                                  <  1:2928258/100‑1 (MQ=255)
                                                    tCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGt                                                 <  1:2649609/100‑1 (MQ=255)
                                                     cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                <  2:1024846/100‑1 (MQ=255)
                                                     cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                <  1:997412/100‑1 (MQ=255)
                                                     cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                >  1:737051/1‑100 (MQ=255)
                                                     cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                >  2:360179/1‑100 (MQ=255)
                                                          tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                           <  2:670274/100‑1 (MQ=255)
                                                          cgCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                           <  2:403728/99‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  1:903153/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:2668323/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:3216007/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:1695408/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:2040532/100‑1 (MQ=255)
                                                                  gTTAATTTACGTTAATGTTGTGTGACGGGTGCAATGCATCCGTCACACAACATCACACAATATCCAGCGCAGTTTTTCCGTCTGGTGCCGGATATGCCtt                                   >  1:3336989/1‑100 (MQ=255)
                                                                   ttAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGTCGGATATGCCtt                                   >  1:2136703/1‑99 (MQ=255)
                                                                           cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                          >  2:960831/1‑100 (MQ=255)
                                                                           cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                          >  2:3052205/1‑100 (MQ=255)
                                                                           cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGAAGTTTTTCCTTTCGGTGCCGGATATGCCTTATCAAGCAt                          >  1:838863/1‑100 (MQ=255)
                                                                              tAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGc                       >  2:1262707/1‑100 (MQ=255)
                                                                               aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                      >  1:1795137/1‑100 (MQ=255)
                                                                               aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                      <  1:1963030/100‑1 (MQ=255)
                                                                                 tgttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTaa                    >  1:1782795/1‑100 (MQ=255)
                                                                                  gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  2:3328245/1‑100 (MQ=255)
                                                                                    tgtgtgCCGGGTGCAATGCATCCGGCACACATCATCATACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                 <  1:2336008/100‑1 (MQ=255)
                                                                                    tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                 >  2:840009/1‑100 (MQ=255)
                                                                                     gtgtgCCGGATGCATTGCACCCGGCACACAACATTACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTcc                >  1:374784/1‑100 (MQ=37)
                                                                                      tgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCg               >  2:566833/1‑100 (MQ=255)
                                                                                       gtgCCGGATGCATTGCACCCGGCACACAACATTACGCCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  2:1900982/1‑100 (MQ=37)
                                                                                       atgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  2:3305404/99‑1 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1063249/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:940624/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            <  1:1701686/100‑1 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:2581236/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            <  1:323665/100‑1 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:2683446/1‑100 (MQ=255)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:1363859/1‑100 (MQ=38)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:874267/1‑100 (MQ=38)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1369807/1‑100 (MQ=38)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1073685/1‑100 (MQ=38)
                                                                                            gggTGGAATGCATCCGGCACACCACATCACCCCATATCGAGCGCAGTGTCTACTTTTAGTGCCGGATTTGCCTTCTCCAGCATACCTTATTCCTCTtaa          >  1:1264370/1‑99 (MQ=255)
                                                                                            gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:874319/100‑1 (MQ=255)
                                                                                            gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:51870/100‑1 (MQ=255)
                                                                                            ggatGCATTGCACCCGGCACACAACATTACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:374787/97‑1 (MQ=38)
                                                                                              gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagaa       >  2:2182828/1‑100 (MQ=255)
                                                                                               tGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGaaa      >  1:166943/1‑100 (MQ=255)
                                                                                                gCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAg     <  2:3171066/100‑1 (MQ=255)
                                                                                                 cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGt    <  1:3344208/100‑1 (MQ=255)
                                                                                                  aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt   >  1:3239496/1‑100 (MQ=255)
                                                                                                   aTGCATCCGGCACACCACATCACACCATATCCAGCGCAGATTTTCCTTTTGGTGCCGGATATGCCTTATCCCGCATACCTAATTCCGCTGAAGAAAGtt   >  2:500814/1‑99 (MQ=255)
                                                                                                   aTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTTTTTTTTTTTCCGGATATGCCTTATCCAGCATAGCTACTTCCGCTGAAGAAAGTTc  >  1:2580595/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864664‑1864862

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