Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,937,436 G→C 6.6% noncoding (61/90 nt) tRNA Ser anticodon TGA_ Cove score 61.04 ← tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,937,4360GC 6.6% 314.9 / 6.1 81noncoding (61/90 nt)tRNA Ser anticodon TGA_ Cove score 61.04tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base C (3/3);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACCCGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGAACCTGCGA  >  BX571857/1937340‑1937532
                                                                                                |                                                                                                
cGAACCCGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAg                                                                                                  <  1:860665/97‑1 (MQ=255)
  aaCCCGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGcc                                                                                                <  1:746257/97‑1 (MQ=255)
    cccGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGt                                                                                              >  2:1097889/1‑97 (MQ=255)
    cccGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGt                                                                                              >  1:745514/1‑97 (MQ=255)
     ccGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGtt                                                                                             <  2:1051762/97‑1 (MQ=255)
     ccGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGtt                                                                                             >  1:659072/1‑97 (MQ=255)
      cgcgATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTa                                                                                            >  1:543160/1‑97 (MQ=255)
      cgcgATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTa                                                                                            <  1:559331/97‑1 (MQ=255)
      cgcgATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTa                                                                                            <  2:177780/97‑1 (MQ=255)
        cgATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCCTTAGTAAAACGGAGGAAGCGGGATTCGAACCCCCGCGCGCCGTTCAg                                                                                          >  2:378947/1‑97 (MQ=25)
          aTCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGcc                                                                                        <  1:873235/97‑1 (MQ=255)
           tctcCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGccc                                                                                       >  1:1481254/1‑97 (MQ=255)
            ctcCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGcccc                                                                                      >  2:90316/1‑97 (MQ=255)
            ctcCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGcccc                                                                                      >  1:648095/1‑97 (MQ=255)
             tcCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCt                                                                                     >  2:455365/1‑97 (MQ=255)
                    tGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGtt                                                                              <  2:1047938/97‑1 (MQ=255)
                     gACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGttt                                                                             <  2:1186786/97‑1 (MQ=255)
                     gACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGttt                                                                             >  1:207110/1‑97 (MQ=255)
                       caggcaggCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTc                                                                           >  1:720923/1‑97 (MQ=255)
                       caggcaggCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTc                                                                           >  1:512279/1‑97 (MQ=255)
                        aggcaggcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCa                                                                          >  1:557453/1‑97 (MQ=255)
                         ggcaggcGTGCTGACCGCTCCCCCACGAGACCATTAGTAAGGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCaa                                                                         <  2:1218432/97‑1 (MQ=17)
                          gcaggcGTGTTAACCGCTCCACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAg                                                                        <  2:1565663/97‑1 (MQ=255)
                          gcaggcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAg                                                                        <  1:807367/97‑1 (MQ=255)
                            aggcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGAc                                                                      >  1:981377/1‑97 (MQ=255)
                             ggcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGAcc                                                                     <  1:619974/97‑1 (MQ=255)
                             ggcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGAcc                                                                     <  1:245188/97‑1 (MQ=255)
                              gcGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCg                                                                    <  2:136362/97‑1 (MQ=255)
                               cGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGa                                                                   <  2:828177/97‑1 (MQ=255)
                               cGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGa                                                                   >  1:422260/1‑97 (MQ=255)
                                  gtTAACCGCTACACTGCGAGACCAATAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATcc                                                                <  1:1293443/97‑1 (MQ=255)
                                   ttAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATccc                                                               >  1:96330/1‑97 (MQ=255)
                                   ttAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATccc                                                               >  1:1323961/1‑97 (MQ=255)
                                   ttAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATccc                                                               >  1:455660/1‑97 (MQ=255)
                                      aCCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTc                                                            <  1:438074/97‑1 (MQ=255)
                                      aCCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTc                                                            <  1:438094/97‑1 (MQ=255)
                                       ccGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCa                                                           <  1:921617/97‑1 (MQ=255)
                                        cGCTACACTACGAGACCATTAGTAAAGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAg                                                          <  1:1497195/97‑1 (MQ=255)
                                        cGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGGTTTCAAGACCGATCCCTTCAg                                                          >  2:207488/1‑97 (MQ=255)
                                            acacTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCgg                                                      >  2:213920/1‑97 (MQ=255)
                                                     gaCCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGGCGGTGTTCAAGACCGATCCCCTCAGCCGGACTAgggca                                             >  2:1118577/1‑95 (MQ=14)
                                                      aCCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGCCGGTTTTCAAGACCGATCCCTTCCGCCGGCCTTGGGTac                                            >  2:1315679/1‑96 (MQ=14)
                                                       ccATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTCTTCAAGACCGATCCCTTCAGCCGGACTTGGGTAtt                                           >  2:1351004/1‑97 (MQ=255)
                                                        cATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtc                                          >  1:1203137/1‑97 (MQ=255)
                                                         aTTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtcc                                         >  2:1409012/1‑97 (MQ=255)
                                                          ttAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtcct                                        >  1:781793/1‑97 (MQ=40)
                                                           tAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtcctc                                       >  2:148685/1‑97 (MQ=39)
                                                              tAAAGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCaa                                    <  1:1283856/97‑1 (MQ=18)
                                                                aaaCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCaaaa                                  >  2:1467355/1‑97 (MQ=38)
                                                                 ggcGGAGGAAGAGGGATTCGAACCCCCGCGGGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAAt                                 <  2:1287474/95‑1 (MQ=16)
                                                                  aCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAAtt                                >  2:986027/1‑97 (MQ=37)
                                                                  aCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAAtt                                <  2:1333436/97‑1 (MQ=38)
                                                                     gaggaAGAGGGATTCGAACCCCCGCGGCCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATtata                             <  2:477186/97‑1 (MQ=255)
                                                                     gaggaAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATtata                             >  2:905748/1‑97 (MQ=40)
                                                                      aggaAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTatat                            <  2:113571/97‑1 (MQ=255)
                                                                       ggaAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATg                           <  1:611665/97‑1 (MQ=255)
                                                                       gcaggagGGATTCGAACCCCCGCGACCCGTTAAGCCCCTGTCGGTTTTCAAGACCGAGCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATg                           <  1:877804/93‑1 (MQ=14)
                                                                         aaGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGa                         <  1:137903/97‑1 (MQ=255)
                                                                          agagGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGaa                        >  2:746388/1‑96 (MQ=11)
                                                                          agagGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGAc                        >  1:801179/1‑97 (MQ=255)
                                                                           gagGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGAcc                       >  1:127987/1‑97 (MQ=255)
                                                                           gagGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGAcc                       >  1:918635/1‑97 (MQ=255)
                                                                           gagGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGAcc                       <  1:480466/97‑1 (MQ=255)
                                                                            agGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCt                      >  2:790677/1‑97 (MQ=255)
                                                                             gggATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCtt                     >  1:906838/1‑97 (MQ=11)
                                                                             gggATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCtt                     >  2:468896/1‑97 (MQ=255)
                                                                               gATTCGAACCCCCGCGGCCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTCGc                   >  2:878150/1‑97 (MQ=255)
                                                                               gATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGc                   >  2:685136/1‑97 (MQ=255)
                                                                                  tCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAgg                <  2:1275854/97‑1 (MQ=255)
                                                                                   cGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGa               >  2:111033/1‑97 (MQ=255)
                                                                                      aCCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTc            >  1:963105/1‑97 (MQ=255)
                                                                                       cccccGCGCGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCg           <  2:1293953/97‑1 (MQ=255)
                                                                                       cccccGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGGCTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCg           >  2:1290672/1‑97 (MQ=255)
                                                                                       cccccGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCg           <  1:69231/97‑1 (MQ=255)
                                                                                       cccccGCGAGCCGTTAAGCCCCTGTCGGCTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCg           >  1:1039674/1‑97 (MQ=255)
                                                                                        ccccGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGa          >  2:194338/1‑97 (MQ=255)
                                                                                        ccccGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGa          >  1:364077/1‑97 (MQ=255)
                                                                                         cccGCGGGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGaa         <  2:1498760/97‑1 (MQ=255)
                                                                                         cccGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGaa         <  1:1206863/97‑1 (MQ=255)
                                                                                           cgcgAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGCTCCCTTCAGCCGGACTTGGGTATTCCTCCACAATTATATGGACCTTGCAGGACTCGAAcc       >  1:1405542/1‑97 (MQ=255)
                                                                                               aGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGAACCTGCg   <  1:960511/97‑1 (MQ=255)
                                                                                               aGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGAACCTGCg   >  1:1246730/1‑97 (MQ=255)
                                                                                                gCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGAACCTGCGa  >  1:909656/1‑97 (MQ=255)
                                                                                                |                                                                                                
CGAACCCGCGATCTCCTGCGTGACAGGCAGGCGTGTTAACCGCTACACTACGAGACCATTAGTAAAACGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATTATATGGACCTTGCAGGACTCGAACCTGCGA  >  BX571857/1937340‑1937532

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: