Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,250,210 T→A 100% V165E (GTG→GAG)  yeiH → UPF0324 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,250,2100TA95.4% 242.4 / 1.2 80V165E (GTG→GAG) yeiHUPF0324 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (36/39);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGTTGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTT  >  NC_000913/2250117‑2250306
                                                                                             |                                                                                                
gCTGGTTGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGACGCTg                                                                                             >  2:1109690/1‑99 (MQ=255)
  atcttGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTg                                                                                          <  2:1701567/97‑1 (MQ=255)
   ggTTGATCGTTGGCGGTAGCACTAGCTGTGTTGCTGCCGCGGTACTGGTGACTAAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTg                                                                                             <  2:1377162/96‑1 (MQ=255)
    gTTGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCa                                                                                        >  2:996949/1‑100 (MQ=255)
    gTTGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCa                                                                                        >  2:1827081/1‑100 (MQ=255)
      tGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTgagat                                                                                      >  1:1701579/1‑96 (MQ=255)
       gATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAAcc                                                                                     >  2:1465239/1‑100 (MQ=255)
        aTCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCg                                                                                    >  1:1528143/1‑100 (MQ=255)
         tCGGGGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGTCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgt                                                                                   <  2:1706546/100‑1 (MQ=255)
           ggTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttg                                                                                 <  1:1840651/100‑1 (MQ=255)
           ggTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttg                                                                                 <  1:422430/100‑1 (MQ=255)
            gTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgt                                                                                <  1:672363/100‑1 (MQ=255)
            gTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgt                                                                                <  2:1754907/100‑1 (MQ=255)
            gTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgt                                                                                <  1:2076999/100‑1 (MQ=255)
            gTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgt                                                                                <  1:972417/100‑1 (MQ=255)
             tGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgtt                                                                               >  2:623258/1‑100 (MQ=255)
             tGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCgttgtt                                                                               <  2:2517489/100‑1 (MQ=255)
              gCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTa                                                                              <  1:858635/100‑1 (MQ=255)
              gCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTa                                                                              >  2:1760500/1‑100 (MQ=255)
              gCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTa                                                                              <  1:2280790/100‑1 (MQ=255)
                    aGCAGAATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCg                                                                        >  1:1777383/1‑100 (MQ=255)
                      cAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCggg                                                                      >  1:2489291/1‑100 (MQ=255)
                              gtggtgCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTcgc                                                              <  1:13212/100‑1 (MQ=255)
                               tggtgCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCTGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTcgcg                                                             <  1:1688489/100‑1 (MQ=255)
                                 gtgCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGAt                                                           <  1:73627/100‑1 (MQ=255)
                                 gtgCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGAt                                                           <  1:2475565/100‑1 (MQ=255)
                                  tgctgcCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGAtt                                                          >  1:1676876/1‑100 (MQ=255)
                                      gcCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTcc                                                      <  1:2353928/100‑1 (MQ=255)
                                      gcCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTcc                                                      <  1:2099608/100‑1 (MQ=255)
                                      gcCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTcc                                                      >  2:1937841/1‑100 (MQ=255)
                                        cgcgGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAGGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCctc                                                    <  2:2493832/100‑1 (MQ=255)
                                        cgcgGTACTGGCCACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCctc                                                    >  1:1681479/1‑100 (MQ=255)
                                            gTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTAcc                                                >  1:120838/1‑100 (MQ=255)
                                            gTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTAcc                                                >  1:120858/1‑100 (MQ=255)
                                             tACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTAccc                                               <  1:1220366/100‑1 (MQ=255)
                                               cTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCgg                                             >  1:1026266/1‑100 (MQ=255)
                                               cTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCgg                                             <  2:698626/100‑1 (MQ=255)
                                                 ggCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCg                                           <  2:456561/100‑1 (MQ=255)
                                                 ggCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCg                                           <  1:730188/100‑1 (MQ=255)
                                                 ggCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCg                                           >  1:2059063/1‑100 (MQ=255)
                                                 ggCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTCCGGGACCGTCGCGATTTTCCTCTACCCcgcg                                           >  2:735651/1‑100 (MQ=255)
                                                   cGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGat                                         <  1:550260/100‑1 (MQ=255)
                                                   cGACTGAGCCCGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTCCGTGACCGTCGCGATTTTCCTCTACCCGGCGat                                         <  2:1681467/100‑1 (MQ=255)
                                                    gACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGata                                        <  1:1959474/100‑1 (MQ=255)
                                                    gACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGata                                        >  2:1789819/1‑100 (MQ=255)
                                                      cTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGatata                                      >  2:507537/1‑100 (MQ=255)
                                                      cTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGatata                                      >  2:2449388/1‑100 (MQ=255)
                                                      cTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAACTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGatata                                      <  2:1798571/100‑1 (MQ=255)
                                                         aGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATcc                                   >  1:823401/1‑100 (MQ=255)
                                                          gCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCg                                  >  1:2358219/1‑100 (MQ=255)
                                                          gCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCg                                  >  1:1668039/1‑100 (MQ=255)
                                                           cggtgaGTTAAGGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGGTATATTCGGGACCTTCGCGATTTTCCTTTACCCGGCGATATATCCg                                  <  2:107371/94‑1 (MQ=255)
                                                           ccGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGc                                 <  1:2449402/100‑1 (MQ=255)
                                                           ccGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCACGATTTTCCTCTACCCGGCGATATATCCGc                                 <  1:1410179/100‑1 (MQ=255)
                                                              gTAGTGAAAGCGGAAGCCAGTAAAGAAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGa                              <  2:2107560/100‑1 (MQ=255)
                                                               tAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCt                                <  2:447965/97‑1 (MQ=255)
                                                                 gTGAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGt                           <  2:1528132/100‑1 (MQ=255)
                                                                   gAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtct                         >  2:594072/1‑100 (MQ=255)
                                                                   gAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtct                         >  2:1031142/1‑100 (MQ=255)
                                                                   gAAAGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtct                         >  2:2207514/1‑100 (MQ=255)
                                                                    aaaGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtctc                        >  2:134349/1‑100 (MQ=255)
                                                                    aaaGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGtctc                        <  2:2512760/100‑1 (MQ=255)
                                                                      aGCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCaa                      >  2:1547890/1‑100 (MQ=255)
                                                                       gCGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAAt                     >  2:869906/1‑100 (MQ=255)
                                                                        cGGAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATg                    >  2:1332180/1‑100 (MQ=255)
                                                                         gaaaGCCAGTACATTAACCGAAGCTGTTGCAACCGTAGTTATCGTCAGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATgg                   <  2:1517/98‑1 (MQ=255)
                                                                          gaatcgcggaaagtgaTCGGGGCTGTTTCAACGCTTGTTATTTTCGGGACCGGCGCGGTTTTCCTCTGCCCGGCGATATATCCGCTGATGTCTCACTgg                   <  2:1503/90‑1 (MQ=255)
                                                                          gAAGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGt                  >  2:123553/1‑100 (MQ=255)
                                                                           aaGCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGtt                 <  2:1664738/100‑1 (MQ=255)
                                                                             gCCAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGttt                <  1:1833658/99‑1 (MQ=255)
                                                                              ccAGTAAAGTAACCGAGGCTGTTGCAACCGTTGTTAGCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAg              <  2:1413094/100‑1 (MQ=255)
                                                                                 gtaaagtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTcc           >  1:1257174/1‑100 (MQ=255)
                                                                                 gtaaagtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTcc           <  1:1357824/100‑1 (MQ=255)
                                                                                   aaagtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCgg         <  1:875040/100‑1 (MQ=255)
                                                                                   aaagtaaCCGAGGCTGTTGCAACCGTTGTTATATTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATCTCCGCTGATGTCTCAATGGTTTAGTc            >  2:478442/1‑97 (MQ=255)
                                                                                    aagtaaCCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCAGGa        >  1:2487363/1‑100 (MQ=255)
                                                                                    aagtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGTACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGa        <  2:2547394/100‑1 (MQ=255)
                                                                                      gtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGaaa      >  1:2469275/1‑100 (MQ=255)
                                                                                      gtaaCCGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGaaa      <  1:869914/100‑1 (MQ=255)
                                                                                          ccGAGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACttt  >  2:104796/1‑100 (MQ=255)
                                                                                             |                                                                                                
GCTGGTTGATCGGTGCCGGTAGCAGTATCTGTGGTGCTGCCGCGGTACTGGCGACTGAGCCGGTAGTGAAAGCGGAAGCCAGTAAAGTAACCGTGGCTGTTGCAACCGTTGTTATCTTCGGGACCGTCGCGATTTTCCTCTACCCGGCGATATATCCGCTGATGTCTCAATGGTTTAGTCCGGAAACTTT  >  NC_000913/2250117‑2250306

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: