Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 770,931 G→A 9.1% L204L (CTG→CTA)  PP_0662 → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947770,9310GA 9.1% 144.0 / 3.2 52L204L (CTG→CTA) PP_0662threonine synthase
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base A (3/2);  total (27/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 = 1.00e+00

CCCTGCAGCAGGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTGGGTTCAGCATCCCCACGG  >  NC_002947/770784‑771079
                                                                                                                                                   |                                                                                                                                                    
cccTGCAGCAGGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCgctgct                                                                                                                                                    <  2:743585/150‑1 (MQ=255)
 ccTGCAGCAGGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGctgctg                                                                                                                                                   >  1:922753/1‑150 (MQ=255)
    gcagcaGGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGt                                                                                                                                                <  1:452267/150‑1 (MQ=255)
     cagcagGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGAGGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTg                                                                                                                                               >  1:772735/1‑150 (MQ=255)
      agcagGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGa                                                                                                                                              <  1:764963/150‑1 (MQ=255)
        cagGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACt                                                                                                                                            >  2:737857/1‑150 (MQ=255)
          ggCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCATGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTgg                                                                                                                                          >  1:227825/1‑150 (MQ=255)
                 tGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCAGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGcc                                                                                                                                   >  2:69987/1‑150 (MQ=255)
                     gcTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCgg                                                                                                                               <  1:737857/150‑1 (MQ=255)
                     gcTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCgg                                                                                                                               <  2:532513/150‑1 (MQ=255)
                          aCCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGa                                                                                                                          >  2:301/1‑150 (MQ=255)
                            cccAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGAtt                                                                                                                        >  1:969575/1‑150 (MQ=255)
                             ccAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCggg                                                                                                                              >  2:877864/1‑143 (MQ=255)
                             ccAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCggg                                                                                                                              <  1:877864/143‑1 (MQ=255)
                               aGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCg                                                                                                                     <  1:700879/150‑1 (MQ=255)
                                 cgcAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCgga                                                                                                                   <  1:994510/150‑2 (MQ=255)
                                        aCGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTACTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTc                                                                                                            >  1:582348/1‑150 (MQ=255)
                                                cgACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTcca                                                                                                     >  1:821180/1‑149 (MQ=255)
                                                cgACCGCGCGGCAGTGAGGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTcca                                                                                                     <  2:821180/149‑1 (MQ=255)
                                                      cgcgGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTACTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTg                                                                                              >  2:432707/1‑150 (MQ=255)
                                                       gcgGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCATGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTgg                                                                                             <  2:227825/150‑1 (MQ=255)
                                                            aGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCgg                                                                                        >  1:363370/1‑150 (MQ=255)
                                                             gTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGgt                                                                                       <  2:922753/150‑1 (MQ=255)
                                                                  gcgcAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCt                                                                                  >  1:641500/1‑150 (MQ=255)
                                                                        tGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCa                                                                            <  2:574731/150‑1 (MQ=255)
                                                                               tCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATtgtg                                                                     <  2:363370/150‑1 (MQ=255)
                                                                               tCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCAACCAACTGGGTCGGTGTGCTGGCGCAGATtgtg                                                                     <  1:360425/150‑1 (MQ=255)
                                                                                  gACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTc                                                                  >  1:882805/1‑150 (MQ=255)
                                                                                   aCAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCt                                                                 >  1:139364/1‑150 (MQ=255)
                                                                                      gCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTAtt                                                              >  2:941745/1‑150 (MQ=255)
                                                                                        gTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTc                                                            >  2:1176241/1‑150 (MQ=255)
                                                                                            gCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATg                                                        >  1:1129774/1‑150 (MQ=255)
                                                                                             cTGCGTGGCCGTGGATGCCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGCTTGTGTTCTATTTCCATgc                                                       >  1:101010/1‑150 (MQ=255)
                                                                                              tGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATgcg                                                      >  1:780069/1‑150 (MQ=255)
                                                                                                    gCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGgc                                                                               <  1:600688/119‑1 (MQ=255)
                                                                                                    gCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGgc                                                                               >  2:600688/1‑119 (MQ=255)
                                                                                                    gCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTg                                                >  2:462060/1‑150 (MQ=255)
                                                                                                    gCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTg                                                >  1:1094722/1‑150 (MQ=255)
                                                                                                      cGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGGTGCGTGACTGGCCATGCCCCGGGCTGATTTTGGTCAGCTTCAACTCGACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTgaa                                              >  1:486432/1‑148 (MQ=255)
                                                                                                      cGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGt                                                                                                                   >  2:398407/1‑81 (MQ=255)
                                                                                                      cGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGt                                                                                                                   <  1:398407/81‑1 (MQ=255)
                                                                                                              gCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTg                                      <  1:941745/150‑1 (MQ=255)
                                                                                                                    ttcactaTTGCCACTCGCTGTTTTCGGCGCTGCTGCGTGACTGCCCATGCCCCGGGCTTATTCCGGTCAGCTTCAACTCCACCAACCGGGTCGGTGTGCTCTCGCAGATTGTGTTCTATCTCCATCCGGCTCTGCACCTCGGTGGCGGCa                                <  2:486432/144‑1 (MQ=255)
                                                                                                                     tGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCAt                               <  2:1129774/150‑1 (MQ=255)
                                                                                                                     tGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCAt                               <  1:462060/150‑1 (MQ=255)
                                                                                                                      gACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATa                              >  2:870333/1‑150 (MQ=255)
                                                                                                                          aTTGCCAGTCGCTGGTGTCTGCTCTACGGCGTGCCTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGcc                          >  1:1147158/1‑150 (MQ=255)
                                                                                                                              ccAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGt                      >  2:24789/1‑150 (MQ=255)
                                                                                                                              ccAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGt                      <  1:69987/150‑1 (MQ=255)
                                                                                                                               cAGTCGCTGGTGTCGGCGCTACTGCGTGACTGGTCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGtt                     <  1:432707/150‑1 (MQ=255)
                                                                                                                                 gTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTgg                   <  2:780069/150‑1 (MQ=255)
                                                                                                                                  tCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTggg                  <  2:969575/150‑1 (MQ=255)
                                                                                                                                  tCGCTGGTGTCGGCGCTACTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTggg                  <  2:582348/150‑1 (MQ=255)
                                                                                                                                       ggTGTCGGCGCTGCTGCGTGACTGGCCCTGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTGGGTTCAg             <  2:570925/150‑1 (MQ=255)
                                                                                                                                                gctgctGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTGGGTTCAGCATCCCCAc    >  1:314602/1‑150 (MQ=255)
                                                                                                                                                  tgctgcGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTGGGTTCAGCATCCCCACgg  <  2:139364/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCTGCAGCAGGCTTGCTGGCGCTATACCCCAGCGCAGGCACGCCCGCCGACCGCGCGGCAGTGATGCGCAGTGCCGCGTCCGACAGCGTGAGCTGCGTGGCCGTGGATGGCAGCTTTGACGATTGCCAGTCGCTGGTGTCGGCGCTGCTGCGTGACTGGCCATGCCCCGGGCTGATTCCGGTCAGCTTCAACTCCACCAACTGGGTCGGTGTGCTGGCGCAGATTGTGTTCTATTTCCATGCGGCGCTGCAGCTCGGTGGCGGCATACGCCCGGTTGGGTTCAGCATCCCCACGG  >  NC_002947/770784‑771079

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