Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,799 C→G 9.0% T517T (ACG→ACC tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,7990CG 9.0% 139.6 / 7.9 56T517T (ACG→ACCtynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base G (3/2);  total (31/25)
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.04e-01

CTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTA  >  NC_000913/1449651‑1449946
                                                                                                                                                    |                                                                                                                                                   
cTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGccc                                                                                                                                                     <  1:348146/149‑1 (MQ=255)
    aTCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAg                                                                                                                                                 >  2:952682/1‑149 (MQ=255)
         gATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                            <  2:869341/149‑1 (MQ=255)
             gAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTGATCTTTCGCCGTCTCATCGTGCAGGGTTTTCGCTTTTACACCTTTCACCGCTTCGATGCCCGTAGAACCGGCAt                                                                                                                                        >  1:583391/1‑149 (MQ=255)
               ggCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCg                                                                                                                                      <  1:597784/149‑1 (MQ=255)
                gCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGa                                                                                                                                     <  1:246583/149‑1 (MQ=255)
                         ataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCCTCTTTCGCCGGCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                            >  2:977535/1‑149 (MQ=255)
                         ataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTAGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                            >  1:255927/1‑149 (MQ=255)
                            aaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGTTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGAATCGATGCCAATAGTGc                                                                                                                         >  1:322107/1‑149 (MQ=255)
                                tgtgtgtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCAGGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGtt                                                                                                                     <  1:709572/144‑1 (MQ=255)
                                   tGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTCGCGCGGGTCATCTTTCGCCGTCTCATCGTGCTTGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTc                                                                                                                  <  1:493177/149‑1 (MQ=255)
                                     gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGACGTTTTCAt                                                                                                                >  2:37243/1‑149 (MQ=255)
                                        tgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGa                                                                                                             >  2:773969/1‑149 (MQ=255)
                                            gtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGAt                                                                                                         >  2:89891/1‑149 (MQ=255)
                                               cccACGATATTGTTATCGATACGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGATTTCGCTTTAACAACTTTCACCGCTTCGATTCCCGTAGAACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                      <  1:836596/149‑1 (MQ=255)
                                               cccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                      >  2:73499/1‑149 (MQ=255)
                                               cccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTATAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                      >  2:876986/1‑149 (MQ=255)
                                                ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAg                                                                                                     <  2:990656/149‑1 (MQ=255)
                                                 cACGATTTTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCCTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGt                                                                                                    <  1:361611/149‑1 (MQ=255)
                                                 cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTACCACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGt                                                                                                    >  1:987179/1‑149 (MQ=255)
                                                 cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCTTCGATGCCAATAGGGCCGTTTTCATGGAAGATCCAGt                                                                                                    >  2:116923/1‑149 (MQ=255)
                                                  aCGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGACATCTTTCGCCGTCGCAGCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCCCCGGCATCGATGCCAATAGTGACGTTTTAATGGAAGATcaattc                                                                                                   >  1:979160/1‑144 (MQ=255)
                                                            gATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAg                                                                                         <  1:952213/149‑1 (MQ=255)
                                                             aTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGt                                                                                        >  1:710428/1‑149 (MQ=255)
                                                              tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCCTCGTGCATGGTTTTCGCTTTAACACCTTTCACCGGTTCGATGCCCGTAGCACCGGCATCGATGGCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                       <  1:936917/149‑1 (MQ=255)
                                                                           cGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAcg                                                                          >  1:141272/1‑148 (MQ=255)
                                                                                cgcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                     <  1:119153/149‑1 (MQ=255)
                                                                                  cgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGa                                                                   <  2:88696/149‑1 (MQ=255)
                                                                                      tCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGGCAAAAATGTAGTCACAGTTACCCACTGTACTGATccc                                                               >  2:234359/1‑148 (MQ=255)
                                                                                       cATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAg                                                              <  1:73499/149‑1 (MQ=255)
                                                                                        aTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgc                                                             <  1:338906/149‑1 (MQ=255)
                                                                                        aTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCAAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAGTGATCCAgc                                                             >  2:349683/1‑149 (MQ=255)
                                                                                        aTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGAAGCCCGTAGCACCGGCAGCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGGAGTCATAGTTACCCACTGTACTGATCCAgc                                                             <  1:206367/149‑1 (MQ=255)
                                                                                        aTCTTTCGCCGTCTCATCGAGCATGGTTTTCGCTTTAACACCTTACACCGCTTCTATGCCCGTAGCCCCGGCATCGATGCCAATAGTGACGTTGTCATGGAAGATACAGTCTAAAATGTAGTCATAGTTACCCACAGATCTGATCCAgc                                                             >  1:294715/1‑149 (MQ=255)
                                                                                          cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGCTGCCAATAGTGCGGTTTTCATGGAAGATCCAGTCAAAAATGTAGACATAGTTACCCACTGTACTGATCCAgcgc                                                           >  2:713444/1‑149 (MQ=255)
                                                                                          cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgcgc                                                           <  1:37243/149‑1 (MQ=255)
                                                                                           tttCGCCGTCTCATCGTGCAGGGATATCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                          <  2:853894/149‑1 (MQ=255)
                                                                                            ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcac                                                         <  2:652368/149‑1 (MQ=255)
                                                                                              cGCCGTATCATCGTTCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacca                                                       <  1:321789/149‑1 (MQ=255)
                                                                                                             cATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCGCGTAGCACCGGCATCGATGCCAATTGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAATGCTCCAGCGCACCACTAACTCCCGGCGtt                                        >  1:280521/1‑149 (MQ=255)
                                                                                                               tGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCg                                      >  2:182393/1‑149 (MQ=255)
                                                                                                               tGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCg                                      <  1:471768/149‑1 (MQ=255)
                                                                                                                  ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGCTGCCAATAGTGCCGTTTTCAAGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACAGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTa                                   >  2:322295/1‑149 (MQ=255)
                                                                                                                     tCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTg                                >  2:797578/1‑149 (MQ=255)
                                                                                                                       gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                              <  1:952682/149‑1 (MQ=255)
                                                                                                                       gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCAGTTTTCATGGAAGATCCAGTCAAAATTGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACGAACTCCCGGCGGTCGGTACTGAc                              >  2:425675/1‑149 (MQ=255)
                                                                                                                            aaCACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTgg                         >  2:659789/1‑149 (MQ=255)
                                                                                                                                 cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTgg                    >  2:449415/1‑149 (MQ=255)
                                                                                                                                  tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTGCAGGCCTTCGGTACTGACTTTGGGCTGGc                   >  2:200231/1‑149 (MQ=255)
                                                                                                                                  tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                              >  2:872038/1‑138 (MQ=255)
                                                                                                                                  tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                              <  1:872038/138‑1 (MQ=255)
                                                                                                                                          cTTCCATGCCCGTAGCACCTGCATCGATGCCAATCTTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTATTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTcc           <  1:208644/149‑1 (MQ=255)
                                                                                                                                           atCGATGCCGGTGCTACGGGCATCGCTGCCAATAGTGGCTTTTTCATGGAAGATCCAGTCAAAAATGTAGACAGAGTACCCCACTGTACTTAACCAGCGCACCACAACCACCCGGCCGTCCCTACAGACTATGGGCTGGCCCATTTCCt          >  1:569575/2‑149 (MQ=255)
                                                                                                                                           atCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          >  1:895513/2‑149 (MQ=255)
                                                                                                                                           atCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt          >  2:394626/2‑149 (MQ=255)
                                                                                                                                            tCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACAAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTg         <  2:858607/149‑1 (MQ=255)
                                                                                                                                               aTGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTGATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATg      <  2:500008/149‑1 (MQ=255)
                                                                                                                                                   ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta  <  1:660345/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTA  >  NC_000913/1449651‑1449946

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