Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,879 (T)8→7 100% intergenic (‑159/‑232) uspB ← / → uspA universal stress (ethanol tolerance) protein B/universal stress global response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,8720T.95.6% 657.1 / 15.7 155intergenic (‑152/‑239)uspB/uspAuniversal stress (ethanol tolerance) protein B/universal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (74/74);  total (80/75)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 176 total aligned reads displayed.

ATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTC  >  NC_000913/3639775‑3639967
                                                                                                 |                                                                                               
aTACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATACTGTAGGAGATCAGGTTGttt                                                                                               <  2:2385492/100‑1 (MQ=255)
 tACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttt                                                                                              <  2:1127879/100‑1 (MQ=255)
 tACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttt                                                                                              <  2:2715270/100‑1 (MQ=255)
  aCGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGttttt                                                                                             <  1:3737566/100‑1 (MQ=255)
   cGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttttt                                                                                            <  1:3700130/100‑1 (MQ=255)
    gCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGttttttt                                                                                           <  1:3713681/100‑1 (MQ=255)
     cTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATAAGGTTGtttttttc                                                                                          >  1:2620545/1‑99 (MQ=255)
      tGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttttttcc                                                                                         <  1:1997035/100‑3 (MQ=255)
      tGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttttttcc                                                                                         <  1:592107/100‑3 (MQ=255)
       gCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttttttcca                                                                                        <  1:35354/100‑4 (MQ=255)
         cAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGa                                                                                     <  2:3294255/100‑1 (MQ=255)
         cAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGa                                                                                     >  1:1229722/1‑100 (MQ=255)
          aaTCATTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGTTAACCAACATACGAGCGGCTTTATAGATAGTGTATGAGATCAGGTTGTTTTTTTCCAGa                                                                                     <  2:3294715/99‑1 (MQ=255)
          aaTCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttttttccag                                                                                       <  1:1899074/98‑5 (MQ=255)
           aTCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGttttttt                                                                                           >  2:236609/1‑93 (MQ=255)
            tCAGTTAACACCAGGTCCTTGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAgg                                                                                  <  1:1597154/100‑1 (MQ=255)
            tCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATTGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAgg                                                                                  <  2:2086158/100‑1 (MQ=255)
            tCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAgg                                                                                  >  1:1274018/1‑100 (MQ=255)
            tCAGTTAACAACAGGTCCTGGAGAAACCTCTGTGGTGATGACCAACACACGAGCGGCTATGTAGATAGTGTAGAAGATCAGGTTCtttattcgcaca                                                                                      >  2:1528392/1‑91 (MQ=255)
             cAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGt                                                                                 >  2:3299028/1‑100 (MQ=255)
             cAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAAATCAGGTTGTTTTTTTCCAGAAGGt                                                                                 >  2:1749114/1‑100 (MQ=255)
              aGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGtt                                                                                >  1:1927741/1‑100 (MQ=255)
               gTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTa                                                                               >  1:3879311/1‑100 (MQ=255)
                  aaCACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAAcc                                                                            <  2:1596704/100‑1 (MQ=255)
                    cacCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTATAGGAAATCAGGTTGTTTTTTTCCAGAAGGTTAACCa                                                                           >  1:645189/1‑99 (MQ=255)
                      ccAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTa                                                                        <  2:210053/100‑1 (MQ=255)
                      ccAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTa                                                                        <  2:442467/100‑1 (MQ=255)
                        aGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGTCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATc                                                                      <  1:2383177/100‑1 (MQ=255)
                            ccTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCTGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAata                                                                  <  1:3300842/100‑1 (MQ=255)
                            ccTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAata                                                                  >  1:1868897/1‑100 (MQ=255)
                            ccTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAata                                                                  <  2:1134523/100‑1 (MQ=255)
                            ccTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGATTTTTTCCAGAAGGTTAACCACTATCAata                                                                  >  1:1483513/1‑100 (MQ=255)
                                 ataaaCCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTc                                                             <  1:763534/98‑1 (MQ=255)
                                  gaAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATAttta                                                            >  2:2213652/1‑98 (MQ=255)
                                  gaAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCa                                                            >  1:1452541/1‑100 (MQ=255)
                                  gaAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCa                                                            <  1:3594344/100‑1 (MQ=255)
                                   aaaCCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCAt                                                           >  1:2694598/1‑100 (MQ=255)
                                    aaCCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATg                                                          >  2:918130/1‑100 (MQ=255)
                                      tccgatctgtgGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTc                                                        <  2:2667576/93‑1 (MQ=255)
                                      ccGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTc                                                        <  1:1238734/100‑1 (MQ=255)
                                       cGCTTTTGTGGTGACCAACATACGTGCGGCTCTATAGAGAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATg                                                          <  2:1693816/97‑1 (MQ=255)
                                       cGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGATTAACCAATATCAATATATTCATGTCg                                                       >  1:1115684/1‑100 (MQ=255)
                                        gCTTTTGTGGTGACCAACATACGAGCGGCTCTATAAATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTACCCACTATCAATATATTCATGTCg                                                       <  2:2545315/99‑1 (MQ=255)
                                             tgtgGTGACCAACATACGAGCGGCTCTCTAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCATAAGGTTAAACACTATCAATATATTCATGTCGaaa                                                    >  1:2602696/1‑97 (MQ=255)
                                              gtggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAAAAGGTTAACCACTATCAATATATTCATGTCGCAAAtt                                                 >  2:3703985/1‑99 (MQ=255)
                                               tggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTATCCACTATCAATATATTCATGTCGAAAATTTg                                               <  2:1370483/100‑1 (MQ=255)
                                               tggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTg                                               <  1:1960890/100‑1 (MQ=255)
                                               tggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTg                                               <  2:3407356/100‑1 (MQ=255)
                                                ggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGt                                              >  2:1728758/1‑100 (MQ=255)
                                                ggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGt                                              >  1:3063728/1‑100 (MQ=255)
                                                ggtgACAAACATACGAGCGGCTCTATAGATATTGTAGGAGATCATGTTGTTTTTTTCCAGAATGTTAACCACTATCAATATATTCATGTCGAAAATTTGt                                              <  2:2980916/100‑1 (MQ=255)
                                                 gtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt                                             <  2:2689930/100‑1 (MQ=255)
                                                 gtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt                                             <  2:1886114/100‑1 (MQ=255)
                                                 gtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt                                             <  1:237541/100‑1 (MQ=255)
                                                 gtgAACAACATACGAGCGGCTATATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt                                             <  1:3526852/100‑1 (MQ=255)
                                                 ccgatctaCATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt                                             <  1:2849199/93‑1 (MQ=255)
                                                   gACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGACGGTTAACCACTATCAATATATTCATGTCGAAATTTTGTTTa                                           >  2:1696369/1‑100 (MQ=255)
                                                   gACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTa                                           <  1:3368947/100‑1 (MQ=255)
                                                    aCCAACATACGTGGTGCTCTATAGTTAGTGTAGGAGTTTAGGTTGTTTTTTTCCAGAAGGTTTACCTTTATCATTATATTCCTGTCGAAAATTTg                                               <  2:2422640/95‑1 (MQ=255)
                                                    aCCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTAt                                          <  1:2233415/100‑1 (MQ=255)
                                                      cAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCATGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATc                                         <  1:498001/99‑1 (MQ=255)
                                                      cAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCATAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTAtat                                        >  1:2549265/1‑98 (MQ=255)
                                                       aaCATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTa                                       >  2:600293/1‑100 (MQ=255)
                                                        ncATACGAGCGGCTCTATAGATAGAGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTTTTTATCt                                        >  2:1329320/2‑98 (MQ=255)
                                                        aCATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGttagatcgg                                      >  2:2849173/1‑93 (MQ=255)
                                                         cATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAAc                                     <  2:3051364/100‑1 (MQ=255)
                                                          aTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTTGAAAATTTGTTTATCTAACg                                    <  2:422854/100‑1 (MQ=255)
                                                          aTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACg                                    <  1:2175183/100‑1 (MQ=255)
                                                           tACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGa                                   >  1:1981424/1‑100 (MQ=255)
                                                            aCGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAg                                  >  1:3408774/1‑100 (MQ=255)
                                                              gAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTa                                <  1:3830494/100‑1 (MQ=255)
                                                               aGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTaa                               >  1:314200/1‑100 (MQ=255)
                                                                gCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAg                              <  1:2521874/100‑1 (MQ=255)
                                                                  ggCTCTATAGATAGTGTAGGAGATCAGTTTTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGt                                 <  1:129149/95‑1 (MQ=255)
                                                                   gCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCaa                           >  1:3750773/1‑100 (MQ=255)
                                                                      ctATAGATAGTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGc                        >  1:1101931/1‑100 (MQ=255)
                                                                       tataGATAGTGTAGGAGATCAGGTGGTTTTTTTCTAGAAGGGGAACCACTATCAATATATTCATGTCGAAAATTTGTTGATCTAACGAGTAAGCAAGGCg                       >  2:321557/1‑100 (MQ=255)
                                                                         tagatagTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGa                     <  2:3597327/100‑1 (MQ=255)
                                                                          agatagTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGAt                    >  1:1994453/1‑100 (MQ=255)
                                                                          agatagTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGAt                    <  1:3013529/100‑1 (MQ=255)
                                                                             tagTGTAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTAACTAACGAGTAAGCAAGGCGGATTGa                 >  1:546232/1‑100 (MQ=255)
                                                                               gtgtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGATAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACg               <  2:581215/100‑1 (MQ=255)
                                                                               gtgtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGAca               >  2:3398393/1‑99 (MQ=255)
                                                                               gtgtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACg               >  2:297705/1‑100 (MQ=255)
                                                                               gtgtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACg               <  2:3893007/100‑1 (MQ=255)
                                                                                 gtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGa             >  2:2008610/1‑100 (MQ=255)
                                                                                 gtAGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAAAGAGTAAGCAAGGCGGATTGACGGa             >  1:915005/1‑100 (MQ=255)
                                                                                  tAGGAGATCAGGGTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCAAAAATTTGTTTAACTAACGAGTAAGCAAGGCGGATTGACGGat            >  1:3274169/1‑100 (MQ=255)
                                                                                   aGGAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGatc           <  2:2063876/100‑1 (MQ=255)
                                                                                    ggAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGatca          >  1:1770354/1‑100 (MQ=255)
                                                                                    ggAGATCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGatca          >  2:717618/1‑100 (MQ=255)
                                                                                     gagaTCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTTACGAGTAAGCAAGGCGGATTGACGGatca          <  2:2823041/99‑1 (MQ=255)
                                                                                      agaTCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGatcatc        >  1:475108/1‑100 (MQ=255)
                                                                                        aTCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATATAAAGGGTAAGCAAGGCAAATTGACGGATCATcc       >  1:1529133/1‑100 (MQ=255)
                                                                                         tCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCgg     <  1:3647150/100‑1 (MQ=255)
                                                                                         tCAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCgg     <  2:3408736/100‑1 (MQ=255)
                                                                                          cAGGTTGTTTTTTTCCAGAAGGTTAACCACTATCTATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCggg    <  1:243613/100‑1 (MQ=255)
                                                                                           aGGTTGTTTTGTTCCAGAAAGTTAAAAACTATCAATATATAAATGTCGAAAATTTGTTTAACTAAAGAGTAAGAAAGGCGGATTGAAGAatcatcag      >  2:3498537/1‑95 (MQ=255)
                                                                                            ggTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTc  <  2:809478/100‑1 (MQ=255)
                                                                                            ggTTGTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATATAACGAGTAAGAAAAGCGGATTGACGGAGCATCCgggcc  >  1:3012813/1‑98 (MQ=255)
                                                                                                 |                                                                                               
ATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTC  >  NC_000913/3639775‑3639967

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