Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,993,895 C→T 7.5% E173K (GAA→AAA)  PP_2617 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,993,8950CT 7.5% 147.9 / 2.0 53E173K (GAA→AAA) PP_2617hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CAATCCGCGTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCTGTACCTGAG  >  NC_002947/2993748‑2994035
                                                                                                                                                    |                                                                                                                                            
cAATCCGCGTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTACAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTcc                                                                                                                                             <  8:141844/149‑1 (MQ=255)
 aaTCCGCGTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCa                                                                                                                                            <  1:134947/149‑1 (MQ=255)
 aaTCCGCGTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCa                                                                                                                                            <  4:322519/149‑1 (MQ=255)
     cgcgTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGaa                                                                                                                                         <  3:281156/148‑1 (MQ=255)
         tGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACGCCTGCAGCCCTGGCACCGCCCGCACCGAAACGCTGGGGGGGTAGAACGCGTCGTCCAGTTGCAGGCTGCCGGAGGGGCGCCGGTCCAGAATACGt                                                                                                                                    >  1:274577/1‑149 (MQ=255)
         tGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGt                                                                                                                                    <  7:229491/149‑1 (MQ=255)
             aaGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTgcca                                                                                                                                <  5:345755/149‑1 (MQ=255)
             aaGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTgc                                                                                                                                  >  7:152617/1‑147 (MQ=255)
               ggCCGCGCGCCCGCTCGCGCAAGGCATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGTACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTgccagg                                                                                                                              <  1:215183/149‑1 (MQ=255)
                    cgcgcCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTGGGGGCGCCGTTCCAGAATACGTgccagggccag                                                                                                                         <  8:152617/149‑1 (MQ=255)
                     gcgcCCGCTCGCGCAAGGTATGGGTGATATCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCCCACCGAAACGCTGGTGGGCTAGAACGCTACCTCCAGAAGCAGCCTGCCGTCGGGGCGCCCTTCCAGAATACGTgacagggccagg                                                                                                                        >  1:347447/1‑149 (MQ=255)
                           gCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGgtg                                                                                                                  >  7:132921/1‑149 (MQ=255)
                            cTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGgtgt                                                                                                                 <  3:108168/149‑1 (MQ=255)
                               gcgcAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAgg                                                                                                              >  3:30963/1‑149 (MQ=255)
                                 gcAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATa                                                                                                                                       <  6:181294/122‑1 (MQ=255)
                                    aGGGATGGGTGATGTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCgctg                                                                                                         <  3:261743/149‑1 (MQ=255)
                                     ggTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTTCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCgctgc                                                                                                        >  2:207840/1‑149 (MQ=255)
                                     ggTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTTCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCgctgc                                                                                                        >  2:213212/1‑149 (MQ=255)
                                      gTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCgctgct                                                                                                       <  3:161315/149‑1 (MQ=255)
                                                     tctAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCg                                                                                        >  2:174893/1‑149 (MQ=255)
                                                       tAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGgtg                                                                                                                  <  3:318708/121‑1 (MQ=255)
                                                       tAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGgtg                                                                                                                  >  4:318708/1‑121 (MQ=255)
                                                         gAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGttt                                                                                     <  8:132921/148‑1 (MQ=255)
                                                           aaCCGCTGCAGCGCTGGCACCGCCCTCACCGTCACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGGTGTCTTCGGTCTTGCGTTTGa                                                                                   <  2:347447/148‑1 (MQ=255)
                                                            aCCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGc                                                                                 <  2:274577/149‑1 (MQ=255)
                                                                cTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTTCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCa                                                                             <  1:213212/149‑1 (MQ=255)
                                                                cTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTTCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCa                                                                             <  1:207840/149‑1 (MQ=255)
                                                                 tGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAg                                                                            >  2:336868/1‑149 (MQ=255)
                                                                      cgcTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCCGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTgg                                                                       >  4:245690/1‑149 (MQ=255)
                                                                                 cccgcACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAgg                                                            <  6:34767/149‑1 (MQ=255)
                                                                                      aCCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTAc                                                       >  5:156981/1‑149 (MQ=255)
                                                                                          aaaCGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTg                                                    >  4:46855/1‑148 (MQ=255)
                                                                                             cGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCAcc                                                <  1:147066/149‑1 (MQ=255)
                                                                                              gCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGAGTAGGCGCTGGTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGTGCACCGCAAGGTCTAGCTGCAc                                                 >  8:148701/1‑147 (MQ=255)
                                                                                               cTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTggg                                                                      <  2:20532/125‑1 (MQ=255)
                                                                                               cTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTggg                                                                      >  1:20532/1‑125 (MQ=255)
                                                                                                tggtggGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGCTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAg                                              >  1:84610/1‑149 (MQ=255)
                                                                                                     gggTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCAcc                                        <  2:84610/149‑1 (MQ=255)
                                                                                                      ggTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGAAAGGTCTACCTGCACCAGGTCACCg                                       >  6:224585/1‑149 (MQ=255)
                                                                                                         aGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGt                                      >  4:40316/1‑147 (MQ=255)
                                                                                                          gAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAg                                   >  2:15703/1‑149 (MQ=255)
                                                                                                               cTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAg                              >  5:18058/1‑149 (MQ=255)
                                                                                                                         ttGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCt                    >  4:61033/1‑149 (MQ=255)
                                                                                                                           gCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCtt                   >  2:280188/1‑148 (MQ=255)
                                                                                                                            cAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTg                  <  3:245690/148‑1 (MQ=255)
                                                                                                                            cAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGa                 >  8:100439/1‑149 (MQ=255)
                                                                                                                               gCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAgg                                                            >  3:141156/1‑103 (MQ=255)
                                                                                                                               gCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAgg                                                            <  4:141156/103‑1 (MQ=255)
                                                                                                                               gCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTACAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGAttt              >  4:351159/1‑149 (MQ=255)
                                                                                                                                   ccGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCTGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCTg          >  4:70687/1‑149 (MQ=255)
                                                                                                                                   ccGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCt           <  6:156981/148‑1 (MQ=255)
                                                                                                                                     gTCGGGGCCCCGTTCCAGAATACGTGCCACGGCCAGGCGGGTGTATGCTCTGCTGTCTTCGGACTTGCCTTTGAGCTGCAGGTTGGGCACCGCAAGCTCTACCTGCACCAGGTCACCGACAGCGGAGAGGGTATCCTTGATATCCTGt         >  4:283286/1‑148 (MQ=255)
                                                                                                                                        ggggCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCTGTAcc      <  1:174893/148‑1 (MQ=255)
                                                                                                                                           gcgcCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCTGTACCTGAg  <  4:30963/149‑1 (MQ=255)
                                                                                                                                                    ccAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTgg                           <  6:127022/116‑1 (MQ=255)
                                                                                                                                                    ccAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTgg                           >  5:127022/1‑116 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CAATCCGCGTGGCAAGGCCGCGCGCCCGCTCGCGCAAGGTATGGGTGATTTCCTCTAGAAACCGCTGCAGCGCTGGCACCGCCCGCACCGAAACGCTGGTGGGGTAGAACGCTTCGTCCAGTTGCAGGCTGCCGTCGGGGCGCCGTTCCAGAATACGTGCCAGGGCCAGGCGGGTGTAGGCGCTGCTGTCTTCGGTCTTGCGTTTGAGCTGCAGGTTGGGCACCGCAAGGTCTACCTGCACCAGGTCACCGTCAGCGGAGTGGGTATCCTTGATTTCCTGTACCTGAG  >  NC_002947/2993748‑2994035

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: