Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG54.4% 52.3 / 100.5 68V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (6/25);  new base G (37/0);  total (43/25)
Fisher's exact test for biased strand distribution p-value = 2.78e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-12
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 121 total aligned reads displayed.

GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGC  >  REL606/358‑506
                                                                            |                                                                             
gTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt                                                                               <  1:1398066/76‑1 (MQ=255)
gTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt                                                                               <  1:3850728/76‑1 (MQ=255)
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  2:1176095/76‑1 (MQ=255)
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  2:1309635/76‑1 (MQ=255)
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  2:2937470/76‑1 (MQ=255)
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  1:2410189/76‑1 (MQ=255)
  aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                             >  2:3795832/1‑76 (MQ=255)
  aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGttg                                                                             >  1:2500041/1‑74 (MQ=255)
  aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg                                                                             <  2:3966830/76‑1 (MQ=255)
    aTCAGCGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATAATCTGGAGAGCAATGCCAGGCAGGGGCCGGTggtg                                                                           >  1:680500/1‑74 (MQ=255)
     tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                         <  2:36150/76‑1 (MQ=255)
      cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc                                                                       <  1:30126/76‑1 (MQ=255)
       aGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAcc                                                                      <  1:2769225/76‑1 (MQ=255)
        gTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCg                                                                     <  2:4988602/76‑1 (MQ=255)
         tGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCAcct                                                                     >  1:1149783/1‑75 (MQ=255)
          ggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTc                                                                   <  2:220871/76‑1 (MQ=255)
           gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTcc                                                                  <  2:1923071/76‑1 (MQ=255)
           gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTcc                                                                  <  1:4195112/76‑1 (MQ=255)
           gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTcc                                                                  >  2:3574863/1‑76 (MQ=255)
           gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCGCACCgcc                                                                   >  1:2728854/1‑74 (MQ=255)
             aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTcccc                                                                >  2:1885988/1‑74 (MQ=255)
             aaaTGCAGAACGTTTTCTGCGGGTTACCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCGCCctc                                                                >  2:2184725/1‑76 (MQ=255)
              aaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctct                                                               <  1:254305/76‑1 (MQ=255)
               atacAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctc                                                              <  2:453809/73‑1 (MQ=255)
                tGCAGAACGTTTTCTGCTGGTTGCCGATATTCTGGAATGCAATGCCAGGCAGGGGCAGGTGGCCACCGGCctctct                                                            <  2:5051930/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                            <  2:4218558/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                            <  1:4748327/76‑1 (MQ=255)
                 gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                           <  1:252990/76‑1 (MQ=255)
                  cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCGCCCCTCccctcttc                                                          >  2:585070/1‑74 (MQ=255)
                   aGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCCTCCCCGcc                                                         >  2:3208303/1‑76 (MQ=255)
                      actttttGTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                      <  2:4443808/73‑1 (MQ=255)
                      aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                      <  2:2013791/76‑1 (MQ=255)
                       cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGGACCCTCCTCTCTccccccc                                                     >  2:4118282/1‑75 (MQ=255)
                       cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGcccccc                                                     >  2:2614370/1‑75 (MQ=255)
                        gTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCGc                                                    >  2:4977898/1‑76 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTCCCCCCGccc                                                  >  1:239888/1‑75 (MQ=255)
                          tgtCTGGGTGTTGCAGATATTCTGGAAAGCAAGGCCAGCCAGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCa                                                 <  2:2376989/74‑1 (MQ=255)
                            tCTGCGGTTTCCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCaaa                                               <  1:3315160/76‑1 (MQ=255)
                            tCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGTCCCTCTCCCCCCCCCcaa                                               >  2:2826741/1‑76 (MQ=255)
                              tGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTTGCCCCCCCCaaaa                                              >  2:4588835/1‑76 (MQ=255)
                                cGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGGCCTCTCTTCCCCCGCCCAAACCa                                           >  1:1262367/1‑76 (MQ=255)
                                cGGGTTGCCGATATTCTGGAAAGAAAGGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCa                                           <  1:338567/76‑1 (MQ=255)
                                 gggTTGCATATATGCTGAAAACCAATGCACGGCCGGGGCGGGGGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCac                                          <  2:706371/76‑1 (MQ=255)
                                  ggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCTCTGCCCCCGCCCAAAACacc                                         >  1:3227987/1‑76 (MQ=255)
                                  ggTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCCCTCTTCCCCCCCCACAAACaac                                         >  1:4030814/1‑74 (MQ=255)
                                   gTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTGCCCCCCCCCCAAACacca                                        >  1:537416/1‑76 (MQ=255)
                                    ttGCCGATATTCGGAAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                       <  1:493329/76‑1 (MQ=255)
                                     tgtcGATGTGCTGGAAAGAAATGCCGGGCGGGGGCAGGTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCaccaac                                      <  2:3829895/73‑1 (MQ=255)
                                     tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCTTCTCTCTGCCCCCGCCCAAATTcccaac                                      >  2:764122/1‑76 (MQ=255)
                                     tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCTCCCCCCCGCCCCCGCCAAAaacaacaac                                      >  1:1834032/1‑76 (MQ=255)
                                     tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCGTCTCTCTTCCCCCCCCCAAAACaccaac                                      >  2:4948449/1‑76 (MQ=255)
                                      gcagATATTCTGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaacc                                     <  2:4488471/73‑1 (MQ=255)
                                      gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTGCTCTCCGCCCCCGCCCAAATCaccaaac                                     >  2:2294844/1‑74 (MQ=255)
                                      gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCTCTGCCCCCGCCCAAAACaccaacc                                     >  2:2290871/1‑76 (MQ=255)
                                      gCCAATATTCGGGAAAGCAAGGCCAGGCGGGGGCGGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaacc                                     <  2:1392991/76‑1 (MQ=255)
                                        cGATATTCTGAAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAaccac                                   <  1:3710184/76‑1 (MQ=255)
                                          atatTGTGAAAAGCAAGGCCAGGCGGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                 <  2:1036499/76‑1 (MQ=255)
                                          atatTCTGGAAAGCAATGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                 <  1:4555251/76‑1 (MQ=255)
                                          atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTGCCCCCCCCAAAATCACCAacccccc                                 >  1:4864627/1‑75 (MQ=255)
                                           tatTCTGAAAAGCAATCCCAGGCGGGGGCGGGTGGCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCtg                                <  1:2765005/76‑1 (MQ=255)
                                             ttCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCCCCGCCCCTCCGCCCCCGCCCAAAACGCCAACCACCtgcg                              >  2:4588735/1‑74 (MQ=255)
                                               cTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTTCCCCCGCCAAAAATAACAAACACCtcgtgg                            >  2:1723852/1‑76 (MQ=255)
                                                tgaaaaGCATTGCAGGGCGGGGGAGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                           <  2:3087279/73‑1 (MQ=255)
                                                tGGAAAGCAATGCCAGGCAGGGGCAGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                           <  1:545355/76‑1 (MQ=255)
                                                 ggAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                          <  2:944602/76‑1 (MQ=255)
                                                 ggAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCCAAATCACCAACAACCTGGTGGCg                          >  2:2697044/1‑76 (MQ=255)
                                                 gaaaaGCAATGCCAGGCGGGGGCAGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                          <  2:2182663/74‑1 (MQ=255)
                                                 gaaaaCCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                          <  2:2220031/74‑1 (MQ=255)
                                                  gAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCGCCTCTCTTCCCCCCCCCAAAACACCCACCACCCCGTGGCga                         >  1:4075791/1‑76 (MQ=255)
                                                  gAAAGCAATGCCAGGCAGGGGCAGGGGGGCCACGTCCCCTCTGCCCCCCCCCAAAACACCAACCCCCCGGTGGCgg                         >  2:1001959/1‑75 (MQ=255)
                                                   aaaGCAACGCCAGGAAGGGGCGGGTGGCCCCCTTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgat                        <  1:4471220/76‑1 (MQ=255)
                                                   aaaCCAATGCGCGGCGGGGGCAGGTGACCCCGGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgat                        <  1:3699176/76‑1 (MQ=255)
                                                    aaGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                       <  2:2410043/76‑1 (MQ=255)
                                                    aaGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                       >  2:4171638/1‑76 (MQ=255)
                                                      gCATTGCCAGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                     <  1:603043/76‑1 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCAGGGGGGCCCCCCCCCCTCTGCCCCCCCCACAAATACCAACCACCTCGTGGGgatggt                     >  2:1311210/1‑74 (MQ=255)
                                                      gCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTTCCCCCCCCAAAAACACCAACCACCTGGTGGCggtgat                     >  2:3934273/1‑76 (MQ=255)
                                                       cAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCTCTCTTCCCCCGCCAAAATTAACAACCCCCCCGTGGCGAGGAtt                    >  2:1133118/1‑76 (MQ=255)
                                                        aaggCCGGGCGGGGGCGGTTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                   <  1:3214970/73‑1 (MQ=255)
                                                        aaTGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCCCCCAAAACACCAACCACCTGGGGGCGGTGATTg                   >  2:2387901/1‑76 (MQ=255)
                                                          tGCCAGGCAGGGGCAGGGGGCCCCCCGCCTCCCTCCCCCCCCCAAAAAAACCAACCAACCGGTGGCGCTGATTgca                 >  2:3927829/1‑74 (MQ=255)
                                                           gCCAGGCAGGGGCAGGGGGGCCCCCTCCCCTCTGCCCCCCCCAAAAATAACCAACACCCCGGGGCGATGAGTGgaa                >  2:685683/1‑76 (MQ=255)
                                                            ccAGGCGGGGGCGGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaa               <  1:4372958/76‑1 (MQ=255)
                                                            ccAGGCAGGGGCAGGGGGCCACCGTCCTCTCTTCCCCCCCCCAAATCACCAACCACCTGGTGGCGATGATTGaaaa               >  2:3521935/1‑76 (MQ=255)
                                                             caggcaggGGGAGGTGCCCCCCCTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaa              <  1:4984286/76‑1 (MQ=255)
                                                             caggcaggGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAAACACCAAACACCTGGTGGCGAAGATTGaaaaa              >  2:1362823/1‑76 (MQ=255)
                                                             caggcaggGGCAGGGGGGCACCGGCCCCTCTTCCCCCCCCCAAAACAACCAACCCCCGGTGGCGATGATTGgaaaa              >  2:4197521/1‑76 (MQ=255)
                                                             caggcaggGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATAACCCACCACCCGGGGGGGAGGATTTgaaaa              >  2:2923693/1‑76 (MQ=255)
                                                             caggcaggGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAAATACCCAACACCTTGGGGCGATGTTTTaaaaa              >  2:1967916/1‑76 (MQ=255)
                                                              aggcaggGGCAGGGGGCCAGCGGCCGCTCTGGCCCCGCCAAAATAACCAACCACCCGGTCGCGATGAGTTcaaaaa             >  2:2532010/1‑76 (MQ=255)
                                                                gcgggggcgggGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAcc           <  1:747253/73‑1 (MQ=255)
                                                                    gggcgggTCGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       <  1:2262892/76‑1 (MQ=255)
                                                                    gggcaggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       >  2:1264427/1‑76 (MQ=255)
                                                                    gggcaggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       >  1:3938100/1‑76 (MQ=255)
                                                                     gggaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAg      <  2:2328819/73‑1 (MQ=255)
                                                                       caggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg    >  2:1605516/1‑76 (MQ=255)
                                                                        aggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   <  2:2249482/76‑1 (MQ=255)
                                                                        aggTGGCCACCGTCTTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  1:1312622/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCTCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  2:4168026/1‑76 (MQ=255)
                                                                          ggTGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGc  <  1:1281837/76‑1 (MQ=255)
                                                                            |                                                                             
GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGC  >  REL606/358‑506

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