Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,757,688 T→A 8.8% intergenic (+31/+33) lpp → / ← ldtE murein lipoprotein/murein L,D‑transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,6880TA 8.8% 186.3 / 8.0 57intergenic (+31/+33)lpp/ldtEmurein lipoprotein/murein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base A (2/3);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GTTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGA  >  NC_000913/1757541‑1757832
                                                                                                                                                   |                                                                                                                                                
gTTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAtt                                                                                                                                                 >  2:308112/1‑149 (MQ=255)
        gCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCa                                                                                                                                         >  4:295129/1‑149 (MQ=255)
        gCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                         <  2:441388/149‑1 (MQ=255)
        gCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                         <  2:219134/149‑1 (MQ=255)
          tGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtt                                                                                                                                       >  3:24250/1‑149 (MQ=255)
              cAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   <  2:274998/149‑1 (MQ=255)
               aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttttttt                                                                                                                                  >  4:96093/1‑149 (MQ=255)
               aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttttttt                                                                                                                                  <  3:460002/149‑1 (MQ=255)
                acgacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   <  3:148002/147‑1 (MQ=255)
                acgacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   >  4:148002/1‑147 (MQ=255)
                  gacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGt                                                                                                                               >  1:358832/1‑149 (MQ=255)
                       gAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGcc                                                                                                                          >  3:187645/1‑149 (MQ=255)
                        aaCGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGc                                                                                                                           <  3:543094/147‑1 (MQ=255)
                        aaCGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGc                                                                                                                           >  4:543094/1‑147 (MQ=255)
                        aaCGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCg                                                                                                                         <  4:23425/149‑1 (MQ=255)
                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTAc                                                                                                                    >  4:324125/1‑149 (MQ=255)
                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTAc                                                                                                                    >  1:109698/1‑149 (MQ=255)
                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAttttt                                                                                                                                    >  3:425252/1‑133 (MQ=255)
                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAttttt                                                                                                                                    <  4:425252/133‑1 (MQ=255)
                               tGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCg                                                                                                                  <  4:277292/149‑1 (MQ=255)
                                      cGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGc                                                                                                           <  2:441077/149‑1 (MQ=255)
                                          gTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCa                                                                                                       <  1:243736/149‑1 (MQ=255)
                                           ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCAc                                                                                                      >  1:240876/1‑149 (MQ=255)
                                           ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCAc                                                                                                      <  4:551625/149‑1 (MQ=255)
                                                  tgctAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAAc                                                                                               <  3:163538/149‑1 (MQ=255)
                                                     tAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACata                                                                                            >  1:324116/1‑149 (MQ=255)
                                                     tAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACata                                                                                            >  4:361246/1‑149 (MQ=255)
                                                     tAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACata                                                                                            >  4:560861/1‑149 (MQ=255)
                                                         gATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTcc                                                                                        >  2:580839/1‑149 (MQ=255)
                                                         gATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTcc                                                                                        >  1:450347/1‑149 (MQ=255)
                                                            gACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCtt                                                                                     <  3:324125/149‑1 (MQ=255)
                                                                   cTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTTCCTTGCTCTgg                                                                              <  1:432666/149‑1 (MQ=255)
                                                                   cTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTgg                                                                              >  4:271959/1‑149 (MQ=255)
                                                                            aCCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCAtt                                                                     >  4:92301/1‑149 (MQ=255)
                                                                             ccAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTc                                                                    >  1:20350/1‑149 (MQ=255)
                                                                              cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACat                                                                                             >  1:525018/1‑123 (MQ=255)
                                                                              cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACat                                                                                             <  2:525018/123‑1 (MQ=255)
                                                                              cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCt                                                                   <  1:580839/149‑1 (MQ=255)
                                                                                      ggACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAc                                                           >  3:526967/1‑149 (MQ=255)
                                                                                      ggACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAc                                                           <  1:563684/149‑1 (MQ=255)
                                                                                        acaacaTGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGActc                                                         <  4:187645/149‑1 (MQ=255)
                                                                                                  tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                               >  3:559181/1‑149 (MQ=255)
                                                                                                  tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                               >  4:238703/1‑149 (MQ=255)
                                                                                                  tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                               >  3:145068/1‑149 (MQ=255)
                                                                                                  tactaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                               >  1:163133/1‑149 (MQ=255)
                                                                                                           ccGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTgctg                                      <  3:295129/149‑1 (MQ=255)
                                                                                                             gCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTgctggc                                    <  3:96093/149‑1 (MQ=255)
                                                                                                                 gTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCg                                >  2:525264/1‑149 (MQ=255)
                                                                                                                     tAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAg                            <  2:324116/149‑1 (MQ=255)
                                                                                                                         aCCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCGCTTGCTCTGGTTCACCACTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGc                        >  4:337636/1‑149 (MQ=255)
                                                                                                                         aCCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGc                        >  1:332872/1‑149 (MQ=255)
                                                                                                                           cTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGa                             <  2:63757/142‑1 (MQ=255)
                                                                                                                           cTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGa                             >  1:63757/1‑142 (MQ=255)
                                                                                                                              tgaagtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCAc                   <  4:145068/149‑1 (MQ=255)
                                                                                                                              tgaagtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCAc                   <  2:163133/149‑1 (MQ=255)
                                                                                                                                 agtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGc                <  2:109698/149‑1 (MQ=255)
                                                                                                                                  aaaaaaaaTGGCGCACAATGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCt               >  3:122557/4‑149 (MQ=255)
                                                                                                                                               cacaTTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGa  <  3:361246/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGA  >  NC_000913/1757541‑1757832

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