Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,763 A→T 18.2% 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,7630AT18.2% 316.1 / 63.4 137intergenic (+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 (56/55);  new base T (11/14);  total (68/69)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 146 total aligned reads displayed.

CAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864664‑1864862
                                                                                                   |                                                                                                   
cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                     >  2:2189543/1‑100 (MQ=255)
cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                     >  2:226685/1‑100 (MQ=255)
cAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCaa                                                                                                     >  1:2007668/1‑100 (MQ=255)
     aCGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCAt                                                                                                <  2:1423310/100‑1 (MQ=255)
       gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                              >  2:1230336/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             >  1:1108871/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             <  2:2206206/100‑1 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             >  1:673233/1‑100 (MQ=255)
        cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCg                                                                                             <  1:4395723/100‑1 (MQ=39)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCgg                                                                                            <  2:3570147/100‑1 (MQ=255)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                            <  2:3979183/100‑1 (MQ=39)
         aaaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                            <  1:742523/100‑1 (MQ=39)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:2457322/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  1:2170596/100‑1 (MQ=38)
          aaaaGTGACCAAAGAGAAACCTGCACATCGGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:1324082/100‑1 (MQ=38)
           aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                          <  1:1369508/100‑1 (MQ=255)
           aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGca                                                                                          <  1:3645545/100‑1 (MQ=255)
            atgTGGCCTAAGAGAATCCTTCACTTCTGGCTCTATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                         <  1:2180562/98‑1 (MQ=255)
            aaGTGCCCAAAGCGAAACCTGGATATTTGGCTGAATCCATTCGCGTTGCGAAACGTCAATTTATGTTAATGTTGTGTGCCGGGTGCAATg                                                                                                   <  1:2640392/90‑1 (MQ=255)
            aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                         <  2:552298/100‑1 (MQ=255)
            aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  1:3197135/100‑1 (MQ=38)
              gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                       <  1:2223026/100‑1 (MQ=255)
              gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcacac                                                                                       <  2:4084202/100‑1 (MQ=38)
              gTGACCAAAGAGAAACCTGCACATCTAGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                       >  2:2678376/1‑100 (MQ=255)
               tGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacaca                                                                                      <  2:1253491/100‑1 (MQ=255)
                gACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACacaa                                                                                     <  2:3768138/100‑1 (MQ=38)
                   cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                  >  2:713640/1‑100 (MQ=255)
                    aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATc                                                                                 >  2:4375169/1‑100 (MQ=255)
                    aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATc                                                                                 >  1:377884/1‑100 (MQ=255)
                    aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGAGAAACGTTAATTTACGTTACTGTTCTGTGCCGGGTGCAATGCATCCGGCACcaaaca                                                                                   >  1:1075773/1‑98 (MQ=255)
                     aaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATca                                                                                >  2:4142735/1‑100 (MQ=255)
                       gagaACCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcaca                                                                              <  2:2044206/100‑1 (MQ=255)
                         gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            >  2:2826297/1‑100 (MQ=255)
                         gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            >  2:2835779/1‑100 (MQ=255)
                         gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGCGCAATGCATCCGGCACACAACATCACAcc                                                                            >  1:2422064/1‑100 (MQ=255)
                           aaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGACGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            >  1:2060057/1‑98 (MQ=255)
                            aCCTGCACATCTCGAGCAATCCATTCGCGTTGCGAAACGTTCATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCata                                                                         <  1:2561513/100‑1 (MQ=255)
                             ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCacac                                                                             <  1:2870509/95‑1 (MQ=255)
                             ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCatat                                                                        >  1:2110433/1‑100 (MQ=255)
                               tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTCCGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCatatac                                                                      <  1:2853014/100‑3 (MQ=255)
                               tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTTCCGGGTGCACTGCATCCGGCACACAACATCACACCATATcc                                                                      >  2:2765004/1‑100 (MQ=255)
                                     tCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCa                                                                >  1:1272265/1‑100 (MQ=255)
                                       tGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAg                                                               <  1:3924505/99‑1 (MQ=255)
                                        ggCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGGATCCGGCACACAACATCACACCATATCCAGCGCAGtt                                                             >  2:618004/1‑100 (MQ=255)
                                        ggCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGtt                                                             <  1:3582072/100‑1 (MQ=255)
                                         gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttt                                                            <  2:34745/100‑1 (MQ=255)
                                          cgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGtttt                                                           <  2:3904923/100‑1 (MQ=255)
                                            cAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTc                                                         <  2:4299137/100‑1 (MQ=255)
                                             aaTGCATTCGCTTTCCGAAACGTTAATTTACGTTAATGTTGTGTCCCGGTTGGAATGCTTCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        <  1:567882/100‑1 (MQ=255)
                                             aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        <  2:3536784/100‑1 (MQ=255)
                                              aTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCt                                                       <  2:1690603/100‑1 (MQ=255)
                                                 cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtttt                                                    <  2:4198719/100‑1 (MQ=255)
                                                   ttCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                  <  1:4375138/100‑1 (MQ=255)
                                                   ttCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGAAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                  <  2:767712/100‑1 (MQ=255)
                                                     cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                >  2:1054198/1‑100 (MQ=255)
                                                       cgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGcc                                              >  2:2253331/1‑100 (MQ=255)
                                                         ttGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCgg                                            <  2:2634460/100‑1 (MQ=255)
                                                          tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCAGCCGGCACACATCATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                   <  1:3612640/92‑1 (MQ=255)
                                                            ngAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                          >  2:3359976/2‑99 (MQ=255)
                                                                 cGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCt                                    <  2:1521764/100‑1 (MQ=255)
                                                                   ttAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATTCGGCACACAACATCACACCATATCCAGCGCAGTTGTTCCTTTTGGTGCCGGATATGCCTTa                                  >  1:3026572/1‑100 (MQ=255)
                                                                      aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                               >  1:777781/1‑100 (MQ=255)
                                                                       tttACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCa                              >  1:2872134/1‑100 (MQ=255)
                                                                         tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                            <  2:1272270/100‑1 (MQ=255)
                                                                         tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                            <  2:297310/100‑1 (MQ=255)
                                                                           cGTTAATGTTGTGTGCCGGGTGCAATTCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGGCGGATATGTCTTATCCAGCAt                          <  1:3072422/100‑1 (MQ=255)
                                                                           cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                          <  2:4048193/100‑1 (MQ=255)
                                                                                 tgttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTaa                    >  2:3416916/1‑100 (MQ=255)
                                                                                  gttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACAATATCCAGCGCAGTTTCTCCTTATAGCGCCCGATACTCCTTATCGAGAAtctct                      >  1:3038629/1‑93 (MQ=12)
                                                                                  gttgtGTGCCGGATGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  2:2885410/1‑100 (MQ=255)
                                                                                     gtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTcc                >  1:4165458/1‑100 (MQ=255)
                                                                                      tgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCg               >  2:2962286/1‑100 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  2:1507581/100‑1 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  1:3328326/100‑1 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  1:2127382/1‑100 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  2:1473983/1‑100 (MQ=255)
                                                                                       gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGACGTATATGCCTTATCCAGCATAGCTAATTCCGc              <  1:2885394/100‑1 (MQ=255)
                                                                                       gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCga              >  2:778982/1‑99 (MQ=38)
                                                                                       gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  2:3759071/1‑100 (MQ=38)
                                                                                       gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGACGGATATGCCTTATCCAGCATAGCTAATTCAGc              >  1:3659694/1‑100 (MQ=37)
                                                                                        tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             <  1:4235220/100‑1 (MQ=255)
                                                                                        tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             <  1:2372345/100‑1 (MQ=255)
                                                                                         gCCGGGTTCAATTCATCCGGCACACAACATCACACCATATCCAGCGCATTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             <  1:346209/99‑1 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             >  1:1744401/1‑99 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1647149/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1000667/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:48360/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:1199396/1‑100 (MQ=255)
                                                                                         gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            <  2:2860981/100‑1 (MQ=255)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:309485/1‑100 (MQ=38)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:3032210/1‑100 (MQ=38)
                                                                                         gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCCGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:1637200/1‑100 (MQ=38)
                                                                                          ccGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGATGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTga           >  2:3503992/1‑100 (MQ=38)
                                                                                           cGGGTGCAATGCATCCGGCACAAAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa          >  1:1674990/1‑100 (MQ=255)
                                                                                            gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  1:4243262/100‑1 (MQ=255)
                                                                                                gCAATGCATCCGGCACACAACATCACACCATATCCAGCGCCGTTTTTCCTTTTGGTGCCGGATCTGCCTTATCCAGCATATCTAATTCCGCTGAAGAAAg     >  2:2507457/1‑100 (MQ=255)
                                                                                                gCAAGGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAg     <  2:391583/100‑1 (MQ=255)
                                                                                                  aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt   >  2:2314462/1‑100 (MQ=255)
                                                                                                   aTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTc  <  1:885831/100‑1 (MQ=255)
                                                                                                   aTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTc  <  1:3236337/100‑1 (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: