Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,438 C→T 6.8% R318H (CGC→CAC)  PP_0928 ← K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,4380CT 6.8% 168.1 / 2.8 59R318H (CGC→CAC) PP_0928K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGAGCATTACA  >  NC_002947/1070290‑1070587
                                                                                                                                                    |                                                                                                                                                     
gATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAggcg                                                                                                                                                      >  1:1039130/1‑150 (MQ=255)
  tGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAggcggc                                                                                                                                                    >  1:913507/1‑150 (MQ=255)
        ccTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgccagcca                                                                                                                                              >  2:213107/1‑150 (MQ=255)
          tAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAAc                                                                                                                                            >  1:125659/1‑150 (MQ=255)
                   gCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATAccc                                                                                                                                   <  1:180529/150‑1 (MQ=255)
                     gTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCAt                                                                                                                                 <  2:1402036/150‑1 (MQ=255)
                     gTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCAt                                                                                                                                 <  2:1305868/150‑1 (MQ=255)
                                cTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAgccg                                                                                                                      <  1:213107/150‑1 (MQ=255)
                                     gcgTTCCCCTTTCATGCCCCTCACTGCGGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCGCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCcaca                                                                                                                 <  2:1039130/150‑1 (MQ=255)
                                     gcgTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCcaca                                                                                                                 >  1:1180365/1‑150 (MQ=255)
                                          ccccTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATAcc                                                                                                                                    <  1:612047/126‑1 (MQ=255)
                                          ccccTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATAcc                                                                                                                                    >  2:612047/1‑126 (MQ=255)
                                          ccccTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAg                                                                                                            <  2:1180365/150‑1 (MQ=255)
                                           cccTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAgg                                                                                                           >  2:1413384/1‑150 (MQ=255)
                                             cTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCCTGCTAAAAGCCGCCACATGAAGGcc                                                                                                         >  2:442795/1‑150 (MQ=255)
                                                   tGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGc                                                                                                   >  1:1039797/1‑150 (MQ=255)
                                                        cTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGt                                                                                              >  1:1207780/1‑150 (MQ=255)
                                                        cTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAAACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAACCAACGGCCTACCCCTGGTAAAAGCCGCCACATGAAGGCCGTAGGCAAGGt                                                                                              >  2:114214/1‑150 (MQ=255)
                                                            aTGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAg                                                                                          <  1:1413384/150‑1 (MQ=255)
                                                                  ccTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGTAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAg                                                                                    >  2:310485/1‑150 (MQ=255)
                                                                  ccTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGTGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAg                                                                                    >  2:445649/1‑150 (MQ=255)
                                                                  ccTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGGCACGTTGAGGcccga                                                                                    >  2:1371325/1‑146 (MQ=255)
                                                                           aGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAg                                                                           <  2:651980/150‑1 (MQ=255)
                                                                            gcacgcacAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAgg                                                                          >  2:405723/1‑150 (MQ=255)
                                                                                          gtACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACgg                                                                                                                                          >  2:1056374/1‑72 (MQ=255)
                                                                                          gtACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACgg                                                                                                                                          <  1:1056374/72‑1 (MQ=255)
                                                                                          gtACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACAt                                                                                                                >  1:147501/1‑98 (MQ=255)
                                                                                          gtACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACAt                                                                                                                <  2:147501/98‑1 (MQ=255)
                                                                                            aCAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGa                                                          <  1:405723/150‑1 (MQ=255)
                                                                                              aGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATg                                                        >  2:898291/1‑150 (MQ=255)
                                                                                               gCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGc                                                       <  2:902253/150‑1 (MQ=255)
                                                                                                cAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGTGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCg                                                      <  1:445649/150‑1 (MQ=255)
                                                                                                cAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACg                                                                                                                                           >  1:225320/1‑65 (MQ=255)
                                                                                                cAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACg                                                                                                                                           <  2:225320/65‑1 (MQ=255)
                                                                                                cAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCg                                                      >  2:533197/1‑150 (MQ=255)
                                                                                                 aaCACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCaaa                                                                                     <  1:834083/118‑1 (MQ=255)
                                                                                                 aaCACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCaaa                                                                                     >  2:834083/1‑118 (MQ=255)
                                                                                                        tAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCgg                                              >  1:1012521/1‑150 (MQ=255)
                                                                                                        tAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCgg                                              >  2:388597/1‑150 (MQ=255)
                                                                                                        tAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCg                                               >  2:1076401/1‑149 (MQ=255)
                                                                                                        tAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCg                                               <  1:1076401/149‑1 (MQ=255)
                                                                                                           cGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTGATGAATTTACAGGCAAAGAAACAGCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAAt                                           >  1:1085634/1‑150 (MQ=255)
                                                                                                           cGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGACGTAGGACAGGTGCAGGCAAAGAAACACCACACCTTCTGCGCGAGTGATGCGGTAACCGGAAt                                           <  1:388597/150‑1 (MQ=255)
                                                                                                                 ccGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAg                                                             <  2:494990/126‑1 (MQ=255)
                                                                                                                 ccGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAg                                                             >  1:494990/1‑126 (MQ=255)
                                                                                                                  cGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAgccgc                                                                                                                     <  1:381535/69‑1 (MQ=255)
                                                                                                                  cGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAgccgc                                                                                                                     >  2:381535/1‑69 (MQ=255)
                                                                                                                                  aGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAacac                                                                              <  1:1080864/92‑1 (MQ=255)
                                                                                                                                  aGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAacac                                                                              >  2:1080864/1‑92 (MQ=255)
                                                                                                                                   gCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGgcag                   <  2:1207780/150‑1 (MQ=255)
                                                                                                                                        aaTCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGaa                                       <  1:840302/125‑1 (MQ=255)
                                                                                                                                        aaTCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGaa                                       >  2:840302/1‑125 (MQ=255)
                                                                                                                                        aaTCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACg              >  2:1186583/1‑150 (MQ=255)
                                                                                                                                         aTCGCTCCAGGTGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGc                                >  2:1043134/1‑131 (MQ=255)
                                                                                                                                         aTCGCTCCAGGTGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGc                                <  1:1043134/131‑1 (MQ=255)
                                                                                                                                         aTCGCTCCAGGCGGCCAGGCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACCCCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGATCACGGGCAGGCGCAGGGCAGCCACGc             <  2:382347/150‑1 (MQ=255)
                                                                                                                                            gCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGa          >  1:276773/1‑150 (MQ=255)
                                                                                                                                             cgccAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGAg         <  2:1056667/148‑1 (MQ=255)
                                                                                                                                              tCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAg                        >  1:496881/1‑134 (MQ=255)
                                                                                                                                              tCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAg                        <  2:496881/134‑1 (MQ=255)
                                                                                                                                                cAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGAGCAt      >  2:133876/1‑150 (MQ=255)
                                                                                                                                                    cggcCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGAGCATTACa  <  2:913507/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACCGGAATAGAACACTGGCAGGCTCAGGGCAGCAACGCCGAGCATTACA  >  NC_002947/1070290‑1070587

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