Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG56.2% 21.0 / 109.1 89V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (10/29);  new base G (50/0);  total (60/29)
Fisher's exact test for biased strand distribution p-value = 2.83e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-05
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCA  >  REL606/358‑508
                                                                            |                                                                              
 tACATCAGTGGCAAATGCGGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                               <  1:4257228/76‑1 (MQ=255)
  aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                              >  1:5440436/1‑76 (MQ=255)
   cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGggggg                                                                             >  2:1905628/1‑75 (MQ=255)
    aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCCGGGGCAGgggggc                                                                            >  1:2067983/1‑74 (MQ=255)
      cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAc                                                                          >  2:1138788/1‑76 (MQ=255)
       aGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAcc                                                                         <  2:548074/76‑1 (MQ=255)
       aGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAcc                                                                         >  2:6749514/1‑76 (MQ=255)
        gTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCGACCg                                                                       >  2:4676644/1‑76 (MQ=255)
          ggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCTAGGCAGGGGCAGGTGGCCACCGTc                                                                     <  1:1117715/76‑1 (MQ=255)
            cAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCct                                                                   <  1:6619510/76‑1 (MQ=255)
            cAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCct                                                                   <  1:3969808/76‑1 (MQ=255)
            aaaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCct                                                                   <  2:3404265/75‑1 (MQ=255)
             aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                                  <  2:6123433/76‑1 (MQ=255)
              aaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCctcc                                                                 >  2:2024973/1‑75 (MQ=255)
               aTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCACGGCctctc                                                               >  2:4753651/1‑76 (MQ=255)
                tGCAGACCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                              <  1:6953231/76‑1 (MQ=255)
                tGCAGACCGTTTGCTGCGGGTGGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGAAGGTGGCCACCGTCctctct                                                              <  2:1422430/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                              <  2:2243553/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                              <  1:4527794/76‑1 (MQ=255)
                 gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                             <  2:4480545/76‑1 (MQ=255)
                  cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGGCccctcgtc                                                            >  2:7178973/1‑73 (MQ=255)
                  cAGAACATTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCctccctcc                                                            >  2:5054981/1‑74 (MQ=255)
                     aaCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCCCCCCTCTTcccc                                                         >  2:4196578/1‑76 (MQ=255)
                     aaCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCATCCTCTCCTcccc                                                         >  2:2615505/1‑76 (MQ=255)
                      aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                        <  1:4869509/76‑1 (MQ=255)
                      aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCTCCCCTCCTccccc                                                        >  2:498835/1‑76 (MQ=255)
                       cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCg                                                       <  2:2161301/76‑1 (MQ=255)
                       cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCCTCCTCTCTGtccccc                                                       >  1:3422962/1‑75 (MQ=255)
                        gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACccgcccccccgccccccc                                                      >  2:3774875/1‑74 (MQ=255)
                        gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCCCTCTGccccccc                                                      >  1:2263309/1‑74 (MQ=255)
                        gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTTCCTTccccccc                                                      >  2:4335437/1‑74 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCa                                                   <  1:367474/76‑1 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCGCCTCTCCTccccccccc                                                    >  1:2195894/1‑75 (MQ=255)
                           ttCTGGGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCaa                                                  <  2:4903561/76‑1 (MQ=255)
                           ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCaa                                                  <  1:1132394/76‑1 (MQ=255)
                             cTGCGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCaaaa                                                <  1:5326562/76‑1 (MQ=255)
                             cTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAACGGCCCCTCCTCCCCCCGCCaaa                                                 >  1:4371873/1‑76 (MQ=255)
                              tGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCCCTCTGCCCCCCCCcaaaa                                               >  2:1996432/1‑75 (MQ=255)
                               gCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCGCCCCTCTTCCCCCCCCcaaaac                                              >  1:5450783/1‑74 (MQ=255)
                                 gggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGGCCTCTCTCCCCCCGCCAAAAacac                                            >  2:1308585/1‑76 (MQ=255)
                                 gggTTGCCGATATTCTGAAAAGCAAGGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCac                                            <  2:3109797/76‑1 (MQ=255)
                                  ggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTGCCCCCGCCCAAATCacc                                           >  1:5201369/1‑76 (MQ=255)
                                     tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCCTTTCCCCCCGCCAAAACAaccac                                         >  2:2205203/1‑75 (MQ=255)
                                       ccGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaaccc                                      >  1:3517191/1‑75 (MQ=255)
                                       ccGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCCCTCTGCCCCCCCCCAAAACaaccacca                                      >  2:6850544/1‑76 (MQ=255)
                                        cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCCCTTCCCCCCCCCAAAACaaccacccc                                     >  1:2099292/1‑74 (MQ=255)
                                        cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAACCTCCTCTCTGCCCCCCCCAAAAACAACAaccac                                     >  2:826558/1‑76 (MQ=255)
                                         gATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCCCCCAaaacaccaaacacc                                    >  2:4600792/1‑76 (MQ=255)
                                         gATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGCCCTCTCTTCCCCCCCCCAAATAACCAaccccc                                    >  1:1605011/1‑76 (MQ=255)
                                          atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCCGCCCCTCTTCCCCCCCCCAAAACACCCaacaccc                                   >  1:1945492/1‑75 (MQ=255)
                                            tttCTGAAAAGCAATCCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtgg                                 <  2:3580705/75‑1 (MQ=255)
                                            atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCGCCAAAAACAACAaacaccccg                                 >  2:2729575/1‑73 (MQ=255)
                                            atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGTTCTCTCCCCCCCCCGCCAAAAACACAAcccaccgg                                  >  1:23032/1‑74 (MQ=255)
                                             ttCTGGAAAGCAATGCCGGGCGGGGGCGGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggt                                <  1:5712482/76‑1 (MQ=255)
                                              tCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCAAAATCACCAACCACCcggtg                               >  2:3264391/1‑76 (MQ=255)
                                               cTGGAAAGCAATGCCGGGCAGGGGCGGGTGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCtggtgg                              <  1:5554924/76‑1 (MQ=255)
                                               cTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCTCCTCTCTGCCCCCGCCAAAAACACAAACCAACtggagg                              >  2:5186152/1‑76 (MQ=255)
                                                tGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCTCCCTCCTCTCTGCCCCCACCAAAATCTCCAACCACCtggaggg                             >  1:6682633/1‑75 (MQ=255)
                                                tGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCAAAAACCACCACCACCttgtggg                             >  2:1585282/1‑75 (MQ=255)
                                                 ggAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                            <  2:3171775/76‑1 (MQ=255)
                                                 ggAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTTCCCCCGCCAAAATCACCAACCACCTGGTGGCg                            >  2:7027724/1‑76 (MQ=255)
                                                 ggAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGCGATCTCTCCCCCCGCCAAAAATTACCACCCCCTGgcggcg                            >  2:3347311/1‑76 (MQ=255)
                                                 gaaaaGCAATGCCAGGCAGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                            <  1:7109602/74‑1 (MQ=255)
                                                 gaaaaGCAAGGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                            <  2:2514475/74‑1 (MQ=255)
                                                  gAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTGCCCCCGCCCAAAACACCAACCACCCGGTGGCgg                           >  2:4245519/1‑75 (MQ=255)
                                                  gAAAGCAATGCCAGGAAGGGGCAGGTGGCCTCCGTCCTCTCTGCCTCCGCCAAAATCACCAACCACCTGGTGGCga                           <  2:2993162/76‑1 (MQ=255)
                                                   aaaGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTGCCCCCCCCCAAATAACCCACCACCTCGGGGCgat                          >  2:2528315/1‑76 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCGGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                       <  2:5918807/76‑1 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTCCCCCCGCCCAAATCACCAAACACCTGGTGGCgatgct                       >  2:769225/1‑74 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCGCCAAAAACACCAACCACCCGGTGGCggtgat                       >  2:145578/1‑76 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                       >  2:7091840/1‑76 (MQ=255)
                                                      ccAATCCCAGGCAGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                       <  1:6852792/75‑1 (MQ=255)
                                                       cAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAATCACCCACCACCTTGTGGCGGAGAtt                      >  2:4655970/1‑76 (MQ=255)
                                                        aaggCCAGGAGGGGGCGGGTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                     <  2:2040198/73‑1 (MQ=255)
                                                        aaTGCCAGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                     <  1:7009193/76‑1 (MQ=255)
                                                         aTGCCGGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                    <  1:3394919/76‑1 (MQ=255)
                                                         aTGCCAGGCGGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                    <  1:6510585/76‑1 (MQ=255)
                                                         aTGCCAGGCAGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                    <  1:4097914/76‑1 (MQ=255)
                                                         aTGCCAGGCAGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                    <  1:5307149/76‑1 (MQ=255)
                                                         aTGCCAGGCAGGGGCAGGGGGCCCCCGCCCTCCCTGACCCCGCCAAAAACACCAACCCCCCGCGGGGGttgattga                    >  2:7121882/1‑76 (MQ=255)
                                                            ccAGGCAGGGGCAGGGGGGCCCCCGCCCCTTTTACCCCCCCCAAAATACCCACCACCCGGTGGCGAAGAGTGgaaa                 >  2:5313970/1‑76 (MQ=255)
                                                            ccAGGCAGGGGCAGGGGGGCACCCTCCCCTCTTCCCCCCCCCAAAAAACCAACCCCCTGGTGGGGGTGATTTaaaa                 >  2:381343/1‑76 (MQ=255)
                                                             caggcaggGGCAGGGGGGCCCCCGCCTCTCTGCCCCCCCCCAAATCACCAAACACCTGGTGGCGCTGGTTGaaaaa                >  2:4738738/1‑76 (MQ=255)
                                                             caggcaggGGCAGGGGGGCACCCTCCTCTCTGGCCCCCCCCAAAACACCAACCACCTGGTGGCGATGGTTGaaaaa                >  2:1141576/1‑76 (MQ=255)
                                                              aggcaggGGCAGGGGGGCACCGGCCTCTCTTCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaaa               >  2:6848368/1‑76 (MQ=255)
                                                               ggcgggggcaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAc              <  1:144930/76‑1 (MQ=255)
                                                               ggcaggggcaggTGCCCCCCGCCCTTTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAc              <  1:3055971/76‑1 (MQ=255)
                                                               ggcaggggcaggGGGGCCCCGTCCTCTCTTCCCCCCCCAAAATCACCAAACCCCTGGTGGCGATGAATGAAAAAAc              >  2:1777314/1‑76 (MQ=255)
                                                               ggcaggggcaggGGGGCACCCCCCTCTCTTCCCCCCCCCAAAACAACAAACACCCGGTGGCGGTGATTTaaaaaaa              >  2:2706786/1‑75 (MQ=255)
                                                                gcgggggcgggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAcc             <  1:6471215/73‑1 (MQ=255)
                                                                 caggggcaggggGGCACCGTCCTCTCTGCCCCCCCCCAAAACACCAACCCCCTGGTGGCGATGGTTTAAAAAAACa            >  2:6562143/1‑76 (MQ=255)
                                                                 caggggcaggggGGCACCGGCCTCTCTGCCCCCCCCAAAAACACCCACCACCTCGGGGCGATGATTGAAAAAAACa            >  2:164334/1‑76 (MQ=255)
                                                                 caggggcaggggGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCa            >  2:5212637/1‑76 (MQ=255)
                                                                   ggggcggtTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:5983440/76‑1 (MQ=255)
                                                                   ggggcgggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:2126710/76‑1 (MQ=255)
                                                                   ggggcaggGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAACAtt          >  2:6442328/1‑76 (MQ=255)
                                                                   ggggcaggGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:3389804/76‑1 (MQ=255)
                                                                       caggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg      >  1:2842696/1‑76 (MQ=255)
                                                                        aggTGCCCCCCGCCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg     <  1:959264/76‑1 (MQ=255)
                                                                           gTGCCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGcc   <  1:558594/76‑1 (MQ=255)
                                                                            |                                                                              
GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCA  >  REL606/358‑508

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: