Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,068,157 T→C intergenic (‑77/+16) yggE ← / ← argO uncharacterized protein YggE/L‑arginine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,1570TC100.0% 224.9 / NA 60intergenic (‑77/+16)yggE/argOuncharacterized protein YggE/L‑arginine exporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/26);  total (34/26)

CGCTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGACG  >  NC_000913/3068024‑3068292
                                                                                                                                     |                                                                                                                                       
cgcTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATc                                                                                                                                   >  1:847724/1‑140 (MQ=255)
cgcTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATc                                                                                                                                   >  2:588634/1‑140 (MQ=255)
cgcTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATc                                                                                                                                   >  2:83534/1‑140 (MQ=255)
   tGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGa                                                                                                                                <  1:101095/140‑1 (MQ=255)
     ccATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATa                                                                                                                              >  2:617230/1‑140 (MQ=255)
        tCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAgcg                                                                                                                           <  1:915406/140‑1 (MQ=255)
        tCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAgcg                                                                                                                           <  1:812158/140‑1 (MQ=255)
          ccGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGcc                                                                                                                         <  2:329586/140‑1 (MQ=255)
               aaTACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACt                                                                                                                    >  2:639639/1‑140 (MQ=255)
               aaTACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACt                                                                                                                    >  2:957368/1‑140 (MQ=255)
                  aCCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGaa                                                                                                                 >  2:927603/1‑140 (MQ=255)
                    ccATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                                               >  2:660238/1‑140 (MQ=255)
                      aTTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAg                                                                                                             >  1:383843/1‑140 (MQ=255)
                             ccgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGc                                                                                                           <  2:690262/135‑1 (MQ=255)
                             ccgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGc                                                                                                           >  1:690262/1‑135 (MQ=255)
                               gccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCa                                                                                                    <  1:927603/140‑1 (MQ=255)
                                ccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCAt                                                                                                   <  1:466819/140‑1 (MQ=255)
                                   gggCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGTTTGTGCATGAg                                                                                                <  1:957368/140‑1 (MQ=255)
                                    ggCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGc                                                                                               >  1:344372/1‑140 (MQ=255)
                                               tGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGtctc                                                                                    >  1:496509/1‑140 (MQ=255)
                                                    ttCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGcc                                                                               <  2:383843/140‑1 (MQ=255)
                                                       aCGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGc                                                                            >  1:109524/1‑140 (MQ=255)
                                                       aCGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGc                                                                            >  1:828131/1‑140 (MQ=255)
                                                         gTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTg                                                                          <  2:109524/140‑1 (MQ=255)
                                                            aGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCaa                                                                       >  1:860748/1‑140 (MQ=255)
                                                                tctcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCa                                                                   <  2:245575/140‑1 (MQ=255)
                                                                  tcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCCAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAAt                                                                 >  2:917005/1‑140 (MQ=255)
                                                                   ccGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTg                                                                          >  1:104053/1‑130 (MQ=255)
                                                                   ccGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTg                                                                          <  2:104053/130‑1 (MQ=255)
                                                                   ccGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATa                                                                >  1:233573/1‑140 (MQ=255)
                                                                    cGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaa                                                               >  2:730777/1‑140 (MQ=255)
                                                                     gtgtTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              <  1:741347/140‑1 (MQ=255)
                                                                     gtgtTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              <  1:617230/140‑1 (MQ=255)
                                                                     gtgtTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAACTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                                              <  2:989870/140‑1 (MQ=255)
                                                                          aaTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACcaca                                                         >  1:559324/1‑140 (MQ=255)
                                                                           aTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACAt                                                        >  1:996644/1‑140 (MQ=255)
                                                                                cTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAacac                                                   >  2:1005050/1‑140 (MQ=255)
                                                                                cTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAacac                                                   <  2:281509/140‑1 (MQ=255)
                                                                                cTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCCAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAacac                                                   <  1:917005/140‑1 (MQ=255)
                                                                                 ttGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAAcaca                                                  >  2:781687/1‑140 (MQ=255)
                                                                                 ttGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAAcaca                                                  >  1:16298/1‑140 (MQ=255)
                                                                                  tGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACAt                                                 >  1:425127/1‑140 (MQ=255)
                                                                                   gCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATc                                                >  1:104415/1‑140 (MQ=255)
                                                                                   gCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATc                                                >  1:423746/1‑140 (MQ=255)
                                                                                        aaaTTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCcaca                                           <  2:425127/140‑1 (MQ=255)
                                                                                        aaaTTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCcaca                                           <  2:686704/140‑1 (MQ=255)
                                                                                          aTTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACaat                                         <  1:83534/140‑2 (MQ=255)
                                                                                              aGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGa                                     <  2:847724/140‑1 (MQ=255)
                                                                                               ggACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGAt                                    <  1:730777/140‑1 (MQ=255)
                                                                                                       ggAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATg                            <  1:588634/140‑1 (MQ=255)
                                                                                                          ggCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATgcgc                         <  1:781687/140‑1 (MQ=255)
                                                                                                             aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCtgt                      >  2:241997/1‑140 (MQ=255)
                                                                                                              aGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCtgtg                     >  2:599367/1‑140 (MQ=255)
                                                                                                              aGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCtgtg                     <  2:423746/140‑1 (MQ=255)
                                                                                                              aGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCtgtg                     >  1:326829/1‑140 (MQ=255)
                                                                                                              aGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCtgtg                     >  1:944915/1‑140 (MQ=255)
                                                                                                                cTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCt                   >  1:429533/1‑140 (MQ=255)
                                                                                                                cTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCt                   >  1:651637/1‑140 (MQ=255)
                                                                                                                         tctGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTg          >  1:610934/1‑140 (MQ=255)
                                                                                                                          ctGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTgc         >  2:961232/1‑140 (MQ=255)
                                                                                                                                 aTTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGAcg  <  2:610934/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CGCTGCCATCCCGCTAATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAATGCGCTGTGCTTTTGCCGTGCGCAGACG  >  NC_000913/3068024‑3068292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: