Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,562,902 A→T 7.4% S308T (TCT→ACT)  PP_3146 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,562,9020AT 7.4% 151.2 / 3.4 54S308T (TCT→ACT) PP_3146oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (20/30);  new base T (2/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GGCGCCCGACCTGAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAGGGTTCGGTGGCCACCACATAG  >  NC_002947/3562758‑3563050
                                                                                                                                                |                                                                                                                                                     
ggCGCCCGACCTGAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAg                                                                                                                                                   >  1:282176/1‑149 (MQ=255)
   gcCCGACCTGAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCCCGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTgg                                                                                                                                                >  4:289114/1‑149 (MQ=255)
         ccTGAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGgcg                                                                                                                                          >  3:280907/1‑149 (MQ=255)
            gAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGTATAGTGGCAGGcgccg                                                                                                                                       >  3:167104/1‑149 (MQ=255)
            gAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGTATAGTGGCAGGcgccg                                                                                                                                       >  3:162830/1‑149 (MQ=255)
             aGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGcgccgc                                                                                                                                      >  7:238767/1‑149 (MQ=255)
               ggAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCgccgcc                                                                                                                                     >  2:90667/1‑148 (MQ=255)
               ggAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGAACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCgccgcc                                                                                                                                     >  6:255039/1‑148 (MQ=255)
                 aaaCGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGaa                                                                                                                                 >  5:49579/1‑149 (MQ=255)
                 aaaCGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGaa                                                                                                                                 >  6:16990/1‑149 (MQ=255)
                                 ggCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTGGGCCTCTTGCGGGATTACTTTCAGCACCTAGGGCCGCATATAGGCAGCGATGGCCTTGGGGTGGCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGcggcggcgc                                                                                                                 <  1:233967/148‑4 (MQ=18)
                                   ccGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGc                                                                                                               >  2:290540/1‑149 (MQ=255)
                                    cGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCg                                                                                                              <  6:49579/149‑1 (MQ=255)
                                    cGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCg                                                                                                              <  7:9079/149‑1 (MQ=255)
                                     gATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCgg                                                                                                             <  5:490/149‑1 (MQ=255)
                                      aTCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTTCGGt                                                                                                            <  1:18792/149‑1 (MQ=255)
                                      aTCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGt                                                                                                            <  4:280907/149‑1 (MQ=255)
                                         aTGCCACCCCACTGGTAGTCGATGCTCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGa                                                                                                          <  5:255039/148‑1 (MQ=255)
                                          tGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCTCACCGCAGACCGCCGCCTGCTGttcgg                                                                                                        <  3:97887/149‑5 (MQ=21)
                                           gCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGa                                                                                                                                                       <  7:122021/102‑1 (MQ=255)
                                           gCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGa                                                                                                                                                       >  8:122021/1‑102 (MQ=255)
                                            ccACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCg                                                                                                      <  5:333139/149‑1 (MQ=255)
                                             cACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCgg                                                                                                     >  5:142697/1‑149 (MQ=255)
                                             cACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCgg                                                                                                     >  3:190559/1‑149 (MQ=255)
                                                cccACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGgtag                                                                                                  <  8:238767/149‑1 (MQ=255)
                                                  cACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGTATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGgtagta                                                                                                <  4:167104/149‑1 (MQ=255)
                                                  cACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGTATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGgtagta                                                                                                <  4:162830/149‑1 (MQ=255)
                                                    cTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGc                                                                                                       <  5:102378/140‑1 (MQ=255)
                                                    cTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGc                                                                                                       >  6:102378/1‑140 (MQ=255)
                                                       gTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTcc                                                                                            <  7:61001/148‑1 (MQ=255)
                                                         aGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGt                                                                                                            >  7:238686/1‑130 (MQ=255)
                                                         aGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGt                                                                                                            <  8:238686/130‑1 (MQ=255)
                                                          gTCTATGCGCTCGTTGGCAAGTTGAGGGTATACTTTCAGCTCCTTGGGCAGCATTTTGGCTGCGTTGTCGTTGGGGTGTCGTCCTGATTAGTGGCAGGGGCGGCCGAACAGGATGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAgg                                                                                         <  2:282176/148‑1 (MQ=255)
                                                          gTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGc                                                                                        >  8:274720/1‑149 (MQ=255)
                                                           tCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGACACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGc                                                                                        <  1:62133/149‑1 (MQ=255)
                                                              aTGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGAc                                                                                    >  1:5981/1‑149 (MQ=255)
                                                                  gttgggtGGCCTGGTGCGGGATTACTTTCAGGACAGGGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTACAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTc                                                                                <  8:285197/143‑1 (MQ=255)
                                                                  ggaCGTGGGCCAGTTGCGGGAATACTTTCAGCACCTTTGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTc                                                                                <  4:190559/147‑1 (MQ=255)
                                                                    aCGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGTGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTg                                                                              <  8:109720/149‑1 (MQ=255)
                                                                     cGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTgg                                                                             <  7:42687/149‑1 (MQ=255)
                                                                       ttGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTc                                                                           <  1:290540/149‑1 (MQ=255)
                                                                          gCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGc                                                                         <  1:90667/148‑1 (MQ=255)
                                                                                tGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCg                                                                          <  2:69430/141‑1 (MQ=255)
                                                                                tGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCg                                                                          >  1:69430/1‑141 (MQ=255)
                                                                                       aTACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTc                                                           >  5:171319/1‑149 (MQ=255)
                                                                                           tttCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTc                                                                                             >  6:320132/1‑111 (MQ=255)
                                                                                           tttCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTc                                                                                             <  5:320132/111‑1 (MQ=255)
                                                                                                             cATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCa                                     <  6:171319/149‑1 (MQ=255)
                                                                                                               tataGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAgc                                   <  3:289114/149‑1 (MQ=255)
                                                                                                                    gCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTc                              <  7:274720/149‑1 (MQ=255)
                                                                                                                       gCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCggg                           <  2:5981/149‑1 (MQ=255)
                                                                                                                             tCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAggg                     >  7:104516/1‑149 (MQ=255)
                                                                                                                               cTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGctgctg                                         >  4:207707/1‑127 (MQ=255)
                                                                                                                               cTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGctgctg                                         <  3:207707/127‑1 (MQ=255)
                                                                                                                                    ggTGTCGTCCAGAATAGTGGCAGGCGCCGCCGATCAGCAGGCGGCGGTCTTCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAGGGTTCGGTg              <  4:227356/149‑1 (MQ=255)
                                                                                                                                    ggTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAGGGTTCGGTg              <  8:104516/149‑1 (MQ=255)
                                                                                                                                                aaTAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAGGGTTCGGTGGCCACCACATAg  <  5:16990/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGCGCCCGACCTGAGGGAAACGGTTGGCGGTGATGCCGATCATGCCACCCCACTGGTAGTCGATGCGCACGTTGGCCAGTTGCGGGAATACTTTCAGCACCTTGGGCCGCATATAGGCAGCGATGTCCTTGGGGTGTCGTCCAGAATAGTGGCAGGCGCCGCCGAACAGCAGGCGGCGGTCTGCGGTGAGGCGGTAGTAGTCCAGGCCGACCTTCTGGTCGCACAGCGCCATGTTCTGCGGGATCAGGCTGCTGGCCAGCGCTTCGGGCAAGGGTTCGGTGGCCACCACATAG  >  NC_002947/3562758‑3563050

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: