Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,763 A→T 15.8% 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.8% 375.1 / 56.4 146intergenic (+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 (49/74);  new base T (7/16);  total (56/90)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 151 total aligned reads displayed.

CAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864664‑1864862
                                                                                                   |                                                                                                   
cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTCATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                     >  2:319751/1‑100 (MQ=255)
 aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATGCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGGGCAAt                                                                                                    >  1:2092520/1‑100 (MQ=255)
 aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCGGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                    <  2:3501932/100‑1 (MQ=255)
  gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATg                                                                                                   <  2:1801495/100‑1 (MQ=255)
    aaCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGc                                                                                                  >  1:1659673/1‑99 (MQ=255)
    aaCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCCGTTCCAATGCa                                                                                                 >  1:1659316/1‑100 (MQ=255)
       gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              <  2:2045816/100‑1 (MQ=255)
       gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              <  1:2856424/100‑1 (MQ=255)
       gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              <  1:4028950/100‑1 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             >  2:3593036/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             >  2:1639407/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGCGTGCAATGCATCCg                                                                                             >  1:491747/1‑100 (MQ=255)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATc                                                                                               <  1:3151300/97‑1 (MQ=255)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCgg                                                                                            <  2:274765/100‑1 (MQ=255)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                           <  1:96582/100‑1 (MQ=255)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:3629727/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:379526/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  1:270228/100‑1 (MQ=38)
           aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                          <  2:2503296/100‑1 (MQ=255)
           aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGTATTGCACCCGGca                                                                                          <  2:113086/100‑1 (MQ=38)
            aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  2:821670/100‑1 (MQ=38)
            aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  2:2385250/100‑1 (MQ=38)
            aaGTGACCAAAGAGAAACCTGCACAGCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  1:3161272/100‑1 (MQ=38)
            aaGTGACCAAAGAGAAACCCGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  2:1924291/100‑1 (MQ=38)
              gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                       <  2:176746/100‑1 (MQ=255)
              gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcacac                                                                                       <  2:1206184/100‑1 (MQ=38)
                 aCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaac                                                                                    <  2:3778331/100‑1 (MQ=255)
                 aCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaac                                                                                    <  2:3892891/100‑1 (MQ=255)
                  ccAAAGAGAAACCGGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  2:2869183/1‑100 (MQ=255)
                  ccAAAGAGAAAACTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  2:3395329/1‑100 (MQ=255)
                   cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                  >  1:1057312/1‑100 (MQ=255)
                   cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                  >  2:489215/1‑100 (MQ=255)
                       gagaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAttaac                                                                              <  2:3447175/100‑5 (MQ=25)
                         gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            <  1:4175476/100‑1 (MQ=255)
                          aaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCa                                                                           <  2:3322749/100‑1 (MQ=255)
                             ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCatat                                                                        >  2:1563951/1‑100 (MQ=255)
                               tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                      <  1:4245585/100‑1 (MQ=255)
                               tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                      >  2:295224/1‑100 (MQ=255)
                               tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                      <  2:2312519/100‑1 (MQ=255)
                               tGCACATCTGGCGAAATCCATTCGAGTTGCGAAACGTTACTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                      <  2:4276816/100‑1 (MQ=255)
                                   caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  <  1:2190356/100‑1 (MQ=255)
                                       tGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAg                                                               <  1:1298717/99‑1 (MQ=255)
                                          cgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGtttt                                                           >  2:2130918/1‑100 (MQ=255)
                                              aTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCt                                                       <  2:1959407/100‑1 (MQ=255)
                                               tCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGACGGGTGCAATGAATCCGGAACACAACATCACACCATATCAAGCGAAGTTTTTACtt                                                      >  1:3228482/1‑100 (MQ=255)
                                                  aTTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                   <  2:3516204/100‑1 (MQ=255)
                                                   ttCGCGTTGCGAAACGTTAATTTCCGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                  <  1:621748/100‑1 (MQ=255)
                                                   ttCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                  <  2:1881324/100‑1 (MQ=255)
                                                    tCGCGTTGGGAAACGTGAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                   <  1:4144385/98‑1 (MQ=255)
                                                      gcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGc                                               >  2:359740/1‑100 (MQ=255)
                                                       cgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGcc                                              <  1:1944459/100‑1 (MQ=255)
                                                        gTTTCGAAACGTTATGGTGCGTGAATGTTATGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATTGCCAGCGCAGTTTTTCCTTTTGGTGCCg                                             <  1:1246535/100‑1 (MQ=255)
                                                        gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                             <  2:4242701/100‑1 (MQ=255)
                                                         ttGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCgg                                            >  2:3102515/1‑100 (MQ=255)
                                                          tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                           >  1:2379815/1‑100 (MQ=255)
                                                           gCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                          >  2:3899999/1‑100 (MQ=255)
                                                            cGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGata                                         <  2:2452528/100‑1 (MQ=255)
                                                              aaaCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                       >  1:3020288/1‑100 (MQ=255)
                                                               aaCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGc                                      >  1:1574926/1‑100 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:2154703/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  1:2090227/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGcc                                     <  2:1354154/100‑1 (MQ=255)
                                                                aCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGc                                      <  2:1152853/99‑1 (MQ=255)
                                                                  gTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCAGCATCTGCCt                                    >  2:2925523/1‑99 (MQ=255)
                                                                      aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                               <  1:1933015/100‑1 (MQ=255)
                                                                      aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                               <  1:320622/100‑1 (MQ=255)
                                                                      aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATc                                >  2:2614498/1‑99 (MQ=255)
                                                                         tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                            <  1:8921/100‑1 (MQ=255)
                                                                         tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                            <  2:2754143/100‑1 (MQ=255)
                                                                             gtgaTGTTGTGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTTCCGGATATGCCTTATCCAGCATAg                        >  1:3717611/4‑100 (MQ=25)
                                                                              tgaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGc                       <  1:2588705/98‑1 (MQ=255)
                                                                              tAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGc                       <  1:3085929/100‑1 (MQ=255)
                                                                                  gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   <  2:630549/100‑1 (MQ=255)
                                                                                  gttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  1:83038/1‑100 (MQ=255)
                                                                                  gttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  1:1320892/1‑100 (MQ=37)
                                                                                  gttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  2:3738703/1‑100 (MQ=38)
                                                                                   ttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAtt                  >  1:4040977/1‑100 (MQ=38)
                                                                                     gtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTcc                <  2:2035967/100‑1 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  2:675655/1‑100 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  2:3805422/100‑1 (MQ=255)
                                                                                       gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTCATTCCGc              >  1:565722/1‑100 (MQ=38)
                                                                                        tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             <  2:310133/100‑1 (MQ=255)
                                                                                         gacgtgTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            <  1:2764625/95‑1 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:2776026/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGAGCAGTTTTTCCTTTTGGTGCCGGATCTGCCTTATCCAGCATAGCTAATTCCGCt             >  1:997107/1‑99 (MQ=255)
                                                                                          ccGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTga           <  2:649645/100‑1 (MQ=255)
                                                                                           caggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa          >  2:2886854/3‑100 (MQ=255)
                                                                                           cGGGTGCAATGTATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa          >  1:1310806/1‑100 (MQ=255)
                                                                                           cGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa          <  2:3244760/100‑1 (MQ=255)
                                                                                            gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         >  1:3836082/1‑100 (MQ=255)
                                                                                            gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:1280050/100‑1 (MQ=255)
                                                                                             ggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagc        >  1:2416819/1‑99 (MQ=255)
                                                                                             ggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaaga        >  2:1214416/1‑100 (MQ=255)
                                                                                              gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTGCCTTGGGGTGGCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagaa       <  1:1302446/100‑1 (MQ=255)
                                                                                                 cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGt    >  1:2820340/1‑100 (MQ=255)
                                                                                                 cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGt    >  2:1634272/1‑100 (MQ=255)
                                                                                                  aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt   >  2:1778616/1‑100 (MQ=255)
                                                                                                  aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt   <  1:1563947/100‑1 (MQ=255)
                                                                                                  aaTGCATCCGGCAAACAACATCACACCATATCCAGCGCAGTTTTCCCTTTTGGTGCCGGATATGTCTTATCTAGCATAGCTAATTCCGCTGAAGATAGt    >  1:2123960/1‑99 (MQ=255)
                                                                                                   aTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTc  >  2:2436395/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: