Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,196,194 G→C 10.5% A147P (GCC→CCC)  PP_5618 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,196,1940GC10.5% 170.1 / 10.5 57A147P (GCC→CCC) PP_5618Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base C (2/4);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTG  >  NC_002947/4196050‑4196339
                                                                                                                                                |                                                                                                                                                 
cGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTgg                                                                                                                                              <  1:1070272/150‑1 (MQ=255)
       ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAggg                                                                                                                                                   >  1:818618/1‑138 (MQ=255)
       ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAggg                                                                                                                                                   <  2:818618/138‑1 (MQ=255)
       ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGt                                                                                                                                       >  2:1053127/1‑150 (MQ=255)
       ttGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGt                                                                                                                                       <  2:608216/150‑1 (MQ=255)
             tGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCa                                                                                                                                 >  2:736674/1‑150 (MQ=255)
              gggAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCAc                                                                                                                                <  1:1102979/150‑1 (MQ=255)
              gggAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTagacttt                                                                                                                                <  2:1133373/150‑7 (MQ=255)
                    aCAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAGGCCACTGCTCg                                                                                                                          >  2:290190/1‑150 (MQ=255)
                      aGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCt                                                                                                                        >  1:419631/1‑150 (MQ=255)
                        cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                      <  2:712563/150‑1 (MQ=255)
                        cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                      <  2:286723/150‑1 (MQ=255)
                        cAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTc                                                                                                                      >  1:712563/1‑150 (MQ=255)
                                    gCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCCAACTGTTTCGTTGCCATTGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCga                                                                                                          <  1:1053127/150‑1 (MQ=255)
                                      tGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGt                                                                                                                                             >  2:1186207/1‑113 (MQ=255)
                                      tGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGt                                                                                                                                             <  1:1186207/113‑1 (MQ=255)
                                          aGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAggg                                                                                                    >  1:633464/1‑150 (MQ=255)
                                            ggaggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTc                                                                                                                           <  1:972582/125‑1 (MQ=255)
                                            ggaggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTc                                                                                                                           >  2:972582/1‑125 (MQ=255)
                                               ggaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATg                                                                                               <  2:663983/150‑1 (MQ=255)
                                                gaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTggcaacg                                                                                                                                         <  2:365436/107‑6 (MQ=255)
                                                gaTCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTggcaacg                                                                                                                                         >  1:365436/1‑102 (MQ=255)
                                                  tCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGcc                                                                                            <  1:546162/150‑1 (MQ=255)
                                                     aaaaGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGa                                                                                         >  1:443526/1‑150 (MQ=255)
                                                            gCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTc                                                                                  <  2:419631/150‑1 (MQ=255)
                                                               cttCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCa                                                                               >  2:15832/2‑150 (MQ=255)
                                                                ttCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCAt                                                                              <  1:86424/150‑1 (MQ=255)
                                                                   aaGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATgag                                                                           <  2:332694/150‑1 (MQ=255)
                                                                    aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAc                                                                                                              >  2:684611/1‑114 (MQ=255)
                                                                    aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAc                                                                                                              <  1:684611/114‑1 (MQ=255)
                                                                    aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTa                                                                                                                                      >  2:97650/1‑90 (MQ=38)
                                                                    aGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTCACCAGGCCCTGGCAACGGTa                                                                                                                                      <  1:97650/90‑1 (MQ=37)
                                                                       aTGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAcgc                                                                                                            >  1:365237/1‑113 (MQ=255)
                                                                       aTGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTAcgc                                                                                                            <  2:365237/113‑1 (MQ=255)
                                                                        ccgatctcTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTCATCGCCCGACAAGCTCGCATGAGATCGc                                                                      <  1:1068436/143‑1 (MQ=255)
                                                                         gCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGcc                                                                     >  2:428875/1‑150 (MQ=255)
                                                                                tAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTAGCTTCTTCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGccagatcgg                                                              >  2:1068436/1‑143 (MQ=255)
                                                                                      ccTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGc                                                        <  1:512375/150‑1 (MQ=255)
                                                                                          aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                                    <  1:736674/150‑1 (MQ=255)
                                                                                          aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                                    >  1:896221/1‑150 (MQ=255)
                                                                                          aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                                    >  2:905884/1‑150 (MQ=255)
                                                                                          aCCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAg                                                    <  1:905884/150‑1 (MQ=255)
                                                                                                  agagaCAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGt                                            >  1:672806/1‑150 (MQ=255)
                                                                                                      aCAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGccc                                        >  2:687323/1‑150 (MQ=255)
                                                                                                        aGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAgg                                                   >  2:293315/1‑137 (MQ=255)
                                                                                                        aGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAgg                                                   <  1:293315/137‑1 (MQ=255)
                                                                                                         gCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTAt                                     <  1:15832/150‑1 (MQ=255)
                                                                                                          ccGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTa                                                                                                               >  1:1170232/1‑75 (MQ=255)
                                                                                                          ccGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTa                                                                                                               <  2:1170232/75‑1 (MQ=255)
                                                                                                           cGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATga                                                                            >  1:842373/1‑109 (MQ=255)
                                                                                                           cGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATga                                                                            <  2:842373/109‑1 (MQ=255)
                                                                                                                       gcTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCg                       >  2:450426/1‑150 (MQ=255)
                                                                                                                        cTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGt                      <  2:494620/150‑1 (MQ=255)
                                                                                                                                 tCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGc             >  1:216841/1‑150 (MQ=255)
                                                                                                                                   tACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCg           >  1:668604/1‑150 (MQ=255)
                                                                                                                                       gTTGCCAGGGCCTGGTGACGGTAGGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGa       <  1:290190/150‑1 (MQ=255)
                                                                                                                                           ccAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGTGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGc                                              >  1:41048/1‑107 (MQ=255)
                                                                                                                                           ccAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGTGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGc                                              <  2:41048/107‑1 (MQ=255)
                                                                                                                                            cAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGgtg  >  1:772575/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGGCAGCTTGAACTGGGAGGACAGCAAGCCCGACACGCTGCCAGGGAGGATCGAAAAGCTGCCATTCAAGCATGCAAGGCTAGGGGCCTCACCCAGCGAGAGACAGCCGCTGAGCTAGCGCTGTCAAAGTCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTG  >  NC_002947/4196050‑4196339

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: ‑