Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG38.3% 147.1 / 36.8 95V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (17/41);  new base G (36/0);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 2.78e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 199 total aligned reads displayed.

TACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCC  >  REL606/359‑507
                                                                            |                                                                             
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg                                                                            <  2:4422939/76‑1 (MQ=255)
  cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGggggg                                                                           >  1:529848/1‑75 (MQ=255)
  cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                                                                           <  2:2537588/76‑1 (MQ=255)
  cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                                                                           <  1:3029573/76‑1 (MQ=255)
   aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                          <  1:3182721/76‑1 (MQ=255)
    tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                         <  2:1599604/76‑1 (MQ=255)
     cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAc                                                                        >  1:6226840/1‑76 (MQ=255)
        tGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGt                                                                     >  1:4045789/1‑76 (MQ=255)
         ggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTc                                                                    <  1:5816771/76‑1 (MQ=255)
          gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGccc                                                                   >  1:6548360/1‑76 (MQ=255)
          gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGGcc                                                                   >  1:593814/1‑76 (MQ=255)
           cAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCct                                                                  >  1:8073089/1‑76 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCACCGTcccc                                                                 >  2:1934349/1‑74 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                                 <  2:8020355/76‑1 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                                 <  2:6804659/76‑1 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTcccc                                                                 >  2:730541/1‑74 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGACctccctcccc                                                                 >  2:6714847/1‑74 (MQ=255)
             aaTGCAGAACGTTTTCTGCGGGTTGCCGATGTTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctct                                                                <  1:3809032/76‑1 (MQ=255)
             aaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctct                                                                <  1:7210178/76‑1 (MQ=255)
             aaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctct                                                                <  1:7934689/76‑1 (MQ=255)
               tGCAGACCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  2:7524847/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  1:7404513/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  1:1572128/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  2:2349060/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCGCCGGCctcttc                                                              >  1:5524545/1‑74 (MQ=255)
                gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                            <  2:112186/76‑1 (MQ=255)
                gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                            <  2:4607660/76‑1 (MQ=255)
                 cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCTCctctctcc                                                           >  2:2311599/1‑74 (MQ=255)
                 cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCccctcttc                                                           >  2:3609230/1‑74 (MQ=255)
                 cAGAACGCTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGc                                                           <  2:8006649/76‑1 (MQ=255)
                  aGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCCCCTCTTcc                                                          >  2:3437336/1‑76 (MQ=255)
                      cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCGCCGCCCTCTCTGcccccc                                                      >  1:1184397/1‑75 (MQ=255)
                       gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCCTccccccc                                                     >  2:5656055/1‑74 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGcc                                                    <  1:3138322/76‑1 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGCCCTCTCTTcccccccc                                                    >  1:5959921/1‑76 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCTCTCCCCCCGcc                                                    >  2:5519190/1‑76 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCCTCTCCGcccccccc                                                    >  2:2683266/1‑76 (MQ=255)
                         tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGGCCTCTCTTCCCCCCCCa                                                   >  1:1751639/1‑76 (MQ=255)
                         tttCTGCGGGTGGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCa                                                   <  2:3080092/76‑1 (MQ=255)
                            cTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCcaaa                                               >  2:6432171/1‑76 (MQ=255)
                            cTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCTCTTCCCCCCCCcaaa                                               >  1:6751768/1‑76 (MQ=255)
                              gggggTGGCCGATATTCTGAAAAGCAATGCAGGGCGGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATc                                             <  2:1442780/74‑1 (MQ=255)
                              gCGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATc                                             <  1:1485918/76‑1 (MQ=255)
                               ggtgTTGCAGATATTCGGGAAAGCAAGGCCGGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCa                                            <  2:658596/73‑1 (MQ=255)
                               cGGTTTGCCGATATCCTGAAAAGAAATGCCAGGCAGGGGCAGGTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCa                                            <  1:7660114/76‑1 (MQ=255)
                                  tttGCCAATATTCGGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCacca                                         <  2:2185152/75‑1 (MQ=255)
                                  gTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAATCacca                                         >  2:5139547/1‑76 (MQ=255)
                                  gTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGGCCTCTCTGCCCCCGCCAAAATTacca                                         >  1:2704965/1‑76 (MQ=255)
                                   ttcccGATATTCTGGAAAGCAATGCCAGCCAGGGGCGGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                        <  1:969141/73‑1 (MQ=255)
                                   ttGCCGATATTCTGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                        <  1:1305874/76‑1 (MQ=255)
                                   ttGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCCCCCCTCCCTCTCCCCCCCCCCAAAacaacac                                        >  2:4757971/1‑75 (MQ=255)
                                   ttGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTTTCTGGCCCCCCCAAAATAcccaa                                        >  2:1245418/1‑76 (MQ=255)
                                     gCCAATATTCTGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCACCGCCCTCTCTGCCCCCGCCAAAATCaccaacc                                      <  2:7804620/76‑1 (MQ=255)
                                      ccGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCGCCCCCCCGCCCCCCCCCAAAACaaccaacc                                      >  2:6486414/1‑76 (MQ=255)
                                       cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCACCGCCAAAACGaccaacccc                                    >  1:1969941/1‑74 (MQ=255)
                                       cGATATTCTGGAAAGCAATGCCAGGCAGGGCCAGGGGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAaccac                                    <  2:8129669/76‑1 (MQ=255)
                                       cGATATTCTGAAAAGCAATGCCGGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAaccac                                    <  1:6800783/76‑1 (MQ=255)
                                           tttCTGGAAAGCAACGCCAGGCAGGGGCAGGTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtgg                                <  2:7909786/75‑1 (MQ=255)
                                           atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCCCCAAAAAAACCAacccccccg                                >  2:7874502/1‑73 (MQ=255)
                                            ttCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggt                               <  1:4018795/76‑1 (MQ=255)
                                             tCTGAAAACAAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtg                              <  2:1751639/76‑1 (MQ=255)
                                              cTGGAAAGCAATGCCGGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtgg                             <  2:5002704/76‑1 (MQ=255)
                                               tGGAAAGCAATGCCAGGCAGGGGCAGGGGGGGCGCGGCCCCTCCGCCCCCCCCAAACAAGACCAAAACCTGGTGGc                            >  2:1421196/1‑76 (MQ=255)
                                               tGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCCTGCCCCCCCCAAAACACCAAACACCCCGTGGc                            >  1:5081575/1‑76 (MQ=255)
                                               gggAAAGCAATGCCAGGCAGGGCCGGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                            <  1:5028524/75‑1 (MQ=255)
                                                 gAAAGCAATGCCAGGCAGGGGCAGGGGGGGACCCGCCTCTCTTGCCCCACCAAAAACACCAACCCCCAGGAGGCgc                          >  2:6271378/1‑75 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCCCCCTGCCCCCGCCAAAATCACCCACCACCCGGTGGCggt                         >  2:2892867/1‑74 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGTGGCCACCCGCCCCTCCGCCCCCCCCCAAAAAACCAACCACCCGGTGGCgat                         >  2:3095333/1‑76 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAAACAACCACCACCTGGTGGGgat                         >  2:6049892/1‑76 (MQ=255)
                                                   aaGCAATGCCAGGCGGGGGAGGGGGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                        <  1:7792333/76‑1 (MQ=255)
                                                   aaGCAATGCCAGGCAGGGGCGGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                        <  1:3220700/76‑1 (MQ=255)
                                                    aGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatga                       >  2:3843620/1‑76 (MQ=255)
                                                     gaaaTCCCGGGCAGGGGCGGGTGGCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                      <  1:4995543/74‑1 (MQ=255)
                                                      cAATGCCAGGCAGGGGCAGGGGGCCACCCTCCCCCCCTCCCCCGCCCAAAACAACAACCCCCTGGTGGCGGAGAtt                     >  2:3668276/1‑76 (MQ=255)
                                                       aaTGCCAGGCAGGGGCAGGTGGCCACCCTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGAGATGATTg                    >  2:2023655/1‑76 (MQ=255)
                                                       aaTGCCAGGCAGGGGCAGGGGGGCACCCTTCCCCCTTCCCCCCCCCAAAACACCCACCACCCGGTGGGggtgaatt                    >  2:276539/1‑73 (MQ=255)
                                                       aaTGCCAGGCAGGGGCAGGGGGCCCCCGCCCTCCCTTGCCCCGCCCAAAATACCCACCACCATGTGGCGGTGGTTg                    >  2:1528516/1‑76 (MQ=255)
                                                        aTGCCGGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                   <  2:1272894/76‑1 (MQ=255)
                                                        aTGCCCGGGAGGGGCCAGTTGGCCCCCGCCCCCTCGGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                    <  1:6981053/76‑1 (MQ=255)
                                                        aTGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTTCGCCCCCGCCCAAATCACCCACCACCTGGTGGCGATGTTTgt                   >  2:1805763/1‑75 (MQ=255)
                                                         tGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTGCCCCCGCCACAAACACCAACCACCCGGTGGCGCTGCTTGaa                  >  2:6896806/1‑76 (MQ=255)
                                                          gcgagggagggGAAGGTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaa                 <  2:7409513/73‑1 (MQ=255)
                                                          gCCGGGCGGGGACAGTTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaa                 <  1:8151333/76‑1 (MQ=255)
                                                          gCCAGGCGGGGGCGGGGGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaa                 <  1:2991379/76‑1 (MQ=255)
                                                          gCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAAACACCAAACACCCGGTGGCGCTGATTGaaa                 >  2:4225803/1‑76 (MQ=255)
                                                           ccAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTTCCCCCGCCCAAATCACCCAACACCTGGTGGGGATGATTGaaaa                >  2:5860796/1‑76 (MQ=255)
                                                            caggcaggGGCAGGGGGCCACCGTCCTCTCTGCCCCCCCCCAAATCACCCACCACCCGGTGGCGATGATTGaaaaa               >  2:100968/1‑76 (MQ=255)
                                                             aggcaggGGCAGGGGGCCACCGGCCTCTCTTCCCCCGCCCAAAACACCAAACACCTGGGGGCGATGCTTGaaaaaa              >  2:6314834/1‑76 (MQ=255)
                                                                ggggggcgggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCa           <  2:468180/74‑1 (MQ=255)
                                                                caggggcaggggGGCACCGTCCTCTTTACCCCCCCCCAAATCACCCACCACCCGGTGGCGCCGATTGAAAAAACCa           >  2:5723349/1‑76 (MQ=255)
                                                                caggggcaggggGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCa           >  1:2422706/1‑76 (MQ=255)
                                                                aagggggaggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCa           <  1:3894664/75‑1 (MQ=255)
                                                                  ggggccggggccCCCCTTCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt         <  1:1176265/76‑1 (MQ=255)
                                                                  ggggcaggTGGTCACCCGCCTCTCTGCCCCCGCCCAAATAACCCACCACCTGGTGGCGATGATTGAAAAAACCAtt         >  2:4274639/1‑76 (MQ=255)
                                                                  ggggcaggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGCGGCGATGATTGAAAAAACCAtt         >  2:7704570/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  1:4557662/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  2:336392/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  2:2081932/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCga    >  2:5298665/1‑75 (MQ=255)
                                                                           gTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGcc  <  2:1396570/76‑1 (MQ=255)
                                                                            |                                                                             
TACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCC  >  REL606/359‑507

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