Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,768 T→C 16.1% intergenic (+34/+14) 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,7680TC16.1% 212.8 / 32.0 93intergenic (+34/+14)yeaD/yeaED‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (40/38);  new base C (6/9);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCT  >  NC_000913/1864672‑1864867
                                                                                                |                                                                                                   
cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                   >  1:1192996/1‑99 (MQ=255)
cAAAAAGTCACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCTAAACGTTAATTTGCGTTAATATTGTGTGCCGGGTGCGATGCATCCg                                                                                                  >  1:737640/1‑100 (MQ=255)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                                >  1:1193996/1‑100 (MQ=255)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                                <  2:1118843/100‑1 (MQ=255)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                                <  1:1618781/100‑1 (MQ=255)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                                <  2:203149/100‑1 (MQ=255)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                                <  1:12059/100‑1 (MQ=38)
  aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                                <  1:163008/100‑1 (MQ=38)
    acgTGACCATAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTTATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                              <  2:2007953/98‑1 (MQ=37)
    aaGTGACCAAAGAGAACCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                              <  1:2351206/100‑1 (MQ=38)
    aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTTGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                              <  1:134403/100‑1 (MQ=38)
    aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                              <  1:2370970/100‑1 (MQ=255)
      gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacac                                                                                            >  1:1834445/1‑100 (MQ=255)
      gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcacac                                                                                            <  2:2000231/100‑1 (MQ=38)
      gTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcaca                                                                                             <  2:2252274/99‑1 (MQ=255)
         aCCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaac                                                                                         <  2:1834438/100‑1 (MQ=255)
          ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                        >  2:2330000/1‑100 (MQ=255)
          ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                        >  2:2494766/1‑100 (MQ=255)
          ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                        >  2:325658/1‑100 (MQ=255)
           cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATTTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                        >  2:566726/1‑99 (MQ=255)
           cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                       >  2:560957/1‑100 (MQ=255)
           cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                       >  1:134022/1‑100 (MQ=255)
           cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACAt                                                                                       >  1:1304024/1‑100 (MQ=255)
           cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAt                                                                                       <  2:160357/100‑1 (MQ=37)
           cAAAGAGAAACCAGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAt                                                                                       <  1:2384724/100‑1 (MQ=37)
            aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATc                                                                                      >  2:1213212/1‑100 (MQ=255)
               gagaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcaca                                                                                   >  1:2536114/1‑100 (MQ=255)
                 gaAACCTGCACATCTGGCGCAATCCATTCTCGGTGCTAAACGTTAATTTACGTTAACATGGTGTGCCGGGTGCAATGCATCCGGCACACAACATCacac                                                                                  <  2:1355385/99‑1 (MQ=255)
                  aaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCa                                                                                <  1:227234/100‑1 (MQ=255)
                    aCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCata                                                                              <  2:1124401/100‑1 (MQ=255)
                     ccTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCata                                                                              >  2:1881155/1‑99 (MQ=255)
                      cTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATc                                                                            >  1:158811/1‑100 (MQ=255)
                       tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCCCACCATATcc                                                                           >  2:1516579/1‑100 (MQ=255)
                       tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                           >  1:524370/1‑100 (MQ=255)
                       tGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATcc                                                                           <  1:2268870/100‑1 (MQ=255)
                               tGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTTCAATGCATCCGACACACAACATCACACCATATCCAGCGCa                                                                     >  1:968703/1‑98 (MQ=255)
                                    cAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTc                                                              <  2:1226510/100‑1 (MQ=255)
                                     aaTCCATTCGCGTTGCGAAACTTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                             <  1:2622875/100‑1 (MQ=255)
                                     aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                             <  1:2398449/100‑1 (MQ=255)
                                     aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                             >  1:1173542/1‑100 (MQ=255)
                                         cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtttt                                                         <  2:2078773/100‑1 (MQ=255)
                                         cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtttt                                                         >  2:547337/1‑100 (MQ=255)
                                          aTTCGCGTTGCGAAACGTTCATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAATATCACCCCATATCCAGCGCAGTTTTTCCttt                                                          >  2:1849099/1‑98 (MQ=255)
                                          aTTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTg                                                        >  1:9178/1‑100 (MQ=255)
                                           ttCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTgg                                                       <  2:1824985/100‑1 (MQ=255)
                                               cgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGcc                                                   <  2:1482807/100‑1 (MQ=255)
                                                gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGccn                                                  <  2:2064488/100‑2 (MQ=255)
                                                  tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                                <  2:47446/100‑1 (MQ=255)
                                                  tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                                <  2:1145831/100‑1 (MQ=255)
                                                  tGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                                >  2:1520298/1‑100 (MQ=255)
                                                   gCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                               >  2:2001233/1‑100 (MQ=255)
                                                   gCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                               >  1:2503001/1‑100 (MQ=255)
                                                     gAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGatat                                             <  1:2224073/100‑1 (MQ=255)
                                                      aaaCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                            <  1:1896767/100‑1 (MQ=255)
                                                         cGTTAATTTACGTTAATGTTGTGCGCCGGGTGCAATGCATCCTGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCt                                         <  2:1392113/100‑1 (MQ=255)
                                                             aaTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTa                                       <  2:2186210/98‑1 (MQ=255)
                                                              aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                                    <  1:2246931/100‑1 (MQ=255)
                                                              aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                                    <  2:2309107/100‑1 (MQ=255)
                                                                ttACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAg                                  <  1:2001245/100‑1 (MQ=255)
                                                                 tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                                 <  1:1213217/100‑1 (MQ=255)
                                                                   cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                               >  1:2421302/1‑100 (MQ=255)
                                                                     ttAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAg                             <  2:2262540/100‑1 (MQ=255)
                                                                       aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                           >  2:1306633/1‑100 (MQ=255)
                                                                        aTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTa                          <  1:560958/100‑1 (MQ=255)
                                                                        aTGTTGTGTGCCCGGTGCAATGCATCCGGCACACAACAGCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTa                          <  2:2512941/100‑1 (MQ=255)
                                                                         tgttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTaa                         >  2:2718498/1‑100 (MQ=255)
                                                                          gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                        >  1:1686696/1‑100 (MQ=255)
                                                                          gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                        <  2:1239789/100‑1 (MQ=255)
                                                                          gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                        >  1:596734/1‑100 (MQ=255)
                                                                          gttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                        >  2:2108340/1‑100 (MQ=255)
                                                                            tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                      >  1:1997137/1‑100 (MQ=255)
                                                                            tgtgtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATACCTAATTc                      >  2:2521339/1‑100 (MQ=38)
                                                                             gtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTcc                     <  1:145110/100‑1 (MQ=255)
                                                                              tgtgCCGGATGCATTGCACCAGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                      >  2:1674441/1‑98 (MQ=38)
                                                                               gtgCCGGGTGCAATGCATCCGGCACACAACATTACACCACATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc                   <  2:1931088/100‑1 (MQ=255)
                                                                                 gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg                 <  1:2330015/100‑1 (MQ=255)
                                                                                 gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg                 >  1:749925/1‑100 (MQ=38)
                                                                                 gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg                 >  2:2234885/1‑100 (MQ=38)
                                                                                 gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg                 >  1:1566185/1‑100 (MQ=38)
                                                                                 gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg                 >  2:303000/1‑100 (MQ=38)
                                                                                   cGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaa               >  1:1982585/1‑100 (MQ=255)
                                                                                      gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagaa            >  2:1003226/1‑100 (MQ=255)
                                                                                         cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGt         >  1:2494905/1‑100 (MQ=255)
                                                                                          aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt        >  1:579886/1‑100 (MQ=255)
                                                                                          aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCATTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt        >  2:2083977/1‑100 (MQ=255)
                                                                                            tGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCg      <  1:1881163/100‑1 (MQ=255)
                                                                                             gCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGa     >  1:2089076/1‑100 (MQ=255)
                                                                                             gCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGa     >  2:1130131/1‑100 (MQ=255)
                                                                                              cATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGAc    >  2:297822/1‑100 (MQ=255)
                                                                                              cATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGAc    <  1:987074/100‑1 (MQ=255)
                                                                                               aTCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGAcc   <  2:299227/100‑1 (MQ=255)
                                                                                               aTCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGAcc   <  1:802457/100‑1 (MQ=255)
                                                                                               aTCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGAcc   <  2:795401/100‑1 (MQ=255)
                                                                                                tCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCt  >  2:1746074/1‑100 (MQ=255)
                                                                                                |                                                                                                   
CAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCT  >  NC_000913/1864672‑1864867

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