Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,332,8230AG31.7% 74.4 / 12.9 42pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/26);  new base G (11/2);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 1.69e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCGAG  >  NC_000913/2332604‑2333043
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
aaCCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAAAAAATAGCGGCTCGGCTTATCGGCGGCCGGGGGGAtttcat                                                                                                                                                                                                                                >  2:1746470/1‑226 (MQ=255)
            gTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGc                                                                                                                                                                                                                  <  1:1337797/229‑1 (MQ=255)
            gTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGTTTAAAGCcacctaaacgc                                                                                                                                                                                                                   >  1:1355782/1‑221 (MQ=255)
                          ttttGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGATTTACGCCCACATAACCCCGGGGGCCGGaac                                                                                                                                                                                                     >  1:560857/1‑226 (MQ=255)
                               aaGATTCCAGCGGGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTCCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGttt                                                                                                                                                                                               <  1:1553681/228‑1 (MQ=255)
                                                        ccaacccGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGaca                                                                                                                                                                            <  1:1835337/216‑1 (MQ=255)
                                                        cacaACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTTTCGGGCGGGGGGGGGATTTAACCCCCCATAACACCTGGGGCGGGACCCCGTTTGCCTGCCTTCCAGCGCAACGGGGa                                                                                                                                                                     >  2:1832215/1‑228 (MQ=255)
                                                           aaCCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGGGTTTAACCCACACGACACCGTGCGGCCGGCACCACGTTGCTGGCTTCCATCGACCCgggggtaagt                                                                                                                                                                >  2:735400/1‑223 (MQ=255)
                                                           aaCCGACTCGCCGCTATGACTGTATTGTGCGGCAGTACTTAATGAGGAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGaa                                                                                                                                                                    <  1:736690/225‑1 (MQ=255)
                                                              cGACTCGCCGCTATTACTGTATTGTGCGTCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACAGGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACCCGGGGAAACGTc                                                                                                                                                               <  2:1633118/227‑1 (MQ=255)
                                                                                 tATTGTGCGGCAGTACTTAATGGGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAc                                                                                                                                         <  1:1082438/229‑1 (MQ=255)
                                                                                    tgtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGATATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGa                                                                                                                                      <  1:1831451/229‑1 (MQ=255)
                                                                                       gCGGCAGTACTTAAATAGTAATGCTCTAGACCGCGTTCAATGATGATCTCTTTGGTGGTGGTGCCGCGATGCGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCGCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCa                                                                                                                                   <  2:884591/229‑1 (MQ=255)
                                                                                           cAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGaa                                                                                                                                                                    <  1:1306761/193‑1 (MQ=255)
                                                                                           cAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGGTTATCGGCGGCGGGGGGGTTTTAAGCCACACGAACGCGGGCGGGGGAGACCGCGTTGCTGCCTCCCACCGAGACcgggga                                                                                                                                                                    >  2:1304506/1‑191 (MQ=255)
                                                                                                   aaTGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCg                                                                                                                                                                                                            <  1:1295406/147‑1 (MQ=255)
                                                                                                        gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTa                                                                                                                                           <  1:421297/205‑1 (MQ=255)
                                                                                                           aTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAACCGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcgcgc                                                                                                              <  1:681291/229‑1 (MQ=255)
                                                                                                               gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAAc                                                                                                          <  1:818997/229‑1 (MQ=255)
                                                                                                                        cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTAACCCCCATTACCGCCCGGGGCGGGAACCCGTTTGCGGGCCTTCCAGGGACCAGGGGAAAGGCCCGGCATACACCAAAACGTTACAGGCAAATGAAAATTGTGGAGCGCGCAAAcatcacctt                                                                                                  >  2:1056572/1‑221 (MQ=255)
                                                                                                                        cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcgc                                                                                                                <  1:1058436/214‑1 (MQ=255)
                                                                                                                                   gATCGCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCggg                                                                                          <  1:1455112/224‑1 (MQ=255)
                                                                                                                                         tttGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGGGGGGGGGATTTACAGCCCAATACCGCTGGCCGCGGGAAACGGTTTGCCGGCTTTCCGGCGCACGGGGGAAACGCCGGGACAAAGCAAAACGCTATCCAACAATGGATAATAGCTGGGGGCCAAACGTAAAACACACGCGGGCGTTTTTTa                                                                                 >  2:1441955/1‑228 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCAt                                                                    <  1:472287/229‑1 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAATTAAAGCCACCTGACCCCCCGCGGCCGGCACCCGTTCGCTTCCTTCCCGCGACCCCGGGGAAACGCCCGCACAAACCCAATCGGTACCGACCCTTGCGAATTGCTGTGCGCCCCAACCTAACCTCACCCGGGCCTTTTTTTCCGGGTCGccc                                                                     >  2:471458/1‑227 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAGCGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCAt                                                                    <  1:675513/229‑1 (MQ=255)
                                                                                                                                                   gtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCAt                                                                    <  2:231386/229‑1 (MQ=255)
                                                                                                                                                               gccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAgc                                                        <  1:1139279/228‑1 (MQ=255)
                                                                                                                                                                                  agaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTTACGCCACATGACCGCTGGCGTCGGaca                                                                                                                                                                                                     >  1:586388/1‑75 (MQ=255)
                                                                                                                                                                                  agaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGaca                                                                                                                                                                                                     <  2:585364/75‑1 (MQ=255)
                                                                                                                                                                                           aCAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGc                            <  1:387897/228‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAACGCCCCCTGACCGCCGGCGTCGGGCACCGTTTTGCTGCCTTTCAGCGACACGGGGGAACGTCCGGCCTAACGCCAATCGTTTCCCCCCCTTGATAAGTGCTGTGCGCGCAAACTTAACTGCACCCGGGCATTTTTTTCCGGCTCGCCCTCCCGGGCGAGCGGTTTCCGTTTGCCGCTTAACCCCTCGcccttgttt                    >  2:271115/1‑222 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATCGGCGGCGGGGGGATTTAAAGCCACATTACCGc                                                                                                                                                                                                                  >  1:531797/1‑47 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGc                                                                                                                                                                                                                  <  2:530854/47‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCCCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAcc                                                                                                                                        <  1:962276/117‑1 (MQ=255)
                                                                                                                                                                                                  gCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAACGCTGGGGGCGGGCACCGTTTCGCTGCCCTTCCGCGACACGGGGAAAAGTCCGGGATAACGGAAATCGTTACCGCCCAATGATAATTTGCTAGCGCGCAAACTTTAACGCCCCCGGGCCTTTTTTTACGGCGCGCCCCTCCGGGGTAGGAGGTGTCAGTTGTCCGGTTACCGCTTCCCCAGtttt                    >  2:699247/1‑229 (MQ=255)
                                                                                                                                                                                                   gggCTCGGTTTATCGGCGGCGGGGGAATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCTCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCa                          <  2:967804/221‑1 (MQ=255)
                                                                                                                                                                                                      cTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGTCCGCTGGCGTCGGACACCGTCTCGCTGCCTTCCATCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAtt                                                                                                                                 <  1:1187920/120‑1 (MQ=255)
                                                                                                                                                                                                              aTCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGaa                 <  1:686231/219‑1 (MQ=255)
                                                                                                                                                                                                                  ggcggcggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGc     <  1:1749447/228‑1 (MQ=255)
                                                                                                                                                                                                                   gcggcggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCg    <  2:559873/228‑1 (MQ=255)
                                                                                                                                                                                                                       gcggGGTGATTTAACGCCACATGACCGCTTGCGTCGGGCACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAAAGCAAATCGTTACCGACCATTGATTATTGCTGAGCGCCCAAAACTAACTGCACGCGGGCATTTTTTACCGGGCCGCCCACCGGGGCGGGCAGGTGCAGTTCTCCGCCGACCCCCTCCCCACGTTTTGAACCGTTTTTTgtgcgg  >  2:1390557/1‑223 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
AACCATTCATACGTAACAGGAACCTGTTTTGAAGATTCCAGCGCGGCATAACGGAACACAACCGACTCGCCGCTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACCGGCTCGCCATCCGGGTAGAGCAGTTGCAGTTTGCCGCTGACCGCTTCGCCAGTTTTGAACTCTTTTTTGGCGAG  >  NC_000913/2332604‑2333043

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