Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,025,676 +C coding (535/570 nt) ygfQ → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,025,6751.C99.3% 414.0 / ‑4.1 151coding (534/570 nt)ygfQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (72/78);  total (72/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

CGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGG  >  NC_000913/3025600‑3025745
                                                                            |                                                                       
cggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTg                                                                          <  2:1298692/76‑1 (MQ=255)
 ggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGc                                                                        <  2:1817445/76‑1 (MQ=255)
 ggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGc                                                                        >  2:4849122/1‑76 (MQ=255)
 ggtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGc                                                                        >  2:5463297/1‑76 (MQ=255)
  gtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGcc                                                                       <  2:3564610/76‑2 (MQ=255)
  gtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGcc                                                                       <  1:5639088/76‑2 (MQ=255)
  gtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGcc                                                                       >  1:720566/1‑75 (MQ=255)
  gtatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGcc                                                                       <  2:2138694/76‑2 (MQ=255)
   tatcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGcca                                                                      >  1:445257/1‑74 (MQ=255)
    atcggGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGccag                                                                     >  1:4727704/1‑73 (MQ=255)
       gggCTGTTTCTGCTGCTGATTGCTGCTAGCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTgg                                                                  <  2:2505778/76‑1 (MQ=255)
       gggCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTgg                                                                  <  1:4229455/76‑1 (MQ=255)
        ggCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGgc                                                                 <  2:453941/76‑1 (MQ=255)
            gTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTc                                                             <  2:2876728/76‑1 (MQ=255)
              ttCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCgg                                                           >  1:2883321/1‑76 (MQ=255)
               tCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGt                                                          <  2:1120065/76‑1 (MQ=255)
                ctgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg                                                         <  1:1479896/76‑1 (MQ=255)
                ctgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg                                                         >  1:2311457/1‑76 (MQ=255)
                ctgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTg                                                         <  2:1535113/76‑1 (MQ=255)
                 tgctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGc                                                        <  2:2339230/76‑1 (MQ=255)
                   ctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGt                                                      >  2:5015325/1‑76 (MQ=255)
                   ctgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGt                                                      >  2:2190978/1‑76 (MQ=255)
                    tgctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCgnn                                                     <  1:1044100/76‑3 (MQ=255)
                     gctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGttt                                                    <  2:2598095/76‑1 (MQ=255)
                     gctgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGttt                                                    >  1:5397718/1‑76 (MQ=255)
                       tgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAc                                                  >  1:2763001/1‑76 (MQ=255)
                       tgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAc                                                  >  1:3577192/1‑76 (MQ=255)
                       tgATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAc                                                  >  2:3866331/1‑76 (MQ=255)
                        gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcc                                                 >  1:3375161/1‑76 (MQ=255)
                        gATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcc                                                 >  2:1648728/1‑76 (MQ=255)
                         gtggctggtAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                <  2:1014019/73‑1 (MQ=255)
                         aTTGCTGCTATCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                <  2:1158910/76‑1 (MQ=255)
                         aTTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                >  2:1783018/1‑76 (MQ=255)
                         aTTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcct                                                >  1:4493857/1‑76 (MQ=255)
                          ttGCTGCTAACGGTGTGGTTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctc                                               <  1:1507133/76‑1 (MQ=255)
                           tgctgcTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcc                                              >  2:1960073/1‑76 (MQ=255)
                           tgctgcTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcc                                              >  2:4675513/1‑76 (MQ=255)
                            gctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcct                                             <  2:3708549/76‑1 (MQ=255)
                            gctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcct                                             <  1:4075831/76‑1 (MQ=255)
                            gctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcct                                             >  1:2095196/1‑76 (MQ=255)
                            gctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTAcctcct                                             >  1:5514497/1‑76 (MQ=255)
                             ctgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCtt                                            <  2:3357099/76‑1 (MQ=255)
                              tgctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTc                                           <  1:1960073/76‑1 (MQ=255)
                               gctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTGCCTCCTTcc                                          <  1:5800627/76‑1 (MQ=255)
                               gctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTcc                                          >  2:4422659/1‑76 (MQ=255)
                                ctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTccc                                         <  1:1847184/76‑1 (MQ=255)
                                ctAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTccc                                         >  2:583447/1‑76 (MQ=255)
                                ctAACGGTGTGGGCATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTccc                                         >  1:1121987/1‑76 (MQ=255)
                                  aaCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCgg                                       <  2:2041822/76‑1 (MQ=255)
                                   aCGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGt                                      >  1:1151833/1‑76 (MQ=255)
                                     ggTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtga                                    >  1:4178222/1‑76 (MQ=255)
                                     ggTGTGGGTATCGTTATCACAAACCCGATTGAGGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGTtga                                    <  1:5207856/76‑1 (MQ=255)
                                      gtgtGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgat                                   <  1:557255/76‑1 (MQ=255)
                                       tgtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatg                                  <  1:1648728/76‑1 (MQ=255)
                                       tgtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatg                                  >  1:3631815/1‑76 (MQ=255)
                                       tgtgGGGATGGTTATCAAAAACCCGAGTGAAGGCTTGCCAGTGGAGCGCGGTGCGTTTACCTCCTTCCCGGcgatg                                  >  1:3979750/1‑76 (MQ=255)
                                        gtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatga                                 >  1:1314508/1‑76 (MQ=255)
                                        gtgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatga                                 <  1:3804369/76‑1 (MQ=255)
                                        gtgGGCACGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGCggatga                                 >  2:2851310/1‑76 (MQ=255)
                                         tgGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgat                                >  2:5125315/1‑76 (MQ=255)
                                         tgGGGATGGTTATAAAAAACCAGATTGAAGGCTTGCCAGTGGCGCTCGGGGCGTTTACCTCCATACCGGcgctgat                                >  1:5716552/1‑76 (MQ=255)
                                          gtgTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgatg                               <  1:5851142/74‑1 (MQ=255)
                                          gggTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgatg                               <  1:876967/76‑1 (MQ=255)
                                           ggTAGGGTTATCAAAAGCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGtgatgatga                              <  2:5231606/76‑1 (MQ=255)
                                            gTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAg                             <  1:2963151/76‑1 (MQ=255)
                                            gTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAg                             >  1:3487518/1‑76 (MQ=255)
                                            gTATGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAg                             <  1:5439569/76‑1 (MQ=255)
                                             tATGGTTATCAAACACCCGATTGAAGGCTTGCGAGTGGCGCTCGGTGCTTTTACCTCCTTCCCGGTGATGATGAGc                            <  2:2675747/76‑1 (MQ=255)
                                              aTGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCt                           >  1:5761400/1‑76 (MQ=255)
                                              aTGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCt                           >  1:4718013/1‑76 (MQ=255)
                                               tGGTTATCGAAAGCCCGCTTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                          <  2:1089849/76‑1 (MQ=255)
                                               tGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                          <  2:1490778/76‑1 (MQ=255)
                                               tGGTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCtt                          <  2:4779191/76‑1 (MQ=255)
                                                ggTTATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTg                         >  2:1305771/1‑76 (MQ=255)
                                                agTTATTAAAAACCCGCTTGAAGGCATGCCAGTGGCGCTCGGAGCCTTTACATCCTGCCCGGCGTTGATGGGCtgg                         >  2:2919742/2‑74 (MQ=255)
                                                   tATCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTg                      >  1:5086105/1‑76 (MQ=255)
                                                    aTCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTgg                     >  2:3804966/1‑76 (MQ=255)
                                                    aTCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTgg                     >  2:518660/1‑76 (MQ=255)
                                                    aTCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTgg                     >  1:4831154/1‑76 (MQ=255)
                                                     tCAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTggc                    >  2:5387621/1‑75 (MQ=255)
                                                      taaaatCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTGCCCGGTGATGATGAGCTTGCTgggg                   <  2:4262101/75‑1 (MQ=255)
                                                      cAAAAACCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTgggg                   <  1:3706152/76‑1 (MQ=255)
                                                      cAAAAACCCGATGGAAGGCTTGCCAGTGGCGCGCGGTGCGTTTACCTCCTTCACGGTGAGGATGAGCTTTCTgggc                   >  2:1880162/1‑75 (MQ=255)
                                                       aaaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGc                  >  1:1564650/1‑76 (MQ=255)
                                                        aaaaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGATCTTGCTGGGGCt                 <  1:3506406/76‑1 (MQ=255)
                                                          aaCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTgg               <  1:4405905/76‑1 (MQ=255)
                                                           gcccGTTTGAATGCTCGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGc              <  1:2851310/75‑1 (MQ=255)
                                                           aCCCGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTggg              >  1:80448/1‑75 (MQ=255)
                                                             ccGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCgg            >  1:1049844/1‑76 (MQ=255)
                                                             ccGATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCgg            <  1:2431521/76‑1 (MQ=255)
                                                               gATTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCgggc          >  1:4474709/1‑74 (MQ=255)
                                                                aTTGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTcc         >  1:1751761/1‑75 (MQ=255)
                                                                 ttGAAGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCAt        >  1:5022928/1‑76 (MQ=255)
                                                                    aaGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCtt     <  1:2393391/76‑1 (MQ=255)
                                                                    aaGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGGCATcct     >  1:722123/1‑74 (MQ=255)
                                                                     aGGCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGGCCTCTTc    >  2:2126604/1‑76 (MQ=255)
                                                                      ggCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCg   <  1:2250839/76‑1 (MQ=255)
                                                                      ggCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCg   <  2:5686362/76‑1 (MQ=255)
                                                                      ggCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGGCATCTTCg   >  1:2995617/1‑76 (MQ=255)
                                                                       gCTTGCCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCgg  <  1:4631568/76‑1 (MQ=255)
                                                                            |                                                                       
CGGTATCGGGCTGTTTCTGCTGCTGATTGCTGCTAACGGTGTGGGTATGGTTATCAAAAACCCGATTGAAGGCTTGCAGTGGCGCTCGGTGCGTTTACCTCCTTCCCGGTGATGATGAGCTTGCTGGGGCTGGCGGTCATCTTCGG  >  NC_000913/3025600‑3025745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: