Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG36.8% 110.2 / 28.3 68V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (10/33);  new base G (25/0);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 7.46e-11
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 133 total aligned reads displayed.

GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCC  >  REL606/358‑507
                                                                             |                                                                             
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  1:4627762/76‑1 (MQ=255)
 tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                              <  1:5360779/76‑1 (MQ=255)
  aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                             >  2:4609472/1‑76 (MQ=255)
   cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGggggg                                                                            >  1:4580942/1‑75 (MQ=255)
   cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                                                                            <  1:4447646/76‑1 (MQ=255)
    aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCagggggaggggggc                                                                           >  2:5195199/1‑74 (MQ=255)
    aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                           <  2:4240439/76‑1 (MQ=255)
    aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                           <  2:1034111/76‑1 (MQ=255)
    aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc                                                                           <  2:6012834/76‑1 (MQ=255)
     tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                          <  1:3500212/76‑1 (MQ=255)
      cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc                                                                         <  1:5148637/76‑1 (MQ=255)
        gTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCg                                                                       <  1:874329/76‑1 (MQ=255)
         tGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCACCgc                                                                      >  2:4700527/1‑75 (MQ=255)
          ggCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTc                                                                    <  2:631207/76‑1 (MQ=255)
           gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTcc                                                                   <  1:6075013/76‑1 (MQ=255)
             aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                                 <  2:5475784/76‑1 (MQ=255)
                tGCGGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  2:3470720/76‑1 (MQ=255)
                tGCAGAACTTTTTCTGCGGGTTGCTGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  2:1654755/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  2:2980699/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  1:1001913/76‑1 (MQ=255)
                tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                             <  1:2841132/76‑1 (MQ=255)
                 gCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTg                                                            <  2:5345964/76‑1 (MQ=255)
                  cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCTGCCCTCTGc                                                           >  1:334446/1‑76 (MQ=255)
                      aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                       <  1:3879939/76‑1 (MQ=255)
                      aCGGTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                       <  2:4951772/76‑1 (MQ=255)
                       cGTTTTCTGCGGTTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCg                                                      <  2:5391975/76‑1 (MQ=255)
                       cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCg                                                      <  1:3644638/76‑1 (MQ=255)
                         ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGcc                                                    <  1:3078399/76‑1 (MQ=255)
                         ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCCCTCTCTCcccccc                                                    >  2:5288353/1‑76 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGGCACCGCCCTTCTTGcccccccc                                                    >  1:1898113/1‑76 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCa                                                   <  1:6008187/76‑1 (MQ=255)
                          tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCCCTCCTccccccccc                                                   >  2:5336554/1‑75 (MQ=255)
                            tcggCGGTTTGCCGATATTCTGAAAAGCAATGCCAGGCGGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCaaa                                                 <  2:5643050/73‑1 (MQ=255)
                              tGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCCCCcaaaa                                               >  2:5436236/1‑75 (MQ=255)
                               gggggTTGGCCGATATTCTGGAAAGCAAGGCCAGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAAt                                               <  2:5803284/74‑1 (MQ=255)
                               gggggTTCCCGATATTCTGAAAAGCAAGGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATc                                              <  2:4415598/74‑1 (MQ=255)
                               gCGGTTTGCGGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGCCCTCTCTGCCCCCGCCAAAATc                                              <  2:1269394/76‑1 (MQ=255)
                                   gTTGCCGATATTCTGGAAAGCAATGCCGGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCacca                                          <  2:4418257/76‑1 (MQ=255)
                                    ttGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaa                                        <  1:2542965/76‑1 (MQ=255)
                                    ttGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCCCTCTTCCCCCCCCAAAAACaccaa                                        >  2:2265296/1‑76 (MQ=255)
                                     tGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCCCCTCTCTGCCCCCGCCAAAAAAaccaac                                       >  2:1134219/1‑76 (MQ=255)
                                       gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCaccaacc                                      <  1:4384926/76‑1 (MQ=255)
                                       gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCTCTTCCCCCGCCAAAAATaccaacc                                      >  2:3032383/1‑76 (MQ=255)
                                       gCCGATATTCTGGAAAGCAAGGCCAGGCAGGGGCAGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCaccaacc                                      <  2:5716829/76‑1 (MQ=255)
                                         cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGGCCTCTCTTCCCCCGCCCAAAACACCAaccac                                    >  2:5995718/1‑76 (MQ=255)
                                          gATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCAGCCTTTCCCCCCCCCGCAAAATCaccacaccac                                    >  1:3048487/1‑76 (MQ=255)
                                           atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCCAAAACACCCaacaccc                                  >  2:5193339/1‑75 (MQ=255)
                                           atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCCCTCTTCCCCCCCCAAaaacaccaaacaccc                                  >  1:1159401/1‑75 (MQ=255)
                                           atatTCTGAAAAGCAATGCCAGGCGGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                  <  1:2763774/76‑1 (MQ=255)
                                            tatTCGGGAAAGCAAGCCCGGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtg                                 <  2:3686258/76‑1 (MQ=255)
                                             atTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGGCCCCGTTCTCTCTGCCCCCGCCCAAAATACCAAACACCttg                                >  2:2735758/1‑74 (MQ=255)
                                             atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCCCCCCTCTCCCCCCGCCCAAATCACCCACCACCtgc                                >  1:5303628/1‑75 (MQ=255)
                                             atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGGCCTCTTCTACCCCCCCCAAACCACcacccaccccg                                >  1:2263171/1‑73 (MQ=255)
                                               tCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtg                              <  2:6176302/76‑1 (MQ=255)
                                               tCTGAAAAGCAATGCCAGGCAGGGGCAGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtg                              <  2:3762569/76‑1 (MQ=255)
                                                cTGGAAAGCAATGCCAGGCAGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtgg                             <  2:4642024/76‑1 (MQ=255)
                                                cTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGCCCTCCCTGCCCCCCCCAAAATCAACAACCCCCcggggg                             >  2:2420357/1‑76 (MQ=255)
                                                 tGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                            >  1:3833688/1‑76 (MQ=255)
                                                 ggaaaaGCAATGCGGGGGGGGGGCGGGTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                            <  1:3943621/73‑1 (MQ=255)
                                                 ggaaaaCCAATGCCGGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                            <  1:5656477/73‑1 (MQ=255)
                                                  ggAACGCAATGCCAGGCGGGGGCGGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                           <  1:563886/76‑1 (MQ=255)
                                                  ggAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCCTCTCTCTTCCCCCCCCAAAAATACCAACCCCCCGGTGGCg                           >  2:880065/1‑76 (MQ=255)
                                                    aaaGCAATGCCAGGCAGGGGCAGGTGGGCCCCGTCCCCTCTGCCCCCCCCACAAACACCCAACACCTCGGGGCgat                         >  1:2025456/1‑76 (MQ=255)
                                                    aaaGCAATGCCAGGCAGGGGCAGGGGGCCACCCTCCTCTCTGCCCCCGCCCCAAACAACAACCCCCCGCGGGCgat                         >  2:4901940/1‑76 (MQ=255)
                                                     aaccAATGCCAGGCAGGGGCAGGGGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                        <  2:1085300/73‑1 (MQ=255)
                                                       aaaaTGCCAGGCGGGGATAGGCGCCCTCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatgat                      <  1:6232480/74‑1 (MQ=255)
                                                        cAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCCCCCAAAACACCAAACACCTGGTGGCGGTGAtt                     >  2:694964/1‑76 (MQ=255)
                                                         aaTGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                    >  2:3109663/1‑76 (MQ=255)
                                                         aaTGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCAAAATCACCCACCCCCCAGCGGCGATGATTg                    >  2:6163545/1‑76 (MQ=255)
                                                          aTGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTTCCCCCCCCAAAAACACCAACCACCTGGTGGCGATGAttaa                   >  2:2935621/1‑74 (MQ=255)
                                                          aTGCCAGGAGGGGGCGGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGa                   <  1:5166167/76‑1 (MQ=255)
                                                           tGCCAGGCAGGGGCAGGGGGCCACCGCCCCCTCTGCCCCCCCCCAAAACAACCACCACCTGGGGGCGCTGATTGaa                  >  2:5596762/1‑76 (MQ=255)
                                                           cgCCGGGCGGGGGCAGGTGCCCCCCGTCCTCTCGCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaa                  <  1:4821401/75‑1 (MQ=255)
                                                            gCCGGGGAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaa                 <  1:1289605/76‑1 (MQ=255)
                                                             ccAGGCAGGGGGAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaa                <  1:482567/76‑1 (MQ=255)
                                                             ccAGGCAGGGGCAGGGGGGCACCCCCCTCTCTGCCCCCCCCCAAAACAACAACCACCTCGGGGGGATGATTGgaaa                >  2:2686815/1‑76 (MQ=255)
                                                             ccAGGCAGGGGCAGGGGGCCACCCCCCTCTCTGCCCCCGCCAAAAACAACAAACACCCGGGGGCGATGATTTaaaa                >  2:2710694/1‑76 (MQ=255)
                                                              caggcaggGGGAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaa               <  1:5293475/76‑1 (MQ=255)
                                                              caggcaggGGCAGGGGGGCCCCGTCCTCTCTGCCCCCGCCAAAATCACCCACCACCTGGTGGCGATGAGTGaaaaa               >  2:1515473/1‑76 (MQ=255)
                                                               gggcgggggttggTGCCCCCTGTCCCCCCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaaa              <  1:5031605/75‑1 (MQ=255)
                                                               aggcaggGGCAGGGGGGCCCCGGCCTCTCTTCCCCCCCCAAAATCACCCACCCACCGGTGGCGGTGAGTGaacaaa              >  2:212674/1‑76 (MQ=255)
                                                               aggcaggGGCAGGGGGCCACCGTCCTCTCCGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaaa              >  2:1557643/1‑76 (MQ=255)
                                                                 gagggggcgggTGCCCCCCGCCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAcc            <  1:116988/73‑1 (MQ=255)
                                                                  caggggcaggggGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCCGGTGGCGAGGATTGGAAAAAACa           >  2:4779790/1‑76 (MQ=255)
                                                                   gggggcgggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAt          <  1:866277/75‑1 (MQ=255)
                                                                   gggggcgggTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAt          <  2:6123609/75‑1 (MQ=255)
                                                                   aggggcaggGGGGCCCCGCCCTCCCCTCCCCCCCCCCAAACACCAACCACCTGGTGGCGAATAATGAAAAAAACac          >  2:431566/1‑75 (MQ=255)
                                                                   aggggcaggGGGGCAACCCCCCCTCCTCCCCCCCCCCAAACACCAACCACCTGGGGGGGATGATTGAAAAAAACaa          >  2:3577094/1‑75 (MQ=255)
                                                                    ggggcgggTGGCCCCGTCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt         <  2:881560/76‑1 (MQ=255)
                                                                    ggggcaggTGCCCCCTGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAtt         <  1:2059439/76‑1 (MQ=255)
                                                                     ggggaggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa        <  2:2048170/76‑1 (MQ=255)
                                                                     gggcaggTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa        <  2:106321/76‑1 (MQ=255)
                                                                      ggcgggGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAg       <  1:6137778/76‑1 (MQ=255)
                                                                       gcaggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGc      >  2:5500030/1‑76 (MQ=255)
                                                                        gcgggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGc      <  2:5419205/73‑1 (MQ=255)
                                                                         aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  2:1388206/1‑76 (MQ=255)
                                                                         aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  2:2314741/1‑76 (MQ=255)
                                                                         aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    >  2:4665242/1‑76 (MQ=255)
                                                                          gtgTGCCCACCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg    <  2:2164645/74‑1 (MQ=255)
                                                                           ggTGCCCCCCGCCCTCCCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGc   <  1:5772142/76‑1 (MQ=255)
                                                                             |                                                                             
GTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCC  >  REL606/358‑507

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