Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,763 A→T 15.9% 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.9% 338.2 / 55.7 139intergenic (+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 (51/65);  new base T (10/12);  total (62/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 144 total aligned reads displayed.

AGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864665‑1864862
                                                                                                  |                                                                                                   
aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                    >  2:2062075/1‑100 (MQ=255)
 gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATg                                                                                                   >  1:812694/1‑100 (MQ=255)
 gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATg                                                                                                   <  1:2265842/100‑1 (MQ=255)
 gAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATg                                                                                                   <  1:914231/100‑1 (MQ=255)
     cGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATc                                                                                               >  1:1535701/1‑100 (MQ=255)
     cGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATc                                                                                               >  2:3044303/1‑100 (MQ=255)
      gCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcn                                                                                              <  2:401088/100‑2 (MQ=255)
       cAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCg                                                                                             <  2:1193493/100‑1 (MQ=255)
         aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                           <  1:1601060/100‑1 (MQ=255)
         aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                           <  1:2507517/100‑1 (MQ=255)
         aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:1673588/100‑1 (MQ=38)
         aaaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGc                                                                                           <  2:2802662/100‑1 (MQ=38)
          aaaGTGACCGAAGAGAAGCCTGCACATCTGGCGCAATCCATTCGCGTTGGGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGc                                                                                           <  1:706846/99‑1 (MQ=255)
          aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGca                                                                                          <  2:3419230/100‑1 (MQ=38)
          aaaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGca                                                                                          <  1:1727786/100‑1 (MQ=38)
           aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcac                                                                                         <  1:394507/100‑1 (MQ=255)
           aaGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcac                                                                                         <  1:3594212/100‑1 (MQ=38)
            aGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGcaca                                                                                        <  1:2551026/100‑1 (MQ=255)
               gACCAAAGCGAAACCTGCACATCTGGCGCAATCCATACGCGTTGCGAAACGTCAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGcacaca                                                                                      >  1:3422652/1‑99 (MQ=255)
               gACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGCGCAATGCATTCGGCACacaa                                                                                     >  1:4204800/1‑100 (MQ=255)
                 ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTACATTCCGTTAATGTTATGTGCCGGGTGCAATGCATCCGGcacccatc                                                                                    >  1:3771092/1‑97 (MQ=255)
                 ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  2:3439108/1‑100 (MQ=255)
                 ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  1:4028641/1‑100 (MQ=255)
                 ccAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACacaaca                                                                                   >  2:3713388/1‑100 (MQ=255)
                  cAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAt                                                                                  <  1:511529/100‑1 (MQ=255)
                   aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATc                                                                                 <  1:1120859/100‑1 (MQ=255)
                   aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCAGGGTGCAATGCATGCGGCACACAACATc                                                                                 >  1:3003067/1‑100 (MQ=255)
                     agagAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcac                                                                               >  2:1238247/1‑100 (MQ=255)
                      gagaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACAGAACATcaca                                                                              >  1:1339861/1‑100 (MQ=255)
                      gagaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATcaca                                                                              >  2:3884725/1‑100 (MQ=255)
                       agaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCGGCACACAACAttaacg                                                                             <  2:2189848/100‑6 (MQ=25)
                        gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            <  2:2455479/100‑1 (MQ=255)
                        gaAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACAcc                                                                            <  1:807614/100‑1 (MQ=255)
                         aaaCCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCa                                                                           <  1:1393795/100‑1 (MQ=255)
                                  caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  <  2:2468083/100‑1 (MQ=255)
                                  caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  >  2:3306407/1‑100 (MQ=255)
                                  caTCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAgcg                                                                  <  2:3709707/100‑1 (MQ=255)
                                       ggCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCATCGCAGtt                                                             <  1:2308367/100‑1 (MQ=255)
                                        gcgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttt                                                            >  2:1685984/1‑100 (MQ=255)
                                         cgcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGtttt                                                           >  2:588918/1‑100 (MQ=255)
                                          gcAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGttttt                                                          <  1:2339737/100‑1 (MQ=255)
                                           cAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTc                                                         <  1:53156/100‑1 (MQ=255)
                                            aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        <  2:1231846/100‑1 (MQ=255)
                                            aaTCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTcc                                                        >  2:1474875/1‑100 (MQ=255)
                                                cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCtttt                                                    <  1:3854392/100‑1 (MQ=255)
                                                   tCGCGTTGCGAAGCGTTAATTTAGGTTAATGTTGTGTGCCGGGTGCGATGCATCCGGCACAGAACATCACACCATCTCCAGCGCAGTTTTTCCTTTTgg                                                  <  1:1093162/99‑1 (MQ=255)
                                                    cgctTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                <  1:3306548/100‑1 (MQ=255)
                                                    cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTg                                                >  1:1946218/1‑100 (MQ=255)
                                                    cgcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACAATATCCAGCGCAGTTTTTCCTTTTgg                                                  >  1:484235/1‑98 (MQ=255)
                                                     gcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGt                                                 >  1:3258118/1‑98 (MQ=255)
                                                     gcgTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGc                                               >  1:1878560/1‑100 (MQ=255)
                                                       gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                             <  2:830269/100‑1 (MQ=255)
                                                       gTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCg                                             <  2:310791/100‑1 (MQ=255)
                                                        ttGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCgg                                            <  2:2877017/100‑1 (MQ=255)
                                                          gctaaaCGTTAATTTACGTTAATGTTGTGTGCCGGGAGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGat                                          <  2:3100317/97‑1 (MQ=255)
                                                          gCGAAACGTTCATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGa                                           <  1:2593207/99‑1 (MQ=255)
                                                          gCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCAGAACACAACATCACACCATATCCAGCGCAGTTTTTCCGTTTGGTgaagga                                           >  1:3707722/1‑94 (MQ=255)
                                                            gAAACGTTAATTTACGTTAATGTTGGGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGatat                                        <  1:533456/100‑1 (MQ=255)
                                                             aaaCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                       <  2:3241016/100‑1 (MQ=255)
                                                             aaaCGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                       <  2:613068/100‑1 (MQ=255)
                                                                cGTTAATTTACGTTAATGTTGTGTGCCGGGTGCACTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTCTTGGTGCCGGATATGCCt                                    >  1:1086295/1‑100 (MQ=255)
                                                                cGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCCGTTTTTCCTTTTGGTGCCGGATATGc                                      >  1:220527/1‑98 (MQ=255)
                                                                cGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCt                                    <  2:2568880/100‑1 (MQ=255)
                                                                 gTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCtt                                   >  1:1283985/1‑100 (MQ=255)
                                                                  ttAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTa                                  <  2:2935229/100‑1 (MQ=255)
                                                                     aTTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATcc                               >  2:2997763/1‑100 (MQ=255)
                                                                       ttACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAg                             <  2:3390569/100‑1 (MQ=255)
                                                                        tACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                            <  2:1405467/100‑1 (MQ=255)
                                                                          cGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                          >  2:3208259/1‑100 (MQ=255)
                                                                          cGTTAATGTCGTGTGCCGGGTGCAATGCATCCGCCACCCAAAATCACACCATATCCAGCCCAGTTTTTCCTTTTTGTGCCGGATATGCCTTATCCAGCAt                          <  1:2705798/100‑1 (MQ=255)
                                                                           gTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATa                         >  1:2119565/1‑100 (MQ=255)
                                                                            ttAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAg                        >  2:2599924/1‑100 (MQ=255)
                                                                            gtgaTGTTGTGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAg                        >  2:1100853/4‑100 (MQ=25)
                                                                              aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                      <  1:2039696/100‑1 (MQ=255)
                                                                              aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCt                      <  2:386248/100‑1 (MQ=255)
                                                                              aaTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCAAATCCATCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGc                       >  1:359377/1‑99 (MQ=255)
                                                                               aTGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTa                     >  2:4159360/1‑100 (MQ=255)
                                                                                 tttgtGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   <  1:3061676/99‑1 (MQ=255)
                                                                                 gttgtGTGCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAAt                   >  2:211080/1‑100 (MQ=38)
                                                                                   tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                 >  2:2020061/1‑100 (MQ=255)
                                                                                   tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                 <  2:2715420/100‑1 (MQ=255)
                                                                                   tgtgtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                 >  2:236275/1‑100 (MQ=255)
                                                                                    gtgtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTcc                >  2:968284/1‑100 (MQ=38)
                                                                                      gtgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  2:1036063/100‑1 (MQ=255)
                                                                                      gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              <  1:1100850/100‑1 (MQ=38)
                                                                                      gtgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTCTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc              >  2:1450698/1‑100 (MQ=38)
                                                                                       tgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             >  2:3329867/1‑100 (MQ=38)
                                                                                       tgCCGGATGCATTGCACCCGACACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCt             >  2:1325441/1‑100 (MQ=38)
                                                                                        gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  2:2264618/1‑100 (MQ=255)
                                                                                        gCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            <  1:605507/100‑1 (MQ=255)
                                                                                        gCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTg            >  1:4214002/1‑100 (MQ=38)
                                                                                           gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:1946225/100‑1 (MQ=255)
                                                                                           gggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaag         <  2:137683/100‑1 (MQ=255)
                                                                                            ggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagt        >  1:892656/1‑99 (MQ=255)
                                                                                            ggTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaaga        >  2:394865/1‑100 (MQ=255)
                                                                                             gTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTgaagaa       <  2:2671172/100‑1 (MQ=255)
                                                                                                cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGATAATTCCGCTGAAGAAAg     >  2:3678530/1‑99 (MQ=255)
                                                                                                cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTCTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGTAGAAAg     >  1:664501/1‑99 (MQ=255)
                                                                                                 aaTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGtt   >  1:2540705/1‑100 (MQ=255)
                                                                                                  aTGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTc  >  2:4175635/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTC  >  NC_000913/1864665‑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: