Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064330TG54.3% 72.0 / 130.0 81V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base T (11/26);  new base G (44/0);  total (55/26)
Fisher's exact test for biased strand distribution p-value = 7.55e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-13
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 156 total aligned reads displayed.

TACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGC  >  REL606/359‑506
                                                                            |                                                                           
tACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg                                                                            <  1:7354319/76‑1 (MQ=255)
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                           >  1:6057085/1‑76 (MQ=255)
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                           >  2:3744044/1‑76 (MQ=255)
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                           >  1:6654286/1‑76 (MQ=255)
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                           >  1:4251425/1‑76 (MQ=255)
 aCATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGgggg                                                                           >  1:3679139/1‑76 (MQ=255)
  cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGggggg                                                                          >  1:5984403/1‑75 (MQ=255)
  cATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc                                                                         <  1:3297727/76‑1 (MQ=255)
    tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                       <  2:1243347/76‑1 (MQ=255)
    tCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa                                                                       <  1:2876416/76‑1 (MQ=255)
     cAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc                                                                     <  2:7313214/76‑1 (MQ=255)
       gTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGccgccg                                                                   >  1:3732938/1‑76 (MQ=255)
        tGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGt                                                                  <  1:4519190/76‑1 (MQ=255)
        tGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAcacc                                                                  >  1:3402430/1‑73 (MQ=255)
          gCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTcc                                                                <  2:2775826/76‑1 (MQ=255)
           cAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGAccc                                                               >  2:1477759/1‑75 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCACCGcccc                                                               >  2:6564304/1‑75 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                              <  1:7611029/76‑1 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctc                                                              <  2:2087890/76‑1 (MQ=255)
            aaaTGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGccccc                                                              >  2:4360938/1‑74 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                           <  2:2833908/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctct                                                           <  1:1797098/76‑1 (MQ=255)
               tGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCGCctctcc                                                           >  1:1971179/1‑75 (MQ=255)
                 cAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGCCctctcttc                                                         >  1:1460280/1‑74 (MQ=255)
                  aGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCACCCTCTCCGcc                                                        >  1:6016605/1‑76 (MQ=255)
                     aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                     <  1:3831932/76‑1 (MQ=255)
                     aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                     <  2:2894190/76‑1 (MQ=255)
                     aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                     <  1:1381632/76‑1 (MQ=255)
                      cGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGTCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCg                                                    <  2:3645051/76‑1 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGcc                                                  <  2:3776374/76‑1 (MQ=255)
                        ttttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCCCCCCCTCTTcccccccc                                                  >  1:6192139/1‑76 (MQ=255)
                         tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCAACCGCCCTTCTGTCccccccc                                                 >  2:850672/1‑75 (MQ=255)
                         tttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTGCCCCCCCCa                                                 >  2:6967077/1‑76 (MQ=255)
                          ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCCCTCTGccccccccca                                                >  2:891033/1‑74 (MQ=255)
                          ttCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTTCTCTCTCCCCCCCCCaa                                                >  2:6497502/1‑76 (MQ=255)
                          tccTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCaa                                                <  2:5663627/74‑1 (MQ=255)
                            cTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCaaca                                              >  2:1485004/1‑74 (MQ=255)
                              gCGGGTTGCTGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATc                                            <  2:7593796/76‑1 (MQ=255)
                              gCGGGTTGCCGATATTCGGAAAAGCAATGCCAGGCGGGGGCAGGGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATc                                            <  2:3088152/76‑1 (MQ=255)
                               ggggTTGCCGATATTCGGAAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCa                                           <  1:4412945/75‑1 (MQ=255)
                               cGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCCACGTCCCCTCTGCCCCCGCCcaaaacc                                           >  2:5479867/1‑73 (MQ=255)
                                  tttGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGCCCTCTCTGCCCCCGCCAAAATCacca                                        <  1:1934713/75‑1 (MQ=255)
                                  gTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCacca                                        <  1:4841890/76‑1 (MQ=255)
                                  gTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCTCTCTTCCCCCCCCCAAAACaccc                                        >  2:5071977/1‑75 (MQ=255)
                                  gTTGCCGATATTCTGAAAAGCAATGCCAGGAAGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCacca                                        <  1:5438588/76‑1 (MQ=255)
                                     gCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGCGCCCTTTGCCCCCGCCCAAAACaccaaac                                     >  1:1146949/1‑74 (MQ=255)
                                       cGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAACGGGCTCACTGCCTCCCTCACAAATACCGaccac                                   >  2:3408410/1‑76 (MQ=255)
                                       cGATATTCGGAAAACCAAGCCCAGGCGGGGGCAGGGGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAaccac                                   <  2:5360435/76‑1 (MQ=255)
                                         atatTCTGGAAAGCAATGCCGGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                 <  1:1339716/76‑1 (MQ=255)
                                         atatTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCCGGCCCTCCTCCCCCCGCCCAAATACCAacccccc                                 >  2:6034849/1‑75 (MQ=255)
                                         atatTCGGAAAAGCAATGCCGGGCGGGGGCGGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCt                                 <  1:2743788/76‑1 (MQ=255)
                                           tttCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtgg                               <  1:7063367/75‑1 (MQ=255)
                                           atTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCCCCGGCCTCTCTGCCCCCGCCAAAATCACCCACCACCtcg                               >  1:5549830/1‑74 (MQ=255)
                                            ttCTGAAAAGCAATGCCAGGAAGGGGGAGGTGCCCACCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggt                              <  2:3552434/76‑1 (MQ=255)
                                             tCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCCTCTCTGCCCCCCCCCCAATCAACAACCACCcgggg                             >  1:109052/1‑74 (MQ=255)
                                             tCTGGAAAGCAAGGCCAGGCGGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtggtg                             <  1:3070515/76‑1 (MQ=255)
                                              cTGGAAAGCAATGCCAGGCAGGGGCAGGTGGTCAACGACCTCGCTACCACCTCCCAAATCAACAATCACCttgtgg                            >  1:4026921/1‑76 (MQ=255)
                                              cTGGAAAGCAATGCCAGGCAGGGGCAGGGGGACCCCCGGCCCTCCGTGCCCGCCCAAAAAACCAACCACCaggggg                            >  2:4362020/1‑76 (MQ=255)
                                               tGGAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCCTCCCCCCCCAAAAACACCAACCACCCGGTGGc                           >  2:1323770/1‑76 (MQ=255)
                                               ggaaaaGCAATGCCAGGCGGGGGCGGGTGGCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                           <  2:1230573/73‑1 (MQ=255)
                                                ggAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCCCCCAAAACAACAAACACCCGGTGGCg                          >  2:6925977/1‑76 (MQ=255)
                                                ggAAAGCAATGCCAGGCAGGGGCAGGGGGGCACCCGCCTCTCTGCCCCCCCCAAAAACACCAACCCCCCGCTGGCg                          >  2:2368299/1‑76 (MQ=255)
                                                ggAAAGAAAGCCCAGGAAGGGGCGGGTGGCCACCGCCCTCGCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                          <  2:3673891/76‑1 (MQ=255)
                                                 gAAAGCCAATGCCAGGCGGGGGAGGGTGCCCCCCGGCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                          <  1:4721404/76‑1 (MQ=255)
                                                 gAAAGCAATGCCAGGCAGGGGCAGGGGGCCCACGTCCCCTCTTTCCCCCCCAAAAACACCAAACACCCGGGGGCga                         >  2:5393090/1‑76 (MQ=255)
                                                 gAAAGCAAGGCGTGGCAGGGGCAGGGGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCga                         <  1:6677731/76‑1 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGGGGGCCCCGTCCCCTCTGCCCCCGCCCAAAACAACAACCACCCGGGGGCgat                        >  2:1112616/1‑76 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGGGGGCACCGTCCTCTCTGCCCCCGCCCAAAACACCCACCACCTGGTGGCgat                        >  2:3458444/1‑76 (MQ=255)
                                                  aaaGCAATGCCAGGCAGGGGCAGGGGGGCACCCTCCCCTCTGCCCCCCCCCAAATCACCCAACACCCGGTGGCggt                        >  2:4273105/1‑74 (MQ=255)
                                                   aggCAATGCCAGGGAGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCgatg                       <  1:2298978/74‑1 (MQ=255)
                                                       cAATGCCAGGCAGGGGCAGGGGGCCCCCGTCCCCTCTTCCCCCGCCAAAATCACCAACCACCTGGTGGCGCTGAtt                    >  2:1875777/1‑76 (MQ=255)
                                                        aaTGCCAGGCGGGGGCGGGTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                   <  2:4718152/76‑1 (MQ=255)
                                                        aaTGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                   <  2:1372616/76‑1 (MQ=255)
                                                        aaTGCCAGGCAGGGGCAGGGGGGCCCCCGCCCCTCTGCCCCCCCCACAAACACCAACCAACTTGGGGGGGTGATTg                   >  2:3976906/1‑76 (MQ=255)
                                                        aaTGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTg                   >  2:589890/1‑76 (MQ=255)
                                                          tGCCAGGCAGGGGCAGGGGGCCCCCCACCCCTCTCCCCCCCCCAAAAACACCCACCACCTGGTGGCGATGATTGaa                 >  2:4052780/1‑76 (MQ=255)
                                                           gCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAAACACCAACCACCTGGTGGCGGTGATTGaaa                >  2:7263598/1‑76 (MQ=255)
                                                           gCCAGGCAGGGGCAGGGGGACACCGTTCTCTCTGCCCCCCCCAAAATCACCAACCACCTGGTGGCGATGTATGaaa                >  2:7636450/1‑76 (MQ=255)
                                                            ccAGGGAGGGGGGGGGCGCCCCGGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaa               <  1:7563518/76‑1 (MQ=255)
                                                             caggcaggGGCAGGGGGCCACCCTCCTCTCTTCCCCCCCCCAAATCACCCACCCCCCGGTGGCGATGATTGaaaaa              >  2:5028078/1‑76 (MQ=255)
                                                              gggcgggggcggtTGCCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaaaaa             <  1:3094292/75‑1 (MQ=255)
                                                              aggcaggGGCAGGGGGGCACCCTCCTCTCTGCCCCCCCCCAAAACACCAACCACCTGGTGGGGATGATTGaaaaaa             >  2:6856691/1‑76 (MQ=255)
                                                                gcaggggcgggTGTCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAAcc           <  1:4969142/76‑1 (MQ=255)
                                                                  gggggcaggTGCCCCCCGTCCTCTCTCCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAt         <  1:6450451/75‑1 (MQ=255)
                                                                  aggggcaggGGGGCCCCGTCCTCTCTGCCCCCCCCCAAAACACCAACCCCCCGGGGGCGGTGGTTTAAAAAAACAt         >  2:2558705/1‑76 (MQ=255)
                                                                   ggggcaggTGGCCCCCCCCCTCTCTACCCCCGCCAAAATCACCAACCACCTGGTGGCCATGATTTTAAAAACCAtt        >  2:2261389/1‑76 (MQ=255)
                                                                   ggggcaggGGGGCCCCGTCCTCTCTGCCCCCCCCAAAATCACCAACCACCTGGTGGCGGTGATTGAAAAAACCAtt        >  2:6834846/1‑76 (MQ=255)
                                                                    ggggaggGGGCCCCCCTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       <  2:3183266/76‑1 (MQ=255)
                                                                    gggcgggTGGCCACCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       <  2:7114248/76‑1 (MQ=255)
                                                                    gggcgggGGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTa       <  1:1585607/76‑1 (MQ=255)
                                                                     ggcgggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAg      <  2:5665691/76‑1 (MQ=255)
                                                                      gcagtTGCCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGc     <  1:5521932/76‑1 (MQ=255)
                                                                       gaggTGCCCCCCGCCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg    <  2:4275850/75‑1 (MQ=255)
                                                                       caggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCg    >  1:3338881/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  2:5717691/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  2:5433312/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  2:4036539/1‑76 (MQ=255)
                                                                        aggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   >  1:1229275/1‑76 (MQ=255)
                                                                         gggTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCgg   <  2:801724/75‑1 (MQ=255)
                                                                          ggTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGc  <  2:1843873/76‑1 (MQ=255)
                                                                            |                                                                           
TACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGC  >  REL606/359‑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: