Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,868,859 Δ1 bp 16.7% intergenic (+17/‑96) 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,8590C.16.7% 291.0 / 58.2 102intergenic (+17/‑96)yeaG/yeaHprotein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base . (4/13);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 116 total aligned reads displayed.

GACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGA  >  NC_000913/1868759‑1868958
                                                                                                    |                                                                                                    
   cGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTacac                                                                                                    <  2:1563201/100‑1 (MQ=255)
   cGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTacac                                                                                                    <  1:1637193/100‑1 (MQ=255)
    ntATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTacac                                                                                                    >  2:2484965/2‑99 (MQ=255)
    gTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAcaca                                                                                                   <  2:2219113/100‑1 (MQ=255)
    gTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAcaca                                                                                                   <  2:1702835/100‑1 (MQ=255)
    gTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAcaca                                                                                                   <  2:2436232/100‑1 (MQ=255)
    gTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTAcaca                                                                                                   >  1:394308/1‑100 (MQ=255)
      atgatgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATg                                                                                                 <  1:2340800/100‑1 (MQ=255)
        gatgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGt‑t                                                                                              >  2:2023359/1‑100 (MQ=255)
         ntgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑Ta                                                                                             >  1:2512588/2‑100 (MQ=255)
         atgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACAAGT‑Ta                                                                                             >  1:911428/1‑100 (MQ=255)
           tggaaaaGGCTACACCCGTAAACAGGTGCGGTTACTGTGCGAATGGTATGTGCGCGTAAGTAAATCGTCTTAACAACCCTGGCCCGGTACACAGGT‑TAcc                                                                                           >  1:152254/4‑100 (MQ=255)
           ggAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TAcc                                                                                           <  2:2631081/100‑1 (MQ=255)
             aaaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCgg                                                                                         <  2:817317/100‑1 (MQ=255)
             aaaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCgg                                                                                         <  1:838167/100‑1 (MQ=255)
              aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCggg                                                                                        <  2:2143469/100‑1 (MQ=255)
              aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                        <  1:360490/100‑1 (MQ=255)
              aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                        <  1:1929291/100‑1 (MQ=255)
              aaaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                        <  2:1493882/100‑1 (MQ=255)
              aaaaGGCGACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCggg                                                                                        <  2:1934922/100‑1 (MQ=255)
               aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTGAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                       <  1:136978/100‑1 (MQ=255)
               aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGc                                                                                       <  2:2124642/100‑1 (MQ=255)
               aaaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                       <  2:1075232/100‑1 (MQ=255)
                aaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGcc                                                                                      <  2:622123/100‑1 (MQ=255)
                aaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGcc                                                                                      >  2:1015294/1‑100 (MQ=255)
                aaGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcc                                                                                      <  2:2749932/100‑1 (MQ=255)
                 aGGCTACGCCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGcca                                                                                     <  2:2610969/100‑2 (MQ=255)
                 aGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcca                                                                                     <  2:2675887/100‑2 (MQ=39)
                 aGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcca                                                                                     <  2:753900/100‑2 (MQ=39)
                  ggCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTa                                                                                    <  1:2152199/100‑1 (MQ=255)
                   gCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagn                                                                                   <  1:2351434/100‑4 (MQ=38)
                    cTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTCACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACa                                                                                  >  2:2092900/1‑100 (MQ=255)
                     tccaaccGTAAACAGGTGCGTTTCCTGGGCGAATGGTAGTTGCGCGAACGTAAATCGTCTTAATAACCCTGGCCCGGTA‑ACATGTGTACCGGGccaggg                                                                                  <  1:1941088/94‑5 (MQ=21)
                      acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAAc                                                                                <  1:411576/100‑1 (MQ=255)
                      acacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAAc                                                                                >  1:1900449/1‑100 (MQ=255)
                       cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCtgggttn                                                                               <  2:2512603/100‑7 (MQ=255)
                       cacCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACg                                                                               <  2:271935/100‑1 (MQ=255)
                         cccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGa                                                                              >  2:2220268/1‑99 (MQ=255)
                          ccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACa                                                                            >  1:984652/1‑100 (MQ=255)
                          ccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGccagggttgtta                                                                            <  1:2307917/100‑11 (MQ=25)
                            gTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGc                                                                          <  2:2188566/100‑1 (MQ=255)
                            gTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGc                                                                          <  2:2346567/100‑1 (MQ=255)
                              aaaCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGa                                                                        <  2:271924/100‑1 (MQ=255)
                                aCAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAAc                                                                      <  1:2777795/100‑1 (MQ=255)
                                   ggTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGt                                                                   >  2:580277/1‑100 (MQ=255)
                                     tGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTgg                                                                 <  2:78354/100‑1 (MQ=255)
                                      gCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTggg                                                                >  2:2786498/1‑100 (MQ=255)
                                       cGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGc                                                               >  2:2226884/1‑100 (MQ=255)
                                          ttACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTg                                                            <  2:2836442/100‑1 (MQ=255)
                                           tACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTga                                                           >  2:2507137/1‑100 (MQ=255)
                                            aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgag                                                          >  2:2394936/1‑100 (MQ=255)
                                            aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAACGAACCGTGGGCCTgag                                                          >  1:957189/1‑100 (MQ=255)
                                             cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgaga                                                         <  2:144442/100‑1 (MQ=255)
                                             cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgaga                                                         >  1:887364/1‑100 (MQ=255)
                                               gtgCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAg                                                       >  2:2674437/1‑100 (MQ=255)
                                                   gAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                   <  2:636639/100‑1 (MQ=255)
                                                   gAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGc                                                   >  2:2134187/1‑100 (MQ=255)
                                                        ntATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAacac                                              >  1:2539247/2‑100 (MQ=255)
                                                        gTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAacac                                              >  2:412683/1‑100 (MQ=255)
                                                             tGCGCGTACGTAAATCGTCTTAAAAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCg                                         >  2:2471103/1‑100 (MQ=255)
                                                               cgcgTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTa                                       >  2:658464/1‑100 (MQ=255)
                                                               cgcgTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTa                                       <  1:533228/100‑1 (MQ=255)
                                                                 cgtacgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGCCAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGc                                     <  1:2263318/100‑1 (MQ=255)
                                                                  ttacgtaAATGGTCATAACAACCCTGGCCCGGTACAAATGT‑TATCGGGCTTACAACGACAGCCAACCGTGGGCGTGAGAAGCGGCAACACAGGCGt                                        >  1:1183959/2‑96 (MQ=255)
                                                                  gtacgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCa                                    <  1:1086201/100‑1 (MQ=255)
                                                                   tacgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCAt                                   >  2:34348/1‑100 (MQ=255)
                                                                     cgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATAc                                 <  1:2436513/100‑1 (MQ=255)
                                                                     cgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATAc                                 >  1:1446978/1‑100 (MQ=255)
                                                                            cGTCTTAACGGCCCTGGCCCGGTAAACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATCCAGTTGGc                          <  1:551100/100‑1 (MQ=255)
                                                                             gTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGc                          >  2:1735604/1‑99 (MQ=255)
                                                                             gTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCa                         >  1:1047376/1‑100 (MQ=255)
                                                                             gTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGACCTGAGCAGCGGCAACACAGGCGTAGCATACAGTTGGCa                         >  1:2049135/1‑100 (MQ=255)
                                                                              tCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaa                        >  2:465913/1‑100 (MQ=255)
                                                                               cTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCaaa                       <  2:342613/100‑1 (MQ=255)
                                                                                 tAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATg                     >  2:708447/1‑100 (MQ=255)
                                                                                 tAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGAGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTTGCAAAt                      >  1:968707/1‑99 (MQ=255)
                                                                                 tAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTAAAACCACAGCGAAACGTGGGCCTGAGAAGCGGCAACACCGGCGTAGCATACAGTTGGCaaa                       >  1:559284/1‑98 (MQ=255)
                                                                                 tAACAACCATGGCCCGGTACACATGT‑TACCGGGACTACAAAGACAGCGAACCGCGGGCCGGAGAGGCAGCAACAGAGGCGTAGCAGACAGTTg                            >  1:1395691/1‑93 (MQ=255)
                                                                                  aacaacCCTGGCCCGGTACACATGT‑TAGCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTTGCaaa                       <  1:1933346/97‑1 (MQ=255)
                                                                                  aacaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgt                    <  1:1343077/100‑1 (MQ=255)
                                                                                  aacaacCCTGGCCCGGGCCACAGGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATg                     <  1:2718515/99‑1 (MQ=255)
                                                                                   acaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATgta                   <  1:1404401/100‑1 (MQ=255)
                                                                                    caacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACAAAGGCGTAGCATACAGTTGGCAAATgtag                  >  1:1942485/1‑100 (MQ=255)
                                                                                         taggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgg             >  2:1248160/3‑100 (MQ=39)
                                                                                         cTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgg             >  2:1829682/1‑100 (MQ=255)
                                                                                         cTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgg             >  1:2504363/1‑100 (MQ=255)
                                                                                          tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggg            >  1:2318971/1‑100 (MQ=255)
                                                                                          tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggg            >  2:643582/1‑100 (MQ=255)
                                                                                          tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggg            >  1:1103176/1‑100 (MQ=255)
                                                                                           ggCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgggg           <  2:2049146/100‑1 (MQ=255)
                                                                                           ggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgggg           >  2:166425/1‑100 (MQ=255)
                                                                                            gCCCGGTACACATGT‑TACCGGGCCTACAACGCCAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          <  2:818555/100‑1 (MQ=255)
                                                                                            gCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  2:656812/1‑100 (MQ=255)
                                                                                            gCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACggggg          >  1:1091322/1‑100 (MQ=255)
                                                                                            gCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAACTGTAGTACggggg          >  1:497089/1‑100 (MQ=255)
                                                                                              ccGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGc        <  2:77836/100‑1 (MQ=255)
                                                                                               cGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCa       <  2:1777063/100‑1 (MQ=255)
                                                                                                   acacATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg   <  2:1390840/100‑1 (MQ=255)
                                                                                                   acacATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATg   <  1:1355449/100‑1 (MQ=255)
                                                                                                    cacaTGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGa  >  2:1174629/1‑100 (MQ=255)
                                                                                                    |                                                                                                    
GACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGA  >  NC_000913/1868759‑1868958

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