Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 892,022 A→T N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome892,0220AT96.8% 212.7 / ‑3.8 63N196I (AAT→ATT) malRtranscriptional activator of the Mal operon
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (31/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTTGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCCGCGCAT  >  BS168_reference_genome/891891‑892159
                                                                                                                                   |                                                                                                                                          
taatttGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTaaa                                                                                                                                      <  5:205867/138‑1 (MQ=255)
     tGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACa                                                                                                                                <  3:350299/138‑1 (MQ=255)
     tGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTCTCCGGAGAAATACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACa                                                                                                                                <  5:288625/138‑1 (MQ=255)
        tGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACaa                                                                                                                               <  2:328450/136‑1 (MQ=255)
        tGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACaa                                                                                                                               >  1:328450/1‑136 (MQ=255)
          cgcgATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATg                                                                                                                             <  8:341947/136‑1 (MQ=255)
                       aaGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGAc                                                                                                             <  2:105988/139‑1 (MQ=255)
                        aGTGAACGATCTGTTTATCATCTTTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTg                                                                                                               <  8:46811/136‑1 (MQ=255)
                           gAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGCCTCCg                                                                                                         >  3:260639/1‑139 (MQ=255)
                            aaCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCg                                                                                                         >  2:133015/1‑138 (MQ=255)
                             aCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTcc                                                                                                          <  7:20430/136‑1 (MQ=255)
                             aCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                       <  4:29347/139‑1 (MQ=255)
                                        atcatcATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATg                                                                                            >  6:208672/1‑139 (MQ=255)
                                           atcatTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCg                                                                                          <  5:282581/138‑1 (MQ=255)
                                               tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAg                                                                                     >  4:63363/1‑139 (MQ=255)
                                               tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAg                                                                                     >  4:363032/1‑139 (MQ=255)
                                                  ctctTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACaa                                                                                      <  8:196381/135‑1 (MQ=255)
                                                  ctctTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAg                                                                                     >  8:59154/1‑136 (MQ=255)
                                                     tttCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCaaa                                                                                 <  5:47502/137‑1 (MQ=255)
                                                       tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACt                                                                                                            >  1:12073/1‑108 (MQ=255)
                                                       tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACt                                                                                                            <  2:12073/108‑1 (MQ=255)
                                                       tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCtg                                                                             >  1:310853/1‑139 (MQ=255)
                                                       tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTACCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCt                                                                              <  6:56745/138‑1 (MQ=255)
                                                        ccgcctaaaCACCGCTGCTGATTCCTCAGGCAATAATCCTTTATGAAAAGGGAATTCCCTTAATTTCCATAACCATTTTAAAAAACAATGCCCCCCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAAcc                                                                               <  7:369368/129‑1 (MQ=255)
                                                                 cacCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAg                                                                   >  1:278059/1‑139 (MQ=255)
                                                                    cGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGa                                                                <  2:211338/139‑1 (MQ=255)
                                                                     ggAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGAt                                                               <  3:181604/139‑1 (MQ=255)
                                                                      gAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGa                                                                >  5:307553/1‑137 (MQ=255)
                                                                        gCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGa                                                            <  6:307553/139‑1 (MQ=255)
                                                                                 tCAGGCAAGAATGCTTTCTGCAAAGGGTATTCACTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCAGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAAtt                                                   >  4:200320/1‑139 (MQ=255)
                                                                                  cAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTg                                                  >  3:364818/1‑139 (MQ=255)
                                                                                   aGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGt                                                 <  6:103284/139‑1 (MQ=255)
                                                                                          aaTGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCtttt                                          >  6:251062/1‑139 (MQ=255)
                                                                                              cTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCa                                       <  6:44934/138‑1 (MQ=255)
                                                                                               tttCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGAc                                                                                                             <  2:300238/67‑1 (MQ=255)
                                                                                               tttCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGAc                                                                                                             >  1:300238/1‑67 (MQ=255)
                                                                                                          ggTATTCCCTTTATTTCCATTACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCGTTTTCAACGTTTTTTCTCg                          >  2:143815/1‑139 (MQ=255)
                                                                                                           gTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTtc                             >  8:146145/1‑135 (MQ=255)
                                                                                                               tCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGGACCGTTTTTTCTCGTTGGc                     >  3:203785/1‑139 (MQ=255)
                                                                                                                 ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg                    <  6:166335/138‑1 (MQ=255)
                                                                                                                 ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCATTATGCGACAAGCAAACCTGTGACATAAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg                    <  4:79747/138‑1 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGATCGGAAATTGTCGCTTTTGCAACGTTTTTTATAGTTTGCGAAg                 >  6:201465/1‑138 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCa                                       >  6:162551/1‑116 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCa                                       <  5:162551/116‑1 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  2:141721/1‑138 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  1:222132/1‑138 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  4:310865/1‑138 (MQ=255)
                                                                                                                    ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  5:92611/1‑138 (MQ=255)
                                                                                                                    ctATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAg                                                                   <  6:295734/87‑1 (MQ=255)
                                                                                                                    ctATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAg                                                                   >  5:295734/2‑88 (MQ=255)
                                                                                                                      aTTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCa               <  3:63363/138‑1 (MQ=255)
                                                                                                                      aTTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCAt              >  6:194523/1‑139 (MQ=255)
                                                                                                                       tttCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  7:83696/1‑135 (MQ=255)
                                                                                                                       tttCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  8:67893/1‑135 (MQ=255)
                                                                                                                       tttCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg                 >  7:253476/1‑135 (MQ=255)
                                                                                                                       tttCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTAGTTGGCGAAg                 >  7:33834/1‑135 (MQ=255)
                                                                                                                           cATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTTTGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc         <  2:278059/139‑1 (MQ=255)
                                                                                                                           cATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGtt          <  4:364818/138‑1 (MQ=255)
                                                                                                                           cATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc         <  2:222132/139‑1 (MQ=255)
                                                                                                                             ttACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCcg       >  6:4267/1‑139 (MQ=255)
                                                                                                                             ttACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCcg       <  5:201465/139‑1 (MQ=255)
                                                                                                                             ttACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCcg       <  1:143815/139‑1 (MQ=255)
                                                                                                                              tACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc         >  7:369199/1‑136 (MQ=255)
                                                                                                                                 ctaTTTAAAAAACCAATGGCCGGGCCCCGCTGCCTCCCCAAAAGCTTATGCGACAAGCAAACCTGTGACATTAAGCGATGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCCGCGCat  <  6:25909/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                          
TAATTTGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTCCGCGCAT  >  BS168_reference_genome/891891‑892159

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