Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,360,769 G→A intergenic (+153/‑407) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,7690GA100.0% 228.0 / NA 70intergenic (+153/‑407)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/37);  total (33/37)

CCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAAT  >  NC_000913/3360637‑3360902
                                                                                                                                    |                                                                                                                                     
ccAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGaa                                                                                                                                 <  1:744467/139‑1 (MQ=255)
        tgtCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCa                                                                                                                         >  2:263184/1‑139 (MQ=255)
        tgtCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAACTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAGCAAGTCTGATAAGTGAAGCACATCa                                                                                                                         <  2:54525/139‑1 (MQ=255)
                gAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAAttt                                                                                                                 <  1:634258/139‑1 (MQ=255)
                  aTCAAATTCCTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATGAAGCAATTtct                                                                                                               >  1:854813/1‑139 (MQ=255)
                  aTCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTtct                                                                                                               >  1:272481/1‑139 (MQ=255)
                           cTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACaa                                                                                                      >  2:116139/1‑139 (MQ=255)
                            tGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAAt                                                                                                     <  2:743431/139‑1 (MQ=255)
                              aCGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                                   <  1:116139/139‑1 (MQ=255)
                               cGAAGTGCGCGTTGTTTTCTGGCAAATGCGGAATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCg                                                                                                  <  2:272481/139‑1 (MQ=255)
                                 aaGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGcaca                                                                                                                            <  1:968574/111‑1 (MQ=255)
                                 aaGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGcaca                                                                                                                            >  2:968574/1‑111 (MQ=255)
                                   gTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGtctc                                                                                              >  1:708452/1‑139 (MQ=255)
                                                tCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGTATAATTTTCAATAAATTGTGCTTTGTTCACCTAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                                 <  2:854813/139‑1 (MQ=255)
                                                tCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                                 <  1:263184/139‑1 (MQ=255)
                                                tCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                                 <  1:75924/139‑1 (MQ=255)
                                                    gCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                             >  1:291281/1‑139 (MQ=255)
                                                       aaTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAg                                                                                                                                >  1:301948/1‑85 (MQ=255)
                                                       aaTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAg                                                                                                                                <  2:301948/85‑1 (MQ=255)
                                                       aaTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACctct                                                                          >  2:325942/1‑139 (MQ=255)
                                                              cATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATa                                                                   >  1:1214093/1‑139 (MQ=255)
                                                                  aGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                               <  2:708452/139‑1 (MQ=255)
                                                                       ttGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTg                                                          <  2:554984/139‑1 (MQ=255)
                                                                          aCCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCaa                                                       <  2:1214093/139‑1 (MQ=255)
                                                                                 cTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTc                                                >  2:583336/1‑139 (MQ=255)
                                                                                 cTATAATAGAGATAATTTTCAATAAATTATGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAGTTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTACTc                                                <  2:892312/139‑1 (MQ=255)
                                                                                         gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAGGCTATTCAGTCtttt                                        >  1:880141/1‑139 (MQ=255)
                                                                                         gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCtttt                                        >  2:828776/1‑139 (MQ=255)
                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTTACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                       >  1:604255/1‑139 (MQ=255)
                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCg                                                                <  2:489289/114‑1 (MQ=255)
                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCg                                                                >  1:489289/1‑114 (MQ=255)
                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                       >  1:1094828/1‑139 (MQ=255)
                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                       >  1:1005/1‑139 (MQ=255)
                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                                 <  2:95724/96‑1 (MQ=255)
                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                                 >  1:95724/1‑96 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCAGCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGt                                             <  2:839402/131‑1 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAAt                                                                                                     >  1:1073439/1‑75 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAAt                                                                                                     <  2:1073439/75‑1 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                                   >  1:317313/1‑77 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                                   <  2:317313/77‑1 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                              <  1:618507/130‑1 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                              >  2:618507/1‑130 (MQ=255)
                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGt                                             >  1:839402/1‑131 (MQ=255)
                                                                                                 tttCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                                <  1:583336/139‑1 (MQ=255)
                                                                                                      aTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTaaa                           <  2:291281/139‑1 (MQ=255)
                                                                                                      aTAAATTGTGCTTTGTTCACCAAGTCCGATAAGTGAATCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTaaa                           <  2:836340/139‑1 (MQ=255)
                                                                                                          aTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGt                       <  2:1005/139‑1 (MQ=255)
                                                                                                          aTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGt                       <  1:325942/139‑1 (MQ=255)
                                                                                                            tgtgCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTTAATTACTTGTTAACACCTCTGACAATACCGATTCTACAAGCTATTCAGTCTTTTCATTCTATCTAAATTGtct                     <  2:880141/139‑1 (MQ=255)
                                                                                                            tgtgCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGtct                     <  1:828776/139‑1 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                      <  1:1104846/99‑1 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                      >  2:1104846/1‑99 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTAtt                                                 >  1:1109761/1‑104 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTAtt                                                 <  2:1109761/104‑1 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGa                  <  1:302655/135‑1 (MQ=255)
                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGa                  >  2:302655/1‑135 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAAt                                                                                                                   >  2:531963/1‑36 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAAt                                                                                                                   <  1:531963/36‑1 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCa                                                        <  2:87876/95‑1 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCa                                                        >  1:87876/1‑95 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGatat               >  1:29378/1‑136 (MQ=255)
                                                                                                                     ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGatat               <  2:29378/136‑1 (MQ=255)
                                                                                                                        aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                             <  2:281386/71‑1 (MQ=255)
                                                                                                                        aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                             >  1:281386/1‑71 (MQ=255)
                                                                                                                        aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAatat         >  1:340167/1‑139 (MQ=255)
                                                                                                                               cTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                              <  2:320147/95‑1 (MQ=255)
                                                                                                                               cTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                              >  1:320147/1‑95 (MQ=255)
                                                                                                                               cTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAAtata        <  1:683704/133‑1 (MQ=255)
                                                                                                                               cTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAAtata        >  2:683704/1‑133 (MQ=255)
                                                                                                                               cTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAat  >  1:554442/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATACAAAAT  >  NC_000913/3360637‑3360902

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