Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,868,865:1 +G 21.3% intergenic (+23/‑90) yeaG → / → yeaH protein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,868,8651.G21.3% 246.7 / 83.4 94intergenic (+23/‑90)yeaG/yeaHprotein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein
Reads supporting (aligned to +/- strand):  ref base . (43/31);  new base G (9/11);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

GATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGT  >  NC_000913/1868767‑1868964
                                                                                                   |                                                                                                   
gatgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGt‑t                                                                                                    >  1:910773/1‑100 (MQ=255)
 atgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑Ta                                                                                                   >  2:372370/1‑100 (MQ=255)
  tgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TAc                                                                                                  >  2:576941/1‑100 (MQ=255)
    naaaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCg                                                                                                >  1:1531075/2‑100 (MQ=255)
      aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                              <  2:2563121/100‑1 (MQ=255)
       tctggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                             <  1:2098237/97‑1 (MQ=255)
       aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGc                                                                                             <  2:448368/100‑1 (MQ=255)
       aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGc                                                                                             >  1:486337/1‑100 (MQ=255)
       aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                              <  1:1425417/99‑1 (MQ=255)
       aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                             <  1:2360157/100‑1 (MQ=255)
         aGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcca                                                                                           <  2:1313878/100‑2 (MQ=39)
          ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTa                                                                                          <  2:2085919/100‑1 (MQ=255)
          ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccag                                                                                          <  2:2104057/100‑3 (MQ=39)
          ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcaga                                                                                          >  2:2098251/1‑97 (MQ=255)
           gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTAc                                                                                         <  2:562952/100‑1 (MQ=255)
           gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagg                                                                                         <  2:1453311/100‑4 (MQ=38)
           gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccag                                                                                          <  2:2360815/99‑3 (MQ=38)
            cTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACCTGT‑TACCGGGCCTaaa                                                                                        >  1:1249426/1‑98 (MQ=255)
             tACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGTGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagggt                                                                                       <  2:2493163/100‑6 (MQ=37)
             tACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagggt                                                                                       <  2:132415/100‑6 (MQ=38)
              acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACCAc                                                                                      >  1:2031577/1‑100 (MQ=255)
              acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAAc                                                                                      <  2:2132023/100‑1 (MQ=255)
              acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTCGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagggt                                                                                       <  2:783663/99‑6 (MQ=37)
               cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACg                                                                                     >  2:91665/1‑100 (MQ=255)
               cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCACGGTACACATGT‑TACCGGGCCTACAACg                                                                                     >  1:1064488/1‑100 (MQ=255)
               cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGAGCGTCCGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACCACg                                                                                     >  1:1149676/1‑100 (MQ=255)
                 cccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGAc                                                                                   >  2:531575/1‑100 (MQ=255)
                 cccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGAc                                                                                   >  2:1016016/1‑100 (MQ=255)
                  ccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACa                                                                                  >  2:1925694/1‑100 (MQ=255)
                      aaaCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGa                                                                              >  1:1867197/1‑100 (MQ=255)
                       aaCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGaa                                                                             >  1:34414/1‑100 (MQ=255)
                          aGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCg                                                                          <  1:714808/100‑1 (MQ=255)
                           ggTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCg                                                                          <  1:1128414/99‑1 (MQ=255)
                           ggTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCg                                                                          >  2:2359081/1‑99 (MQ=255)
                             tGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTgg                                                                       <  1:333702/100‑1 (MQ=255)
                              gCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTggg                                                                      <  2:1460080/100‑1 (MQ=255)
                               cGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGc                                                                     <  1:667650/100‑1 (MQ=255)
                                gTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                                    <  1:1316312/100‑1 (MQ=255)
                                gTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                                    <  2:699804/100‑1 (MQ=255)
                                  ttACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTg                                                                  <  1:2473489/100‑1 (MQ=255)
                                   tACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGcc                                                                    >  1:2024308/1‑97 (MQ=255)
                                    aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAGCGACAGCGAACCGTGGGCCTgag                                                                <  1:774688/100‑1 (MQ=255)
                                     cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑GACCGGGCCTACAACGACAGCGAACCGTGGGCCTgaga                                                               <  1:1389094/100‑1 (MQ=255)
                                     cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGA‑TACCGGGCCTACAACGACAGCGAACCGTGGGCATgag                                                                >  1:688329/1‑99 (MQ=255)
                                       gtgCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAg                                                             >  2:949295/1‑100 (MQ=255)
                                          cGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCgg                                                          >  2:989774/1‑100 (MQ=255)
                                           gAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                         <  2:1637298/100‑1 (MQ=255)
                                            aaTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                         <  1:2072693/99‑1 (MQ=255)
                                            aaTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCa                                                        <  1:944057/100‑1 (MQ=255)
                                             aTGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCaa                                                       >  1:1987875/1‑100 (MQ=255)
                                               ggTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAaca                                                     >  2:492509/1‑100 (MQ=255)
                                               ggTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAaca                                                     >  2:85142/1‑100 (MQ=255)
                                                 tATTTGCGCGTAGGTAAATCGTCTCAACAACCCGGGCCCGGTACACACGT‑TACAGGGCCTACTACGATAGCGAACCGTGGGCCTGAGATGCGGCAac                                                      <  1:1152545/97‑1 (MQ=255)
                                                 tATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAAcaca                                                   <  2:486339/100‑1 (MQ=255)
                                                  aTTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAg                                                  <  2:1064492/100‑1 (MQ=255)
                                                    tttcgcgTACGTAAATCGTCTTAACAACCCTGGCCCTGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGACAAGCGGCAACACAGGc                                                >  1:130637/4‑100 (MQ=255)
                                                     tGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCg                                               >  1:343040/1‑100 (MQ=255)
                                                      gcgcGTACGAAAATCGGCTTAACAACCCTGGCGCGGTACACATGA‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGGGGCAACACAgg                                                 >  1:863986/1‑97 (MQ=255)
                                                             cgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATAc                                       >  2:642951/1‑100 (MQ=255)
                                                              gtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                      >  2:488216/1‑100 (MQ=255)
                                                              gtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                      >  2:2241710/1‑100 (MQ=255)
                                                              gtaAATCGTCTTAACAACCCTGGCACGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACa                                      <  1:1025376/100‑1 (MQ=255)
                                                               taAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAg                                     <  1:1086117/100‑1 (MQ=255)
                                                                aaaTCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGt                                    >  2:828887/1‑100 (MQ=255)
                                                                    cGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGc                                >  1:1012590/1‑100 (MQ=255)
                                                                     gTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCa                               <  2:1993035/100‑1 (MQ=255)
                                                                     gTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCa                               <  1:1762085/100‑1 (MQ=255)
                                                                      tCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATCCAGTTGGCaa                              >  2:1171803/1‑100 (MQ=255)
                                                                       cTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaaa                             >  1:1394375/1‑100 (MQ=255)
                                                                        ttAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaact                            >  1:1571304/1‑98 (MQ=255)
                                                                         tAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGCCCGCGAACCGTGGGCCTGAGAAGCGGCCACACAGGCGTCGCATACAGTTGGCaaagg                           >  1:347899/1‑98 (MQ=255)
                                                                         tAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATg                           >  2:792941/1‑100 (MQ=255)
                                                                          aacaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgt                          <  1:1058034/100‑1 (MQ=255)
                                                                           acaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgta                         <  2:2510039/100‑1 (MQ=255)
                                                                              acCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgtagta                      >  2:2031246/1‑100 (MQ=255)
                                                                              acCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgtagta                      <  1:942591/100‑1 (MQ=255)
                                                                                gtaggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACg                    >  2:2265158/4‑100 (MQ=38)
                                                                                  aggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgttcaatcag                  >  2:1949632/2‑92 (MQ=38)
                                                                                    gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg                >  2:1922072/1‑100 (MQ=255)
                                                                                    gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg                >  2:918545/1‑100 (MQ=255)
                                                                                    gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg                >  1:344538/1‑100 (MQ=255)
                                                                                    gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg                >  2:245938/1‑100 (MQ=255)
                                                                                    gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg                >  2:763892/1‑100 (MQ=255)
                                                                                     cccGGTA‑ACATGTGTACCGGGCCTACAACGACCGCGAACCGTGGGCCTGAGACGCGGCATCACAGGCGTAGCATACAGCTGCCACATGTAGTACggggg                >  1:2205765/1‑99 (MQ=255)
                                                                                      ccGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGc              <  2:2140856/100‑1 (MQ=255)
                                                                                        ggTACACATGT‑TACCGGGCCTACAACGCCAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCat            >  2:680203/1‑100 (MQ=255)
                                                                                         gTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCata           <  2:34415/100‑1 (MQ=255)
                                                                                         gTACACATGT‑GACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCata           <  1:1978934/100‑1 (MQ=255)
                                                                                           acacATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg         <  2:2532374/100‑1 (MQ=255)
                                                                                           acacATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg         >  2:613843/1‑100 (MQ=255)
                                                                                            cacaTGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGa        >  2:34439/1‑100 (MQ=255)
                                                                                            cacaTGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGa        >  1:1220018/1‑100 (MQ=255)
                                                                                               aTGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCt     >  1:155514/1‑100 (MQ=255)
                                                                                                tGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGACAAATGTAGTACGGGGGGCATATGACCTg    >  1:910938/1‑100 (MQ=255)
                                                                                                 gT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTgg   >  2:1488441/1‑100 (MQ=255)
                                                                                                  t‑tACCGGGCCTACAACGCCAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGt  <  1:1294001/100‑1 (MQ=255)
                                                                                                  t‑tACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTggn  <  1:1338839/100‑2 (MQ=255)
                                                                                                   |                                                                                                   
GATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGT  >  NC_000913/1868767‑1868964

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: ‑