Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,899,237 T→C E61G (GAG→GGG)  ygcU ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,899,2370TC100.0% 212.3 / NA 59E61G (GAG→GGG) ygcUputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/39);  total (20/39)

AATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTCC  >  NC_000913/2899092‑2899380
                                                                                                                                                 |                                                                                                                                               
aaTGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCtgt                                                                                                                                             <  2:230984/150‑1 (MQ=255)
 aTGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTgg                                                                                                                                           >  1:373151/1‑151 (MQ=255)
      ttgattCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAAcc                                                                                                                                       <  1:240393/150‑1 (MQ=255)
         attCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGt                                                                                                                                   >  2:786881/1‑151 (MQ=255)
         attCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGt                                                                                                                                   >  1:66228/1‑151 (MQ=255)
         attCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGt                                                                                                                                   >  2:222205/1‑151 (MQ=255)
         attCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGtt                                                                                                                                  <  2:337076/151‑1 (MQ=255)
                gCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGa                                                                                                                            <  1:590608/151‑1 (MQ=255)
                 cGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGAc                                                                                                                           <  1:9194/151‑1 (MQ=255)
                   gAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCg                                                                                                                         >  2:518843/1‑151 (MQ=255)
                        gTCGAGCCCCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTg                                                                                                                      >  1:62225/1‑149 (MQ=255)
                            aGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGaa                                                                                                                <  2:692944/151‑1 (MQ=255)
                            aGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGa                                                                                                                 >  1:363356/1‑150 (MQ=255)
                              caccacCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATc                                                                                                              >  1:712427/1‑151 (MQ=255)
                                ccaccGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCg                                                                                                             <  1:554083/150‑1 (MQ=255)
                                    cGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCaaa                                                                                                        <  1:534426/151‑1 (MQ=255)
                                         tttCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGtat                                                                                                    <  2:454036/150‑1 (MQ=255)
                                          ttCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGtata                                                                                                   <  1:40501/150‑1 (MQ=255)
                                               caacaGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGcc                                                                                             <  1:291978/151‑1 (MQ=255)
                                                  caGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATg                                                                                          <  1:824832/151‑1 (MQ=255)
                                                   aGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGa                                                                                         <  2:791290/151‑1 (MQ=255)
                                                    gTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGCTACTTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGGTTTACGACCGCTGCCGGAATCGCCAAAGTATAAATGCCATGa                                                                                         >  1:109695/1‑150 (MQ=255)
                                                       tCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAata                                                                                      <  2:70174/150‑1 (MQ=255)
                                                        ccAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATTCCATGAATATc                                                                                    <  1:794840/151‑1 (MQ=255)
                                                        ccAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATc                                                                                    <  1:58096/151‑1 (MQ=255)
                                                        ccAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATc                                                                                    >  1:799618/1‑151 (MQ=255)
                                                           gCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCgg                                                                                 >  1:134983/1‑151 (MQ=255)
                                                           gCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCgg                                                                                 <  1:206871/151‑1 (MQ=255)
                                                                   tCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGTAAAttt                                                                           <  1:794000/149‑1 (MQ=255)
                                                                   tCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTa                                                                          <  2:548441/150‑1 (MQ=255)
                                                                    cGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTAc                                                                         <  1:632449/150‑1 (MQ=255)
                                                                        ggCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGaaaa                                                                    <  1:504347/151‑1 (MQ=255)
                                                                         gCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAAc                                                                   <  2:256874/151‑1 (MQ=255)
                                                                            gAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGt                                                                >  2:588279/1‑151 (MQ=255)
                                                                            gAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGt                                                                <  1:436365/151‑1 (MQ=255)
                                                                            gAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGt                                                                <  1:368950/151‑1 (MQ=255)
                                                                             aaGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTc                                                               <  1:41151/151‑1 (MQ=255)
                                                                                 aCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAAt                                                            >  1:415259/1‑150 (MQ=255)
                                                                                   cGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACt                                                         <  2:655242/151‑1 (MQ=255)
                                                                                      tACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGtt                                                      >  1:589306/1‑151 (MQ=255)
                                                                                        cgcgGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGttt                                                     <  2:430492/150‑1 (MQ=255)
                                                                                        cgcgGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGttt                                                     >  1:430460/1‑150 (MQ=255)
                                                                                         gcgGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGttttt                                                   <  2:747253/151‑1 (MQ=255)
                                                                                            gCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCt                                                 <  2:471070/150‑1 (MQ=255)
                                                                                                acCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAAt                                            <  2:235908/151‑1 (MQ=255)
                                                                                                  cGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAAt                                            <  2:527402/149‑1 (MQ=255)
                                                                                                     tAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGt                                       >  2:73491/1‑151 (MQ=255)
                                                                                                       aTTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGttt                                     <  1:328807/151‑1 (MQ=255)
                                                                                                            gtgCGCATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCg                                <  1:724094/151‑1 (MQ=255)
                                                                                                              gcgcATTCATAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGt                              <  2:556794/151‑1 (MQ=255)
                                                                                                                       tAAAATTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTAcgcg                     >  1:357776/1‑151 (MQ=255)
                                                                                                                           aTTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATc                  >  1:104915/1‑150 (MQ=255)
                                                                                                                           aTTCAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCt                 <  1:518805/151‑1 (MQ=255)
                                                                                                                              cAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGc               <  2:91136/150‑1 (MQ=255)
                                                                                                                              cAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCa              <  2:614199/151‑1 (MQ=255)
                                                                                                                              cAGCACACGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCa              <  2:470814/151‑1 (MQ=255)
                                                                                                                                 cacaCGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCAcc            >  2:188460/1‑150 (MQ=255)
                                                                                                                                    acGGGATACTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACcaac         >  2:507387/1‑150 (MQ=255)
                                                                                                                                           aCTTGCCCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTcc  <  2:233979/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCATGAATATCCGGAAATTTACGAAAACGGTCAATACTGTTTTTCTTTAATACTGTTTCATCGGTAATTACGCGATCTGCACCAACAATTTCC  >  NC_000913/2899092‑2899380

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