Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,723,474 G→A 6.4% A288A (GCG→GCA)  nepI_1 → Purine ribonucleoside efflux pump NepI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,723,4740GA 6.4% 189.0 / 3.3 62A288A (GCG→GCA) nepI_1Purine ribonucleoside efflux pump NepI
Reads supporting (aligned to +/- strand):  ref base G (23/35);  new base A (2/2);  total (25/37)
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.18e-01

CTGGCATTGGGGCAGGCATCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGCCTGGCTGCT  >  KT2440/2723343‑2723608
                                                                                                                                   |                                                                                                                                      
cTGGCATTGGGGCAGGCATCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGCGATGTGGCGCAg                                                                                                                                     <  4:271114/135‑1 (MQ=255)
   gCATTGGGGCAGGCATCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCg                                                                                                                                  <  8:81398/135‑1 (MQ=255)
    cATTGGGGCAGGCATCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGa                                                                                                                                 <  6:191680/135‑1 (MQ=255)
          ggcaggcaTCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           <  7:199583/135‑1 (MQ=255)
          ggcaggcaTCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           <  4:145163/135‑1 (MQ=255)
          ggcaggcaTCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           <  3:108853/135‑1 (MQ=255)
          ggcaggcaGCGTCGTGGCTTGTGTAGCTGAGCATCGCCTTGGGGGTAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           <  5:372966/135‑1 (MQ=255)
            caggcaTCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTg                                                                                                                         <  6:98601/135‑1 (MQ=255)
                     tcgtTGCTTGTGTACCTTAGCATCGCCTTGTGGGGAATGACCTATGTCGTCTCCCCGACCCTGCTGCAGACAGCGTTCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGtgtg                                                                                                                 <  4:114596/134‑1 (MQ=255)
                        tggggtgtgTTCCTGAGCATCGCCTTGTGGGGCATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGGGGGGATGTGGCGCAGTCGCTGCTGGTGGCGGGGTggc                                                                                                               >  8:26143/6‑132 (MQ=255)
                        tGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGaa                                                                                                              <  4:244832/134‑1 (MQ=255)
                            ttGTGTACCTGAGCATCGCCTTGTGGGGAATGAGCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         <  1:320810/135‑1 (MQ=255)
                            ttGTGTACCTGAGCATCGCCTTGAGGTGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         <  2:173279/135‑1 (MQ=255)
                                  aCCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         >  7:354642/1‑129 (MQ=255)
                                  aCCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         <  8:354642/129‑1 (MQ=255)
                                  aCCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCAtt                                                                                                    <  4:192706/134‑1 (MQ=255)
                                  aCCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTg                                                                                                   >  5:238590/1‑135 (MQ=255)
                                        gCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTgg                                                                                             <  6:346738/135‑1 (MQ=255)
                                          aTCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGc                                                                                           >  3:212269/1‑135 (MQ=255)
                                               cTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGaa                                                                                       >  1:389792/1‑134 (MQ=255)
                                               cTTGTGGGGAATGACCTATGGCGGCGCCCCGACACTGCTGCAGACACCGTGCCCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTCGAACTGCGCCCTTGCTCTGCGCGGaa                                                                                       >  8:204105/1‑134 (MQ=255)
                                                           gACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGCCGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGAtt                                                                          >  4:176365/1‑135 (MQ=255)
                                                           gACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGAtt                                                                          >  6:320798/1‑135 (MQ=255)
                                                             ccTATGGCGGCGCCCCGACCCGGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGGGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGAttgt                                                                        <  2:366984/135‑1 (MQ=255)
                                                                   gcggcgCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           >  1:34627/1‑78 (MQ=255)
                                                                   gcggcgCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTgg                                                                                                                           <  2:34627/78‑1 (MQ=255)
                                                                   gcggcgCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         >  3:26860/1‑96 (MQ=255)
                                                                   gcggcgCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTgcg                                                                                                         <  4:26860/96‑1 (MQ=255)
                                                                       cgcCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGctggctgg                                                              <  5:320798/135‑1 (MQ=255)
                                                                         ccccGACCCTGCTGCAGACAGCGTGCGCCGCCGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTgggg                                                            >  2:83133/1‑135 (MQ=255)
                                                                          cccGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCAtt                                                                                                    >  8:204306/1‑94 (MQ=255)
                                                                          cccGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCAtt                                                                                                    <  7:204306/94‑1 (MQ=255)
                                                                          cccGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGc                                                           >  5:261658/1‑135 (MQ=255)
                                                                          cccGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCACAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGc                                                           >  7:8455/1‑135 (MQ=255)
                                                                          cccGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCACAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGc                                                           >  7:14411/1‑135 (MQ=255)
                                                                                 ctgctgCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGGGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCTGATTGTTGCTGGCTGGGGCCggcgtg                                                    <  8:59226/135‑1 (MQ=255)
                                                                                      gcAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGgcgtggcgtg                                               >  4:26840/1‑135 (MQ=255)
                                                                                       cAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGcgtggcgtgg                                              >  4:285125/1‑135 (MQ=255)
                                                                                       cAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGAGTGTTGCTGGCTGGGGCCGGcgtggcgtgg                                              >  8:80909/1‑135 (MQ=255)
                                                                                         gACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGtt                                            <  6:261658/135‑1 (MQ=255)
                                                                                         gACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGtt                                            <  5:91645/135‑1 (MQ=255)
                                                                                            aGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTgg                                         <  6:238590/135‑1 (MQ=255)
                                                                                              cGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCg                                       <  5:51798/135‑1 (MQ=255)
                                                                                              cGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCg                                       <  3:264207/135‑1 (MQ=255)
                                                                                              cGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCg                                       <  3:146526/135‑1 (MQ=255)
                                                                                              cGTGCGCCCACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCg                                       <  4:107430/135‑1 (MQ=255)
                                                                                               gTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCgg                                      <  6:313852/135‑1 (MQ=255)
                                                                                               gTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCgg                                      <  2:14108/135‑1 (MQ=255)
                                                                                               gTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCACAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCgg                                      <  8:14411/135‑1 (MQ=255)
                                                                                               gTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCACAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCgg                                      <  8:8455/135‑1 (MQ=255)
                                                                                                         gccgcAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTg                            >  6:90396/1‑135 (MQ=255)
                                                                                                           cgcAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTggg                          >  8:365850/1‑135 (MQ=255)
                                                                                                             cAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGt                        <  1:83133/135‑1 (MQ=255)
                                                                                                                      gtggtgATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTgctggctg                                                               <  6:289468/87‑1 (MQ=255)
                                                                                                                      gtggtgATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTgctggctg                                                               >  5:289468/1‑87 (MQ=255)
                                                                                                                        ggtgATGTGGCGCAGTCGATGCTGGTGACGGTGGGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGc                                                                               >  8:374239/1‑69 (MQ=255)
                                                                                                                        ggtgATGTGGCGCAGTCGATGCTGGGGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGc                                                                               <  7:374239/69‑1 (MQ=255)
                                                                                                                          tgATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGc           <  2:87892/135‑1 (MQ=255)
                                                                                                                           gATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGcc          <  3:26840/135‑1 (MQ=255)
                                                                                                                             tgtgGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGCCTg        >  1:117487/1‑135 (MQ=255)
                                                                                                                             tgtgGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTGTGGGGGGGTTGTGCTGGGGTTGATCGCCTTGGCCTg        >  1:288120/1‑135 (MQ=255)
                                                                                                                             tgtgGCGCAGTCGATGCTGGTGACGGTGTCGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGCCTg        >  8:33043/1‑135 (MQ=255)
                                                                                                                                  cgcAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTTGGGGGGTAGTGCTGGGGTTGCTCGCCATGGCCTGgctgc   >  2:365805/1‑135 (MQ=255)
                                                                                                                                   gcAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGAGTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGAACGCCTTGGCCTGgctgct  <  7:187013/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CTGGCATTGGGGCAGGCATCGTCGTGGCTTGTGTACCTGAGCATCGCCTTGTGGGGAATGACCTATGGCGGCGCCCCGACCCTGCTGCAGACAGCGTGCGCCGACGCCGCAGGCCAGGGTGGTGATGTGGCGCAGTCGATGCTGGTGACGGTGTGGAACTGCGCCATTGCTCTGGGCGGAATCATTGGCGGATTGTTGCTGGCTGGGGCCGGCGTGGCGTGGTTTGGCGGGGTAGTGCTGGGGTTGATCGCCTTGGCCTGGCTGCT  >  KT2440/2723343‑2723608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑