Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,404,2150TG31.4% 73.4 / 23.9 51G244G (GGT→GGGyhdHacrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base G (16/0);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 6.75e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404082‑3404352
                                                                                                                                     |                                                                                                                                         
aGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGtttt                                                                                                                                     <  1:72941/140‑1 (MQ=255)
          aTCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGccca                                                                                                                           >  2:800507/1‑138 (MQ=255)
           tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAAc                                                                                                                          <  1:1136315/140‑1 (MQ=255)
                  ctctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                   <  1:692003/140‑1 (MQ=255)
                  ctctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                   <  1:614402/140‑1 (MQ=255)
                    ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGTGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  2:746516/140‑1 (MQ=255)
                    ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  2:1027733/140‑1 (MQ=255)
                    ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                 <  2:1169752/140‑1 (MQ=255)
                     tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCCGTGGTCTGGCGGGGGGGTTTTCTCCTCCAACCACGGGCATg                                                                                                                >  2:1136315/1‑140 (MQ=255)
                     tGGAAAAACAAGTCTGGGCTGGGACAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTg                                                                                                                               <  2:336156/125‑1 (MQ=255)
                     tGGAAAAACAAGTCTGGGCTGGGACAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTg                                                                                                                               >  1:336156/1‑125 (MQ=255)
                         aaaaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCAt                                                                                                            <  2:194453/140‑1 (MQ=255)
                               aGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                        <  2:649219/138‑1 (MQ=255)
                               aGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAt                                                                                                        >  1:649219/1‑138 (MQ=255)
                                 tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGt                                                                                                                    <  2:1212816/124‑1 (MQ=255)
                                 tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGTCATGCCCTTTAttttt                                                                                                    >  2:362613/1‑137 (MQ=255)
                                 tCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGt                                                                                                                    >  1:1212816/1‑124 (MQ=255)
                                        tgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATg                                                                                             <  1:87452/140‑1 (MQ=255)
                                          gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACCGGCATGCCATTTTTTTTGCGGAATGTc                                                                                           >  2:87452/1‑140 (MQ=255)
                                               aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTTTCTGGCAAACACCGTCATTGCCTTTTTTCTGCGGAATATCCGttt                                                                                      >  2:1210717/1‑140 (MQ=255)
                                               aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCCACCACGGGCATGCCCTTTTTTTTGCGTTATGTCCGgtt                                                                                      >  2:659135/1‑140 (MQ=255)
                                                   acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                  <  1:1210717/140‑1 (MQ=255)
                                                      ccGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCCTTTTTTTTGCGGAATGTCCGGTTGGAAggg                                                                               >  2:16040/1‑140 (MQ=255)
                                                        gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                             <  2:646460/140‑1 (MQ=255)
                                                        gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCAACCACGGGCATGCCCTTTTTTTTGCGGAAAGTCCGTTTGCAAagggg                                                                             >  2:861335/1‑139 (MQ=255)
                                                           ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:668906/140‑1 (MQ=255)
                                                           ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGa                                                                          <  1:1050467/140‑1 (MQ=255)
                                                           ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTTTTCTTCGTAATGTCCGTTTGCAAGggggggg                                                                          >  2:22868/1‑139 (MQ=255)
                                                           ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCAAACACGGGCATGCCATTTTTTTTGCGGAATGTTCCTTTGCAAGggggggg                                                                          >  2:358621/1‑139 (MQ=255)
                                                             cGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCCCGGGCATGCCCTTTTTTTTTCGTAATGTTCGTTTGCAAGGGGGGGGtt                                                                        >  2:668906/1‑140 (MQ=255)
                                                               aCAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGTGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCa                                                                      <  2:1010075/140‑1 (MQ=255)
                                                                   aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTTCCAACCACGGTCCTGCCATTTATTTTGCGTAATTTTCGGTTGCAAGGGGGGGGTTCAGTaa                                                                  >  2:14930/1‑140 (MQ=255)
                                                                    gTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAAt                                                                 <  1:358621/140‑1 (MQ=255)
                                                                       cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGAc                                                              <  1:362613/140‑1 (MQ=255)
                                                                       cTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGAc                                                              <  2:965818/140‑1 (MQ=255)
                                                                         ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGc                                                            <  2:1006044/140‑1 (MQ=255)
                                                                         ggCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCAAACACGGGCATGCCCTTTTTTTTGCGTAATGTTCGTTTGGAAGGGGGGGGTTTAGTAATGgcgc                                                            >  2:1039027/1‑140 (MQ=255)
                                                                               aGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTTCCAACCACCGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGccaccaa                                                      >  2:19437/1‑139 (MQ=255)
                                                                                    tGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcgc                                                 >  1:1023053/1‑140 (MQ=255)
                                                                                     ggCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcgcc                                                <  1:1039027/140‑1 (MQ=255)
                                                                                         cAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCgcgc                                            <  2:1104302/140‑1 (MQ=255)
                                                                                         cAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTACTCTTCCAACCCCGGTCATGCCATTTATTCTGCGGAATGTCCGGTTTCAAGGGGTGGGTTCAGGAATGACGCCCCCAGAAAGTCgcgc                                            >  2:659801/1‑140 (MQ=255)
                                                                                          aaaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGTCGCGca                                           <  1:659801/140‑1 (MQ=255)
                                                                                               aaTTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGc                                      <  2:69350/140‑1 (MQ=255)
                                                                                                      gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTTCTCTGCCAACCCCGGGCAAGCCCTTTATTTTTCGGAATGTCCGGTTGCAAGGGGGGGGTTCAGTAAAGACGCCACCAGAACGCCGCGCACAAACCCGGCAg                               >  2:200215/1‑140 (MQ=255)
                                                                                                       cggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGc                              <  1:800507/140‑1 (MQ=255)
                                                                                                         gcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCCTTTTTTTTGCGGAATGTCCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAACCTGGCAGCgg                            >  2:458111/1‑139 (MQ=255)
                                                                                                         gcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTTCAAGGGGGGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGa                            >  2:11527/1‑140 (MQ=255)
                                                                                                               tGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTc                      >  1:1265195/1‑140 (MQ=255)
                                                                                                               tGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTc                      >  1:628074/1‑140 (MQ=255)
                                                                                                                ggcggcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCg                     <  2:1101683/140‑1 (MQ=255)
                                                                                                                   ggcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCg                  <  1:339063/140‑1 (MQ=255)
                                                                                                                              tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:738984/140‑1 (MQ=255)
                                                                                                                              tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa       <  1:16040/140‑1 (MQ=255)
                                                                                                                               ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  1:863486/140‑1 (MQ=255)
                                                                                                                               ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt      <  2:1209649/140‑1 (MQ=255)
                                                                                                                                gCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc     >  1:830128/1‑140 (MQ=255)
                                                                                                                                   ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc            <  1:159995/130‑1 (MQ=255)
                                                                                                                                   ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc            >  2:159995/1‑130 (MQ=255)
                                                                                                                                   ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt  >  2:1024090/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                         
AGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATT  >  NC_000913/3404082‑3404352

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