Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,315,681 A→G 7.0% F160L (TTC→CTC)  PP_2915 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,315,6810AG 7.0% 159.5 / 2.7 57F160L (TTC→CTC) PP_2915membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base G (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTGCGGCGAAAGC  >  NC_002947/3315537‑3315825
                                                                                                                                                |                                                                                                                                                
tGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCACCGCCCCAGGCCGCGGCCGAAAcgc                                                                                                                                              >  3:176785/1‑149 (MQ=255)
  ccGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGa                                                                                                                                            >  3:62361/1‑149 (MQ=255)
    ggCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGa                                                                                                                                            >  6:39331/1‑147 (MQ=255)
    ggCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACa                                                                                                                                          <  6:111847/149‑1 (MQ=255)
         gcgcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcg                                                                                                                                     <  2:77383/149‑1 (MQ=255)
         gcgcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcg                                                                                                                                     <  6:83178/149‑1 (MQ=255)
           gcCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGcc                                                                                                                                   <  5:39331/149‑1 (MQ=255)
             cTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGATAGGGCGCCCg                                                                                                                                 >  6:16917/1‑149 (MQ=255)
                cAGGCCCAGGTGTAGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGt                                                                                                                                <  8:276362/147‑1 (MQ=255)
                       aGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGt                                                                                                                       <  8:108417/149‑1 (MQ=255)
                                 gCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAATCAGGCCGATg                                                                                                             >  1:260745/1‑149 (MQ=255)
                                    gctgggCGCCGCGCTCGGCGCCTAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                                                                                                                         <  2:240092/130‑1 (MQ=255)
                                    gctgggCGCCGCGCTCGGCGCCTAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTaa                                                                                                                         >  1:240092/5‑134 (MQ=255)
                                       aGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAgg                                                                                                       <  7:29588/149‑1 (MQ=255)
                                           gcCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCAt                                                                                                   <  3:124833/149‑1 (MQ=255)
                                              gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGAt                                                                                                              >  5:141502/1‑135 (MQ=255)
                                              gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGAt                                                                                                              <  6:141502/135‑1 (MQ=255)
                                              gAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCa                                                                                                >  3:47206/1‑149 (MQ=255)
                                                gcgcGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCAcc                                                                                              <  3:66190/149‑1 (MQ=255)
                                                 cgcgGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCt                                                                                             <  2:260745/149‑1 (MQ=255)
                                                       gcCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                                                       >  6:300342/1‑149 (MQ=255)
                                                        cccAGCATCATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCAc                                                                                      <  5:300342/149‑1 (MQ=255)
                                                         ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                                                       <  5:231855/147‑1 (MQ=255)
                                                         ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                                                       <  5:235878/147‑1 (MQ=255)
                                                         ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                                                       >  6:231855/1‑147 (MQ=255)
                                                         ccAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTAGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAGACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                                                                       >  6:235878/1‑147 (MQ=255)
                                                           aGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGAc                                                                                   >  1:272773/1‑149 (MQ=255)
                                                                        cccAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                                                      >  1:186045/1‑149 (MQ=255)
                                                                          cAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGcc                                                                    >  2:139790/1‑149 (MQ=255)
                                                                              gCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATg                                                                >  5:100923/1‑149 (MQ=255)
                                                                                        cACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCAAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTgg                                                      >  2:258112/1‑149 (MQ=255)
                                                                                         aCCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgc                                                     >  4:299877/1‑149 (MQ=255)
                                                                                          cccGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgcc                                                    <  1:225667/149‑1 (MQ=255)
                                                                                           ccGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccg                                                   <  3:108139/149‑1 (MQ=255)
                                                                                             gCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccgcc                                                 <  4:62361/149‑1 (MQ=255)
                                                                                                         gccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                                                      <  3:35179/116‑1 (MQ=255)
                                                                                                         gccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATg                                                                      >  4:35179/1‑116 (MQ=255)
                                                                                                          ccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAggg                                     >  2:104674/1‑148 (MQ=255)
                                                                                                            gcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAg                                  >  4:19461/1‑149 (MQ=255)
                                                                                                            gcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAg                                  <  2:272773/149‑1 (MQ=255)
                                                                                                             cTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAgg                                 <  3:19461/149‑1 (MQ=255)
                                                                                                                gCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGt                                <  1:139790/147‑1 (MQ=255)
                                                                                                                gCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTgg                              >  3:294275/1‑149 (MQ=255)
                                                                                                                    agcagcCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCgg                          <  6:100923/149‑1 (MQ=255)
                                                                                                                     gcagcCAACGCCCCAGGCCGCGGCCAAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGc                         <  1:258112/149‑1 (MQ=255)
                                                                                                                           aaCGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGc                   >  4:284020/1‑149 (MQ=255)
                                                                                                                           aaCGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGc                   >  4:279503/1‑149 (MQ=255)
                                                                                                                           aaCGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATCCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGc                   >  7:237903/1‑149 (MQ=255)
                                                                                                                                cccAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAgg                                      <  5:86382/125‑1 (MQ=255)
                                                                                                                                cccAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAgg                                      >  6:86382/1‑125 (MQ=255)
                                                                                                                                cccAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCg                           >  4:101412/1‑136 (MQ=255)
                                                                                                                                cccAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCg                           <  3:101412/136‑1 (MQ=255)
                                                                                                                                 ccAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCt             >  4:289951/1‑149 (MQ=255)
                                                                                                                                    ggCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTgcg          <  1:104674/149‑1 (MQ=255)
                                                                                                                                     gCCGCGGCCGAACCGCGACAGGGCGCCGGTTACGTACGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCTTTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTgcgg         <  4:47206/149‑1 (MQ=255)
                                                                                                                                            ccGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTGCGGCGAAAGc  <  3:289951/149‑1 (MQ=255)
                                                                                                                                             cGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTGCGGCGAAGGc  >  2:228290/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAGCGCAGGGAAGGTGGCCGGCAACGGCCACGCTGCGGCGAAAGC  >  NC_002947/3315537‑3315825

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