Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG52.4% 55.8 / 120.0 82V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (10/29);  new base G (43/0);  total (53/29)
Fisher's exact test for biased strand distribution p-value = 5.05e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-10
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 183 total aligned reads displayed.

ACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCA  >  REL606/360‑508
                                                                         |                                                                                  
aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg                                                                                  <  2:6016971/76‑1 (MQ=255)
  aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggggc                                                                                >  1:7570431/1‑74 (MQ=255)
  aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                                <  1:7419243/76‑1 (MQ=255)
  aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                                <  1:5763543/76‑1 (MQ=255)
   tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                               <  1:5340886/76‑1 (MQ=255)
   tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                               <  2:7625629/76‑1 (MQ=255)
   tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                               <  2:7669612/76‑1 (MQ=255)
     aGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAcc                                                                             >  2:2680072/1‑76 (MQ=255)
      gTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCg                                                                           <  1:6683869/76‑1 (MQ=255)
        ggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCTCCGTc                                                                         >  2:4801847/1‑76 (MQ=255)
          cAAATGCAGAACGTTTTCTGGGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCTccc                                                                       >  1:783902/1‑75 (MQ=255)
          cAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCTccc                                                                       >  2:4315592/1‑75 (MQ=255)
           aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCctc                                                                      >  1:7221596/1‑76 (MQ=255)
            aaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctct                                                                    <  2:7330164/76‑1 (MQ=255)
             aTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctc                                                                   <  2:8568294/76‑1 (MQ=255)
              ttcAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGGACCGCGctctct                                                                  >  2:4686213/3‑76 (MQ=255)
              tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                                  <  1:2472173/76‑1 (MQ=255)
               gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                                 <  1:8707022/76‑1 (MQ=255)
               gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAAGGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                                 <  1:3500271/76‑1 (MQ=255)
                    aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                           <  2:5102326/76‑1 (MQ=255)
                     ctttttCTGCGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCg                                                          <  2:491520/74‑1 (MQ=255)
                      gTTTTCTGCGGGTTTCCGATATTCTGGAAAGCAATGCGAGGCAGGGGCAGGGGGGGCCCGTGCTTTCTGCCCCCGc                                                         >  2:3001008/1‑76 (MQ=255)
                      gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCACTCTCTCCCCCCGc                                                         >  2:3938116/1‑76 (MQ=255)
                      gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCGCCTCTCTGccccccc                                                         >  1:5919850/1‑74 (MQ=255)
                       ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCACCGGCCTCTCTTTCcccccc                                                        >  2:807692/1‑76 (MQ=255)
                       ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGcccccccc                                                        >  1:6483367/1‑76 (MQ=255)
                       ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGCCCCCCCTGCCCCCGcc                                                        >  2:7733077/1‑76 (MQ=255)
                          tCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCCCTCTCCCCCCCCCcaa                                                     >  1:6207131/1‑76 (MQ=255)
                          tCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCCGCCTCCCTGCCCCCCCCcaa                                                     >  1:7743526/1‑76 (MQ=255)
                            tGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCTCTGCCCCCCCCcaaaa                                                   >  2:4018884/1‑75 (MQ=255)
                              cGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGGCCTCTCTTCCCCCCCCcaaaaaa                                                 >  1:8439758/1‑73 (MQ=255)
                              cGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCACGTCCTCTCTGCCCCCGCCCAAATCa                                                 >  2:1819282/1‑76 (MQ=255)
                               gggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCACCCTCTCTGCCCCCCCCACAAACaa                                                >  2:642820/1‑75 (MQ=255)
                                ggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGTACCGTCCTCTCTTCCCCCCCCCAAAACaac                                               >  1:211154/1‑74 (MQ=255)
                                  ttcccGATATTCTGAAAAGCAAGGCCAGGCAGGGGCGGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                            <  1:5594972/73‑1 (MQ=255)
                                  ttGCCGATATTCGGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                            <  1:8042299/76‑1 (MQ=255)
                                   tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCCCCCCTCTTCCCCCGCCAAAAccaccacc                                           >  1:6384591/1‑74 (MQ=255)
                                   tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGGCCTCCCTGCCCCCGCCCAAAACaccaac                                           >  2:3571782/1‑76 (MQ=255)
                                     ccGTTATTCGGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaacca                                         <  1:3823626/76‑1 (MQ=255)
                                      cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCCAAAAcaaccaacac                                         >  2:307882/1‑74 (MQ=255)
                                        atatTCTGGAAAGCAATGCCGGGCAGGGGCAGGTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                     <  1:5970730/76‑1 (MQ=255)
                                        atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                     <  1:2906708/76‑1 (MQ=255)
                                        atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTGCCCCCGCCAAAAACACCCacccccc                                     >  2:3527923/1‑75 (MQ=255)
                                          atTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCCCTCTGCCCCCGCCAAAATCACCAACCCCCtgg                                   >  1:1795381/1‑76 (MQ=255)
                                          atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCTTTTGCCCCCCCAAAATAACCAAaccacccg                                    >  2:5673927/1‑74 (MQ=255)
                                          atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCCGCTCCCTGACCCCCCCAAAAACAACAAACACCcgg                                   >  2:3706137/1‑76 (MQ=255)
                                            tCTGAAAAGCAATGCCAGGCGGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtg                                 <  1:3057228/76‑1 (MQ=255)
                                             cTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCCCTCTGCCCCCCCCCAAAACACCAACCACCcggtgg                                >  1:4658224/1‑76 (MQ=255)
                                              tGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCCCCCAAATCACCAACCACCcggtggg                               >  2:4624061/1‑75 (MQ=255)
                                              tGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGCCCTCTCTGCCCCCCCCCAAAACACCAACCACCTCGGGGc                               >  1:325736/1‑76 (MQ=255)
                                              tGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCCCCTCTCTTCCCCCGCCCAAAACACCCACCACCCggcggc                               >  1:7509451/1‑76 (MQ=255)
                                              tGGAAAGAAAGCCCAGGCGGGGGCGGGGGCCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                               <  1:7788918/76‑1 (MQ=255)
                                              gggAAAGCAATGCCAGGCAGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                               <  1:6973305/75‑1 (MQ=255)
                                                gAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCTCTCCCCCCCCCAAAATCACCAACCACCAGGTGGCgt                             >  2:8749103/1‑75 (MQ=255)
                                                gAAAGCAATGCAGGGCGGGGGGAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCga                             <  1:444289/76‑1 (MQ=255)
                                                 aaaGCAATGCCAGGCAGGGGCAGGGGGCCACCCTCCTCTCTGCCCCCCCCCAAATCACCAACCACCTGGTGGCgaa                            >  1:8145090/1‑75 (MQ=255)
                                                  aaGCAATGCCAGGCAGGGGCAGGGGGGCGCCGTCCTCTCTTCCCCCGGCACAATCACCAACCACCTGGTGGCggtg                          >  2:1062384/1‑76 (MQ=255)
                                                  aaGCAATGCCAGGCAGGGGCAGGGGGCCCCCGTCCCCTCTGCCCCCCCCAAAAACACCAAACACCTGGGGGCgctg                          >  1:3997355/1‑76 (MQ=255)
                                                  aaGCAATGCCAGGCAGGGGCAGGGGGCCACCGGCCTCTCTGCCCCCGCCCAAAACACCAACCACCCGGGGGCgatg                          >  1:899642/1‑76 (MQ=255)
                                                   aGCAATGCCAGGCAGGGGCAGGGGGCCACCCTCCCCTCTCTCCCCGCCAAAAATACCCACCACCCCGGGGCgatca                          >  2:5798188/1‑74 (MQ=255)
                                                     cAATGCCAGGCGGGGGCAGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGAtt                      <  1:3564095/76‑1 (MQ=255)
                                                     cAATGCCAGGCAGGGGCAGGTGGCCACCCCTCTCTCTTGCCCCCCCAACAATAACAACCACCTGGTGGCGGTAAtt                      >  2:1645302/1‑76 (MQ=255)
                                                     cAATGCCAGGCAGGGGCAGGGGGGCACCCGCCGCTCTTCCCCCGCCAAAAAAACCAACCACCTGGGGGCggtggat                       >  2:2449449/1‑73 (MQ=255)
                                                     cAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCCCCAAAATCACCAACCACCTGGTGGCGATGAtt                      >  2:6528450/1‑76 (MQ=255)
                                                     cAATGCCAGGCAGGGGCAGGGGGCCACCGCCCCCTCTTCCCCCCCCCAAATCACCAACCACCTGGTGGCGATGAtt                      >  1:4975685/1‑76 (MQ=255)
                                                      aaTGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTTCCCCCCCCAAAATCACCAACCCCCCGGGGGGGATGAttt                     >  2:3503072/1‑75 (MQ=255)
                                                      aaTGCCAGGCAGGGGCAGGGGGCCCCCGGCCCTCCTTGCCCCCCCCAAAATACCCACCACCCTGGTGGGAAGATTg                     >  2:1481032/1‑76 (MQ=255)
                                                      aaTGCCAGGCAGGGGCAGGGGGCCCCCGGCCTCTCTGCCCCCCCCAAAATCACCAACCACCTGGGGGCGGTGATTg                     >  2:7999082/1‑76 (MQ=255)
                                                      aaTGCCAGGCAGGGGCAGGGGGCCCCCCTGCTCTTTTGCCCCCGCAAAAACAACAACCAACTGGGTGcgacgagtt                     >  2:419281/1‑73 (MQ=255)
                                                      aaTGCCAGGCAGGGGCAGGGGGCCACCCTCCTCTCTGCCCCCGCCAAAATCACCAACCCCCTGGTGGGGCAGATTg                     >  2:4262372/1‑76 (MQ=255)
                                                       aTGCCAGGCAGGGGCAGGGGGGCACCCGCCTCTCTTGCCCCCCCAAAAACACCAACCACCTGGGGGCGATGATTGa                    >  2:3431464/1‑76 (MQ=255)
                                                       aTGCCAGCCAGGGGCAGGTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                    <  1:5827945/76‑1 (MQ=255)
                                                        tGCCGGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaa                   <  1:5185591/76‑1 (MQ=255)
                                                        tGCCAGGCGGGGAATGGTGGCCACCGCCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaa                   <  1:1670153/76‑1 (MQ=255)
                                                        tGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTTCCCCCCCCAAAAACCACAACCACCTGGTGGCGGGGATTGaa                   >  2:130671/1‑76 (MQ=255)
                                                          ccAGGCAGGGGCAGGGGGGCACCCCCCCCTCTGCCCCCCCCCAAAACACCAAACACCCGGCGGCGCTGATTGaaaa                 >  2:6844654/1‑76 (MQ=255)
                                                          ccAGGCAGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaa                 <  1:8925750/76‑1 (MQ=255)
                                                          ccAGGCAGGGGCAGGGGGCCACCGGCCCCCTTCCCCCCCCCAAAATAACCAACCACCTGGTGGCGATGATTGaaaa                 >  2:6667911/1‑76 (MQ=255)
                                                           caggcaggGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAAACACCAACCACCTGGTTGCGATGATTGaaaaa                >  2:2236285/1‑76 (MQ=255)
                                                            aggcaggGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAACAACCAACCACCCGGTGGCGATGATTGaaaaaa               >  2:2864993/1‑76 (MQ=255)
                                                            aggcaggGGCAGGGGGGCACCGTCCCCTCTGCCCCCCGCCAAAACACCCACCACCCTGTGGCGGTGATTGtaaaaa               >  2:3037594/1‑76 (MQ=255)
                                                            aggcaggGGCAGGGGGCCCCCGTCCTCTCTTCCCCCGCCAAAATCACCAACCACCCGGTGGCGATGATTGaaaaaa               >  1:1089777/1‑76 (MQ=255)
                                                             ggcaggggcaggGGGGCCCCCTCCTCTCCGCCCCCGCCAAAATCACCCAACACCCGGGGGCGCGGCATTGaaaaac               >  2:6285960/1‑75 (MQ=255)
                                                              gcaggggcagggGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAcc             <  2:4427675/76‑1 (MQ=255)
                                                              gcaggggcagggGGCCCCCCTCCTCTCTCCCCCCGCCCAAATCACCAACCACCTCGGGGGGCTGGTTTAAAAAAca             >  2:4509737/1‑75 (MQ=255)
                                                               caggggcaggggGGCACCGTCCTCTCTGCCCCCCCCAAAATAACCAACCACCCGGTGGCGCTGGTTGAAAAAAACa            >  2:2987883/1‑76 (MQ=255)
                                                               caggggcaggggGGCACCGGCCTCTCTGCCCCCCCCAAAATCACCCACCCCCTGGTGGCGATGGTTTGAAAAAAca             >  2:2403578/1‑75 (MQ=255)
                                                               caggggcaggTGGGCACCGTCCTCTTTAGCCCCGCCAAAATAAACAACCAACTGGTGGCGAAGATTTTAAAAACCa            >  2:3896050/1‑76 (MQ=255)
                                                                aggggcgggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAt           <  1:8850531/76‑1 (MQ=255)
                                                                aggggcaggGGGGCCCCCTCCCCTCTTCCCCCCCCCAAAACAACAACCACCCGGTGGCGGTGATTGAAAAAAcaaa           >  2:6885823/1‑73 (MQ=255)
                                                                aggggcaggGGGGCAACGGCCCCCCTCTCCCCCCCAAAAAAACCAAACACCACGCGGCGGTGATTGAAAAAAACAt           >  2:455048/1‑76 (MQ=255)
                                                                 ggggcaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:6071927/76‑1 (MQ=255)
                                                                 ggggcaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:476513/76‑1 (MQ=255)
                                                                 ggggcaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt          <  1:2349412/76‑1 (MQ=255)
                                                                  gggccggGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa         <  2:8876078/76‑1 (MQ=255)
                                                                   gcgaggTGGCCCCCTCCCTCTCGCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAg        <  2:6808354/73‑1 (MQ=255)
                                                                     caggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg      >  2:205650/1‑76 (MQ=255)
                                                                     caggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg      >  1:1232575/1‑76 (MQ=255)
                                                                     agggTGCCCTCGGTCCTCTCCGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg      <  2:8181810/74‑1 (MQ=255)
                                                                      gggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg     <  2:325736/75‑1 (MQ=255)
                                                                         |                                                                                  
ACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCA  >  REL606/360‑508

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