Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,159,0760TA66.1% 83.6 / 63.6 62L224Q (CTG→CAG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (12/9);  new base A (21/20);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTGGTTTA  >  NC_000913_3_pae_tpiA/4158930‑4159220
                                                                                                                                                  |                                                                                                                                                
ggCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGa                                                                                                                                               >  3:385750/1‑150 (MQ=255)
        gcgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTgc                                                                                                                                       <  2:758579/150‑1 (MQ=255)
                  aCCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACAtggt                                                                                                                             >  2:967106/1‑150 (MQ=255)
                      aTGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTggcg                                                                                                                         >  1:537087/1‑150 (MQ=255)
                             aTCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGATCAGCCTGGa                                                                                                                                               <  1:682432/121‑1 (MQ=255)
                             aTCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGACGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGa                                                                                                                                               >  2:682432/1‑121 (MQ=255)
                                    aGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAGACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGAt                                                                                                           >  4:222146/1‑150 (MQ=255)
                                       ggAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                                                        >  2:768139/1‑150 (MQ=255)
                                        gAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGAGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCaa                                                                                                                                   <  4:150161/122‑1 (MQ=255)
                                        gAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGAGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCaa                                                                                                                                   >  3:150161/1‑122 (MQ=255)
                                          aaaaCTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACttt                                                                                                     >  1:919530/1‑150 (MQ=255)
                                               tgctgATGCTGGAAGATGGTCACTGTTTGCGCGAGCAGGCAATGGGGTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGc                                                                                                <  1:317652/150‑1 (MQ=255)
                                                        tGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                                                       >  1:66754/1‑150 (MQ=255)
                                                        tGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                                                       >  3:455258/1‑150 (MQ=255)
                                                             gATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCtgtg                                                                                  <  1:881816/150‑1 (MQ=255)
                                                               tGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctggcagtaaa                                                                                <  1:967106/150‑7 (MQ=255)
                                                                  tCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACAtggt                                                                                                                             >  1:595180/1‑102 (MQ=255)
                                                                  tCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACAtggt                                                                                                                             <  2:595180/102‑1 (MQ=255)
                                                                  tCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTgccgcc                                                                             >  1:887348/1‑150 (MQ=255)
                                                                   cACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                                                        <  3:423994/122‑1 (MQ=255)
                                                                   cACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCAc                                                                                                        >  4:423994/1‑122 (MQ=255)
                                                                       gTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgc                                                                        <  2:66754/150‑1 (MQ=255)
                                                                        tttGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgcg                                                                       >  1:479677/1‑150 (MQ=255)
                                                                        tttGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAgcg                                                                       <  1:768139/150‑1 (MQ=255)
                                                                         ttGCGCGATCAGGCAATGGGTTGCTGGTGTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGGGCCGCCGGAgcgc                                                                      <  2:94249/150‑1 (MQ=255)
                                                                           gcgcGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                                                       >  4:211659/1‑131 (MQ=255)
                                                                           gcgcGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCgctgg                                                                                       <  3:211659/131‑1 (MQ=255)
                                                                            cgcgATACGGAAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATATACACTTCCGCGCGACCAGTCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGACGCCGGAGCGCaaa                                                                   <  4:455258/150‑1 (MQ=255)
                                                                                aTCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGCGGCGCAAAcgcg                                                               >  2:791654/1‑150 (MQ=255)
                                                                                     gCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGAACAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATggg                                                          >  2:293122/1‑150 (MQ=255)
                                                                                      cAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATgggg                                                         <  2:919530/150‑1 (MQ=255)
                                                                                           ggTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGgttgtt                                                    <  1:964774/150‑1 (MQ=255)
                                                                                             tttCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTa                                                  >  3:54902/1‑150 (MQ=255)
                                                                                              ttCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAGCGCGATGGGGTTGTTTAt                                                 >  3:317819/1‑150 (MQ=255)
                                                                                                cTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCt                                               >  1:765650/1‑150 (MQ=255)
                                                                                                 tGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTg                                              <  4:385750/150‑1 (MQ=255)
                                                                                                   ttttGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGcc                                            <  2:479677/150‑1 (MQ=255)
                                                                                                     ttGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAActc                                                                                                                                          <  2:674822/54‑1 (MQ=255)
                                                                                                     ttGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCAGGAAActc                                                                                                                                          >  1:674822/1‑54 (MQ=255)
                                                                                                              gggCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                                                                                                      <  2:466442/81‑1 (MQ=255)
                                                                                                              gggCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACtt                                                                                                      >  1:466442/1‑81 (MQ=255)
                                                                                                               ggCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtggtg                                                                                                                            >  2:177299/1‑58 (MQ=255)
                                                                                                               ggCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtggtg                                                                                                                            <  1:177299/58‑1 (MQ=255)
                                                                                                                          aTACACACTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGAGTTGTTTATCTGCCGTGCATTAAGCCGGAACCAcgcc                     <  2:887348/150‑1 (MQ=255)
                                                                                                                            acacacTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAg                                                                         >  1:317335/1‑96 (MQ=255)
                                                                                                                            acacacTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAg                                                                         <  2:317335/96‑1 (MQ=255)
                                                                                                                            acacacTACCGCGCGGCCCGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCAcgccgc                   >  2:998217/1‑150 (MQ=255)
                                                                                                                              acacTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCAc                 >  4:514748/1‑150 (MQ=255)
                                                                                                                              acacTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCg                                           <  3:316031/124‑1 (MQ=255)
                                                                                                                              acacTTCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCg                                           >  4:316031/1‑124 (MQ=255)
                                                                                                                               cacTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACt                >  1:671416/1‑150 (MQ=255)
                                                                                                                                  ttCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTAtt             >  4:584321/1‑150 (MQ=255)
                                                                                                                                   tCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTg            <  2:537087/150‑1 (MQ=255)
                                                                                                                                   tCCGCGCGACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTg            >  4:479420/1‑150 (MQ=255)
                                                                                                                                      gcgcgACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCAcgc                      >  2:396687/1‑137 (MQ=255)
                                                                                                                                      gcgcgACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCAcgc                      <  1:396687/137‑1 (MQ=255)
                                                                                                                                        gcgACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTg       <  2:765650/150‑1 (MQ=255)
                                                                                                                                          gACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTGGt     <  2:476046/150‑1 (MQ=255)
                                                                                                                                          gACCAGCCAGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTGGt     >  1:476046/1‑150 (MQ=255)
                                                                                                                                             cAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTGGTTTa  <  3:584321/150‑1 (MQ=255)
                                                                                                                                               gCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCAtt                                     <  3:557695/113‑1 (MQ=255)
                                                                                                                                               gCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCAtt                                     >  4:557695/1‑113 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGCCAGCGCTGGCTGTGCCGCCGGAGCGCAAACGCGATGGGGTTGTTTATCTGCCGTGCATTAAGCCGGAACCACGCCGCACTATTGGCCTGGTTTA  >  NC_000913_3_pae_tpiA/4158930‑4159220

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